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ABSTRACT 



The external nucleate pool-boiling heat-transfer 
coefficient of both smooth and enhanced horizontal tubes 
in R-114-oil mixtures (0 to 10 mass percent oil) were 
measured for heat fluxes from 1 to 91 kW/m 2 at a saturation 
temperature of 2.2 C. The enhanced tubes tested were a 
Wieland GEWA-T finned tube containing 1.02 fins/mm, a 
Hitachi Thermoexcel -E and a Hitachi Ther moexce 1 -HE tube. 
The Thermoexce 1 -E and -HE tubes with their re-entrant 
cavity designs were found to improve the heat-transfer 
coefficient over the smooth tube value at a constant heat 
flux by a factor of approximately 7 in oil-free R-114, 
while the GEWA-T tube improved the coefficient by a factor 
of about 4. While all of the tubes showed a generally 
decreasing performance with the presence of oil, the GEWA-T 
tube resulted in the minimum reduction. This tube showed 
up to 20% and 35% reduction (compared to the oil-free case) 
at 3 percent and 10 percent oil, respectively. The 
Thermoexce 1 -E and Ther moe x ce 1 - HE tubes showed performance 
reductions of up to 40% and 60%, respectively, with 3 
percent oil. At a practical heat flux of about 30 kW/m 2 
with 3 percent oil, the GEWA-T, Ther moexce 1 -E and 
Ther moexce 1 -HE tubes outperformed the smooth tube by 
factors of 4.8, 4.6 and 4.0, respectively. 
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I ._ INTRODUCTION 



A. BACKGROUND 

In the design of refrigeration systems for use aboard 
Naval vessels, the desire is to make the system as 
efficient as possible. To achieve this higher efficiency, 
one must consider the use of a high pressure refrigerant, 
since this gives a higher energy transfer per unit volume 
of vapor (i.e., E y = h* 9 /v g ) [Ref. 13, and results in less 
of a pressure drop along the vapor path for a given heat 
load. R-12, which is a high pressure refrigerant, has an 
energy transfer per unit volume of vapor of 8570 kJ/m 3 
CRef. 23 and, since it is higher pressure, it requires 
heavier gauge components. R-ll, on the other hand, is a 
low pressure refrigerant with a low energy transfer per 
unit volume of vapor of 2086 kJ/m 3 and requires lighter 
gauge material . 

Table 1 provides a short summary of some of the 
available refrigerants and their respective energy transfer 
per unit volume of vapor, their relative toxicity and their 
required system pressure for optimum operation. As can be 
seen from Table 1, R-12 and R-22 possess the highest values 
of energy transfer per unit volume of vapor, and therefore 
would provide systems of higher efficiency than the other 
refrigerants listed since pressure losses in the system 



wou 1 d 



be 



less. 



However , 



since both R-12 and R-22 are 



high-pressure refrigerants, their respective systems would 
be significantly heavier larger than those of the other 
refrigerants. However, R-114 being a moderate-pressure 
refrigerant with an energy transfer per unit voulume of 
vapor of 3800 kJ/m 3 , it would allow for a relatively 
lighter component system than those of R-12 and R-22. 

TABLE 1 [Ref. 3] 

SUMMARY OF REFRIGERANTS 

Opera t i ona 1 
Sy s tern 



Refrigerant 


k J /m 3 


Toxicity 




Pressure 


R - 1 1 


2086 


moderate 


(Group 


5) 


1 o w 


R-12 


8570 


1 ow 


(Group 


6) 


hi gh 


R-22 


12480 


moderate 


(Group 


5) 


high 


R-113 


1128 


hi gh 


(Group 


4) 


1 ow 


R-114 


3800 


moderate 


( Group 


6) 


moderate 


Add i t i ona 1 1 y , 


R-114 belongs to the 


r ef r i geran t 


group (Group 



6) having the lowest toxicity [ReF. 21 , and it has also 
been proven fairly stable with temperature [Ref. 23. 
Therefore, R-114 is a viable alternative to the higher 
pressure refrigerants. 



B. REVIEW OF REFRIGERANT-OIL MIXTURE BOILING BEHAVIOR 

The presence of oil in the evaporators of 
refrigeration systems is a common occurrence. This is due, 
in part, to the use of lubricants in the hermetically 
sealed refrigeration compressors. Over a period of time, a 
considerable amount of oil is introduced into the 
refrigerant, subsequently changing the physical 
properties of the refrigerant. Jensen and Jackman [Ref. 3] 
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reported that the density and the specific heat behave 
linearly in the refrigerant-oil mixtures, however, the 
viscosity and the surface tension do not. Henrici and 
Hesse [Ref. 4] experimentally determined in 1971 that the 
surface tension of mixtures first decreased for R-114-oil 
mixtures of up to 2.5 percent oil, and at higher oil 
concentrations, the surface tension increased. This non- 
linear behavior of the properties of the mixtures makes it 
difficult to explain the changes in the heat- tr ans f er 
coefficient for these mixtures. 

The most evident example of the addition of oil to the 
refrigerant is the presence of foam in the mixture when 
boiling takes place. With the oil concentration above 1 
percent, there is a substantial amount of foam generated 
from the nucleate boiling surface. This foam is generated 
because the refrigerant in the r e f r i ge r an t - o i 1 mixture is 
more volatile than the oil and therefore can vaporize, 
creating a gas bubble surrounded by an oil-rich layer 
[Ref. 5]. Figure 1.1 shows the ideal model of the growth of 
a bubble in a refrigerant-oil mixture. ‘These generated 
bubbles are at lower density than that of the surrounding 
refrigerant-oil mixture. The bubbles ascend to the liquid- 
vapor interface and collect at the interface and produce a 
foam layer. As reported by Henrici and Hesse [Ref. 4], the 
foaming action is most pronounced in r ef r i ge ran t -o i 1 
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Figure 1.1 Oil Concentration Gradient in a Bubble in 

Refrigerant-Oil Mixtures [Ref. 5]. 
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mi xtures of 



1 to 10 percent oil, and provide a the 



presence of oil results in’ a subsequent decrease in the 
heat- transfer performance of the boiling tube. 

Numerous researchers have tried to explain the reasons 
for the decrease in the hea t - t r ans f er coefficient in 
r ef r i ger ant -o i 1 mixtures of increasing oil content. Thome 
[Ref. 63, after an extensive literature search, found that, 
in 1956, Van Wijk et al. explained the effect as being the 
result of the evaporation of the more-volatile component 
leaving an oil-rich layer which has a higher resultant 
local boiling point. This increased local boiling point 
requires an increase in the amount of superheat necessary 
to continue the vaporization and bubble growth. The result 
is a decrease in the hea t - t r ans f e r coefficient. 

Stephan and Preusser [Ref. 7] determined that the work 
required to form bubbles in the refrigerant-oil mixtures 
was greater than in an equivalent amount of pure 
refrigerant. Therefore, the heat- transfer coefficient was 
lower for the refrigerant-oil mixture than for the pure 
refrigerant. Chongr ungr eong and Sauer [Ref. 8] reported 
that the rate of heat diffusion, which is governed by the 
thermal properties of the oil-rich layer, limits the bubble 
growth and that the effects of surface tension are 
negligible. Thome [Ref. 6] agreed with the above factors 
and concluded that the viscosity variations are important 
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in explaining the decrease in the heat-transfer 
coefficient for some refrigerant-oil mixtures. 

Although researchers have not decided on any one set 
of factors that determine the reasons for the decrease in 
the heat- t r ans f er coefficient as oil content increases in 
refrigerant-oil mixtures, all agreed that the physical and 
thermal properties are important factors in explaining the 
heat - tr ans f er behavior of r e f r i ger an t -o i l mixtures. 

C. REVIEW OF BOILING PERFORMANCE FROM SMOOTH TUBES 

The U. S. Navy attempts to reduce the size of its 
shipboard refrigeration systems by the use of alternate 
refrigerants (such as R-114). At the same time, the Navy 
desires to increase the performance of these systems by the 
use of enhanced surfaces. However, as discussed by 
Wann i ar achch i et al. [Ref. 9], only a limited amount of 
data is available regarding the boiling of refrigerants in 
the presence of oil. 

Stephan [Ref. 10] conducted experiments using a 
smooth horizontal plate as a boiling surface in R-12 
and R-22 refrigerants and reported up to 50 percent 
reductions in the boiling hea t - t r an s f e r coefficient 
with a refrigerant-oil mixture of 9 percent oil (by 
mass), while 50 percent oil resulted in a 90 
percent reduction in the boiling heat- transfer coefficient. 
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In 1971, Henrici and Hesse [Ref. 4] conducted 
experiments involving the boiling of R-114 in the presence 
of oil from a smooth copper tube and reported a maximum 
decrease of 20 percent decrease in the boiling heat- 
transfer coefficient with 1 percent oil by mass. In 1981, 
Stephan and Mitrovic [Ref. 11] reported that for the GEWA-T 
tube (manufactured by the Wieland Company) in R-114-oil 
mixtures, the boiling hea t - t rans f e r coefficient of the 
enhanced surface was significantly altered by the presence 
of the oil. In 1985, Reilly [Ref. 5) reported a drop in 
the heat- transfer coefficient of 0 to 35 percent (depending 
on the heat flux) for a R-114-oil mixture of 10 percent oil 
for a smooth copper tube. 

D. REVIEW OF BOILING PERFORMANCE FROM ENHANCED SURFACES 

Arai et al. [Ref. 121 tested a 200-ton R-12 
centrifugal water chiller. The chiller which used an 
enhanced-sur f ace tube, commercially known as Ther moexce 1 -Em 
made by Hitachi, was 28 percent shorter in length and had 
an approximate improvement in the overall heat-transfer 
coefficient of 50 percent. Additional research by Yilmaz 
and Westwater [Ref. 131, Marto and Lepere [Ref. 14], and 
Carnavos [Ref. 15] on commercially available enhanced 
surface tubes, in various refrigerants other than R-114, 
indicated that a porous-coated surface exhibited the best 
boiling heat-transfer coefficient in a pure refrigerant. 
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Comparisons by Reilly [Ref. 5] indicated that in R-114, 
Union Carbide T s "High Flux" tube with a porous -coated 
surface resulted in an improvement of at least a factor of 
seven over a smooth copper tube at the same operational 
parameters. Reilly’s findings complemented the research of 
Yilmaz and Westwater [Ref. 131. Additionally, Reilly 
reported that for a por ous-coated surface in an R-114-oil 
mixtures of 1 percent oil, the oil was seen to degrade the 
nucleate pool-boiling heat- transfer performance by about 
20 percent. 

1. Surfaces with Reentrant Cavities 

Reentrant cavities have been investigated for a 
considerable period of time. Griffith and Wallis [Ref. 15] 
determined that the reentrant cavity geometry is important 
in the generation of nucleate boiling sites in two ways; 
(1) the diameter of the reentrant cavity mouth determines 
the amount of superheat required to initiate boiling, 
and (2) the cavity shape determines the stability once 
the boiling has begun. Webb [Ref. 16] stated that the 
key to the high performance of the reentrant grooved 
structures can be attributed to three factors? (1) 
reentrant cavity within a critical size range, (2) 
interconnected cavities, and (3) nucleation sites of a re- 
entrant shape. If the cavities are interconnected, 
adjacent cavities can activate each other. Webb [Ref. 16] 
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also stated that reentrant cavities provide a stable vapor 
trap, which can remain active at low values of superheat* 

The GEWA-T (1.02 fins/mm) enhanced - sur f ace tube, 
manufactured by W i e 1 and-Werke AG, is a copper tube with 
spiral fins of a T-shaped profile. It has a reentrant- type 
cavity formed by cold working an integral finned tube. The 
two adjacent T-fins form a spiral cavity with an opening 
of 0.25-mm (see Figure 1.2). 

The Thermoexce 1 -E and Thermoexce 1 -HE enhanced - sur f ace 
tubes are made by the Hitachi Company. Like the GEWA-T 
tube, the Thermoexce 1 -E (Figure 1.3) and Thermoexce 1 -HE 
(Figure 1.4) tubes are also formed by cold working integral 
fins; however, these tubes have low fins which have a 
small space cut out at the fin tips, giving the appearance 
of a sawtooth. These sawtooth fins are bent parallel to 
the tube axis, to a horizontal position, forming tunnels 
with evenly spaced surface pores. This cold working 
procedure results in a high area density of reentrant 
nucleation sites. The Thermoexce 1 -HE tube is a variation 
of the Thermoexce 1 -E tube, designed for improved 
performance in 1 ow-heat- f 1 ux regions. 

Experimental data showing the effect of oil 
concentrations on the heat-transfer coefficient of the 
GEWA-T, the Thermoexce 1 -E and Thermoexce 1 -HE tubes in R-114 
are lacking. In an effort to gain more understanding of 
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i i 
i i 
i 1 
i i 




(a) Schematic of GEWA-T Surface Cross Section 




(b) Cutaway o-f the (TRVA-T Surface 




(c) Micrograph of the GEV/A-T Surface 

Figure 1.2 Schematic, Cut-away Cross Section and 

Micrograph of the GEWA-T Surface. 
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Figure 1.3 



Surface Photograph and Enlarged Outer 
Surface of the Thermoexce 1 -E Surface. 
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F i gure 1 . 4 



Sur face 
Surface 



Photograph and Enlarged Outer 
of the Thermoexce 1 -HE Surface. 
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the pool-boiling performance of these tubes, the David W. 
Taylor Naval Ship Research and Development Center sponsored 
this investigation. The smooth copper tube and the GEWA-T 
tube were provided by the Wieland Company. The Thermoexcel- 
E and Ther moexce 1 -HE tubes were provided by the Hitachi 
Company . 

E. THESIS OBJECTIVE 



The objectives of this thesis are: 



1 . 


take 


data on 


i a 


smooth 


tube 


with the 


use of 




auxiliary heaters 


together 


with 


a convection shield 




to 


enable the 


collection of 


more 


data with a 


mini mum 




of operator effort, 










2. 


take 


boiling 


data 


on a smooth 


copper tube 


in R-114 




to 


provide a 


base 1 i ne 


for 


fol low-on comparison 




for 


enhanced 


tubes at oil 


concentrations 


of 0, 1, 




2, 


3, 6 and 


10 


percent 


(by 


mass) at a 


bo i 1 i ng 




temperature of 


2. 2 


C and a 


heat 


-flux range 


from 500 




W/m 2 


to 91 kW/m 2 . 










3. 


take 


data on 


the 


following tubes for the above- 



-mentioned conditions: 



a. GEWA-T (1.02 fins/mm), 

b. Thermoexce 1 -E, and 

c. Thermoexce 1 -HE . 
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II. DESCRIPTION OF EXPERIMENTAL APPARATUS 



A. OVERALL APPARATUS 

The equipment used for this study consisted of nine 
basic components: (a) a Pyrex-glass tee for the pool 
boiling of the R-114 liquid; Cb) a Pyrex-glass tee for the 
condensing of the R-114 vapor; Cc) an R-114 liquid 
reservoir; (d) a water-ethylene glycol mixture sump; (e) an 
R-12 refrigeration system with an external proportioning 
valve; (f) oil reservoir; (g) a vacuum pump; 
(h) instrumentation and da ta-acqu i s i t i on system; and (i) a 
convection shield. Figure 2.1 shows schematically the 
arrangement of the various components. Complete details of 
the design, construction and operation of the apparatus 
are provided by Karasabun [Ref. 173, Reilly [Ref. 5] 
and Wanniarchchi et al. [Ref. 13. The original design of 
Karasabun was modified by the inclusion of the convection 
shield in the pool-boiling section as shown in Figures 2.2 
and 2.3. 

The operation of the apparatus consisted of the 
boiling of R-114 liquid in the lower glass tee (1), 
and condensing of the R-114 vapor in the upper 
glass tee (2). The resultant condensate was fed to 
the distribution tube within the boiling section by 
gravity. A wa t er -e thy 1 ene glycol mixture from a 30-gallon 
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Schematic of Experimental Apparatus. 



F i gure 2. 1 
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Pyrex Boiling Tee 




Figure 2.2 Schema t i c of Convection Shield in 

Pyrex Boiling Tee. 
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(b) Side View 



Fi gure 2. 3 



Photograph of Convection Shield. 
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sump was used to cool the copper condensing coil, thereby 
condensing the R-114 vapor. This water-ethylene glycol 
mixture sump was maintained between -18 and -14 C by a 
one-half ton R-12 air-conditioning system. The mixture was 
pumped through the condensing coil and controlled manually 
with valve VC . 

Valve VC controls the amount of the water-ethylene 
glycol mixture to the condensing coil, thereby controlling 
the pressure within the test loop. For example, the opening 
of valve VC allows more of the wa ter - e t hy 1 ene glycol 
mixture to pass through the condensing coil, resulting in 
additional condensing of the R-114, and subsequently 
lowering the temperature and pressure within the entire 
sys tern . 

As the subcooled R-114 condensate enters the boiling 
section of the apparatus through the distribution tube, it 
strikes the convection shield. This shield provides 
increased mixing of the subcooled R-114 condensate with the 
remaining R-114 (or R-114/oil mixture) in the lower portion 
of the boiling section. The shield deflects the incoming 
condensate away from the boiling surface and allows for 
further mixing of the condensate and the existing R-114 
liquid by means of the auxiliary cartridge heaters. This 
deflection of the incoming condensate allowed for a more 
even temperature distribution of the liquid R-114 and 
therefore decreased the previously reported temperature 
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gradients within the liquid region. This additional mixing 
of the liquid R-114, coupled with the ability to use 
auxiliary cartridge heaters in the liquid and allowed for 
the timely achievement of steady-state conditions without 
the previously observed thermal gradients in the liquid 
region of the tests done without the aid of the convection 
shield. It was also observed throughout the operation of 
the apparatus that the data obtained using the convection 
shield had a higher degree of repeatability than the data 
obtained without the shield. 

B . BO I L I NG-TUBE CONSTRUCT I ON 

The tubes tested during the course of this study were 
the smooth copper tube, the GEWA-T (1.02 fins/mm), the 
Thermoexce 1 -E and the Thermoexce 1 -HE tubes. All of these 
tubes had an active test section (i.e., heated) of length 
of 203.2-mm (8 in), with an outside diameter of 15.9 mm 
(5/8 in) and an inside diameter of 12.7 mm (1/2 in). The 
remaining 114.3-mm section on either side of the heated 
test section was smooth and unheated and did not nucleate 
under any heat flux or any oil concentration condition. 
Figure 1.2 shows the schematic of the GEWA-T finned 
surface. Karasabun [Ref. 173 described the analysis used 
in the da ta- reduct i on program, allowing for the two end 
surfaces to be treated as extended fins from the center 
section and subsequently accounting for the heat losses. 
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The center section of each tube tested was heated by a 



1 000- W 240 -V stainless steel cartridge heater inserted in 
the center section of each tube. The cartridge heater was 
inserted in the boiling tube and surrounded by a copper 
sleeve with eight 1.3 mm by 1.3 mm thermocouple channels. 
Figure 2.4 shows the layout of the installed tube in the 
boiling section. Figure 2.5 shows the details of the 
thermocouple layout. The thermocouple hot junctions were 
inserted into the copper sleeve and the copper sleeve was 
then tinned; and heat was maintained to keep the solder in 
a molten state as the sleeve was inserted and positioned in 
the center section of the boiling tube. Upon completion of 
the thermocouple installation, calibration of the 
thermocouples was 'achieved. Reilly [Ref. 5] determined in 
the course of his investigation that since the data- 
reduction program utilized the differences between the 
thermocouples in all computations, such as wall 
temperatures minus saturation temperatures, the corrections 
provided by the thermocouple calibration were basically 
unnecessary since the calibration procedures were only 
necessary for items dependent on absolute temperature. 
Appendix A summarizes the procedure used by Karasabun 
[Ref. 17] and Reilly [Ref. 5] to calibrate the 
thermocoup 1 es . 
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Figure 2.4 Sectional View of Boiling Tube 
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(a) View of the boiling tube thermocouple locations as seen 
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(c) End views of the boiling tube. 



Figure 2.5 Boiling Tube Thermocouples Channels. 
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C. CONVECTION SHIELD 

Initial tests on the smooth copper tube revealed an 
axial thermal gradient in the liquid over the length of the 
boiling tube. This thermal gradient increased with 
increasing oil concentration and with increasing heat flux, 
with the thermal gradient ranging from 0.5 K to 1.5 K over 
the length of the boiling section. This thermal gradient 
resulted in questionable data since steady -state conditions 
could not be maintained over the course of the individual 
data runs. 

During the early stages of this investigation, it was 
found that the steady-state conditions were fairly 
difficult to achieve. This was especially true after making 
a step change in the heat flux. For example, if the heat 
flux was increased, the cooling flow rate to the condenser 
must be increased to maintain the system pressure. Since 
the system showed about a 30-second time lag between the 
change of the valve position and the subsequent system 
response, the pressure was often seen to undergo cyclic 
variations. With considerable experience of the operator, 
it was possible to achieve steady-state conditions within 5 
to 10 minutes, as demonstrated by Reilly [Ref. 53. 

In order to minimize the efforts necessary by the 
operator, two 600-W auxiliary cartridge heaters were 
installed within the boiling tee. With this arrangement, 
steady-state conditions could be maintained very easily by 
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maintaining the total heat load in the boiling tee at a 
constant value. For example, if the heat duty for the 



bo i 1 i ng 


tube was increased, the heat duty on the auxiliary 


heaters 


must be decreased by the same magnitude. In this 


manner , 


steady-state conditions could be maintained with 


mini mum 


attention. However, this arrangement required that 



the heat duty on the auxiliary heaters be a maximum when 
the heat duty on the boiling tube is at a minimum. As 
determined by the early data runs (see Figure 2.6), the 
convective effects created by the bubbles of the auxiliary 
heaters created an artificially enhanced hea t - t r ans f er 
coefficient at low heat-flux settings. 

In order to alleviate the problem discussed 



above , 


while maintaining the advantages offered by 


the 


auxiliary heaters, the present investigator 



considered the use of a convection shield. As 
shown in Figure 2.2, this shield was positioned 
around the test tube, thus isolating it from the 
convective bubble effects created by the auxiliary heaters. 

The convection shield also provided increased 
mixing of the refrigerant and condensate without 
interfering with the boiling phenomenon of the boiling 
tube or the circulation of the oi 1 -refrigerant 
mixture in the boiling section by both the boiling 
tube and the auxiliary heaters. 
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Figure 2.6 Convective effects of Auxiliary Heaters 
with and without Convective Shield 



D. 



DATA ACQUISITION AND REDUCTION 



The dat a-acqu i s i t i on system used for this research was 
a Hewlett-Packard 9826A computer and a Hewlett-Packard 
3497A Automatic Data Acquisition Unit. The HP-9826A 
computer was used to control the HP-3497A, to provide hard- 
copy storage capacity and to analyze the data obtained 
during this investigation. 

The da ta-acqu i s i t i on system provided the means by 
which the researcher was able to review system performance 
and stability prior to the taking of any data. System 
parameters were monitored on the video screen following the 
scanning of all of the channels by the HP-3497A. Following 
the automatic scanning of all of the channels, the computer 
would display the desired saturation temperature (2.2 °C) , 
the observed saturation temperature, the liquid temperature 
at both ends of the boiling section and their average, the 
vapor temperature, and the temperature of the water- 
ethylene glycol mixture sump (-15 to -18 °C). Once the 
steady-state conditions were met for the required satura- 
tion temperature within the prescribed variance (±0.1 K) 
of the desired saturation temperature, the da ta-acqu i s i t i on 
system was initiated for the individual data run. 
Since there was a time lag in the request for the 
taking of data and the actual printing of the data due to 
the sampling of each channel twenty times before recording 
the data, the da ta-acqu i s i t i on system would print the data, 
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allowing the researcher to review the results prior to 
storing the data. If, during the sampling of the data, a 
transient trend was observed due to the cycling of the 
R-12 refrigerant system or the temperature difference 
between the thermocouples measuring the liquid R-114 was 
not within the specified temperature range, the researcher 
had the capability to reject the affected data set. This 
capability allowed the researcher to apply a cursory 
analysis of each data point and ascertain if the 
temperatures measured during that individual data set were 
within the prescribed limits of variance. 

The information scanned on each channel of the HP- 
3497A were thermal emf r s provided by coppe r -cons tan tan 
(type-T) thermocouples made of 0 . 25 -mm - d i ame t e r (30 gauge) 
wire. Two additional measurements were provided to the HP- 
3497A by the po we r - s en s i n g device, as described by 
Karasabun [Ref. 173, which converted the AC voltage dialed 
into the variac of the main boiling tube heater into a DC 
signal, allowing the HP-3497A to accurately scan the 
converted value and then allow the HP-9826A to 
automatically apply this value to the analysis calculation. 

Upon the initiation of the command to take data, the 
data-acqu i s i t i on system would sample each channel and then 
compute, according to the step-wise procedure outlined by 
Karasabun [Ref. 173 and print the results. Appendix B 
contains a summary of the da ta-acqu i s i t i on program and a 
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listing of the data-reduction program utilized during the 



course of this research. Appendix C contains 
representative sample of the data printout received 
each data run. 
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111. EXPERIMENTAL PROCEDURES 



A. INSTALLATION OF TUBE IN BOILING APPARATUS 

The external surfaces of each of the boiling tubes 
were cleaned with Nitol (2 percent nitric acid and 98 
percent ethyl alcohol) solution and then rinsed with 
acetone to ensure tube cleanliness. This treatment was also 
effective in the removal of oil and grease residues from 
the tube surfaces. 

After air drying, the individual tube to be tested was 
installed in the boiling section of the apparatus. The 
system was sealed and then evacuated to approximately 29 
inHg by means of a portable mechanical vacuum pump (6), 
shown in Figure 2.1. The apparatus remained under the 
evacuated condition for a period of about one hour to check 
if there were any major leaks. System pressure was measured 
by means of a Marsh pressure gauge (30 inHg to 150 psi 
range with 12.5 inHg and ±0.5 psi accuracy, 
respectively). If there was not a noticeable drop in the 
vacuum reading as observed on the pressure gauge, the 
system was charged with R-114 vapor, which was condensed by 
the wa ter -ethy 1 ene glycol mixture and then the absolute 
pressure was raised to approximately 0.186 MPa (27 psi). 
Upon reaching this pressure, an Automatic Halogen Leak 
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Detector, T1F 5000, with a sensitivity of 3 ppm 
concentration, was used to check for system leaks. 

Upon successful completion of the leak-detection tests 
and taking the necessary steps to fix any leaks found, the 
R-114 liquid level was adjusted to be 20 mm above the top 
of the boiling tube surface. Liquid-level adjustments were 
made by either boiling off the excess refrigerant to the 
R-114 reservoir (3), or drawing the required amount of 
refrigerant liquid into the system from the R-114 
reservoir. The required liquid level for boiling was 
maintained constant for all tubes tested by means of a 
permanent scribe mark on the inside of the oil intake 
endbell of the boiling section. At this initial refrigerant 
level, the system was now ready for boiling tests in pure 
refrigerant. 

B. GENERAL SYSTEM OPERATION 

Table 2 provides a summary of all of the 83 data runs 
made during this investigation. Each data run is annotated 
to indicate the conditions involved in the respective run. 
The data runs were numbered sequentially with an indicator 
code for each tube preceding the number of the data run. 
Since this research is a follow up of previous research, 
the initial data run was indicated as number 168. All data 
runs were undertaken at a saturation temperature of 2.2°C. 
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A listing of the boiling tube indicator codes for Table 2 
is as follows: 

WH -- Wieland Hard Copper Smooth Tube 
GTB -- Wieland Gewa-T (1.02 fins per mm) 

TXE -- Hitachi Thermoexce 1 -E 
THE -- Hitachi Thermoexce 1 -HE. 

Initially, all data runs were taken over a range of 
nine different heat-flux settings, ranging from 2 kW/m 2 to 
91 kW/m 2 . Additional heat fluxes of 0.5, 0.8, 1.2 kW/m 2 

were later included for all tubes with enhanced surfaces to 
study the onset of nucleate boiling in the o i 1 -r ef r i ger ant 
mixtures and to study the effects of hysteresis on these 
tubes. The initiation of all data runs was carried out in 
the same manner to ensure consistency in the obtained data. 
Each run was started by an initial system "warm-up” period 
allowing the R-12 refrigeration system to obtain a minimum 
coolant starting temperature of about -15°C. This was 

accomplished by lighting off of the R-12 system in 
conjunction with the effluent recirculation/cooling pump 
(8) and slight opening of the control valve VC. This 
procedure allowed the system to gradually decrease system 
temperature and pressure, while minimizing the adverse 
effects of • shocking the boiling tube. The system was 
operated in this condition until a stable condition was 
established with the saturation temperature of 2.2°C. This 
operation was most critical when the testing sequence 
involved increasing heat-flux conditions since 
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TABLE 2 

Summary of Data Runs 



Run 

No. 


% 

0i 1 


Heat 

Flux 


No .of Data 
Points 


Aux i 1 i ar y 
Heater 


Convect i on 
Shield 


WH168 


0 


DEC 


20 


off 


without 


WH 169 


0 


INC 


22 


off 


w i t hou t 


WH170 


1 


DEC 


20 


off 


w i thout 


WH 171 


1 


INC 


22 


off 


without 


WH172 


2 


DEC 


20 


off 


without 


WH 173 


2 


INC 


20 


off 


w i thout 


WH 1 74 


3 


DEC 


20 


off 


w i thout 


WH175 


3 


INC 


22 


off 


wi thout 


WH 176 


3 


DEC 


20 


off 


w i thout 


WH177 


6 


DEC 


20 


off 


w i thout 


WH178 


6 


INC 


21 


off 


w i thout 


WH 179 


10 


DEC 


18 


off 


w i thout 


WH 1 80 


10 


I NC 


13 


off 


w i thout 


WH 18 1 


10 


INC 


12 


off 


w i thout 


GTB182 


0 


DEC 


24 


on 


with 


GTB183 


0 


I NC 


23 


on 


with 


GTB184 


0 


DEC 


22 


on 


with 


GTB185 


1 


DEC 


20 


on 


with 


GTB186 


1 


INC 


22 


on 


with 


GTB187 


2 


DEC 


34 


on 


with 


GTB188 


2 


I NC 


18 


on 


-with 


GTB189 


3 


I NC 


22 


on 


with 


GTB190 


3 


DEC 


20 


on 


with 


GTB191 


6 


DEC 


27 


on 


with 


GTB192 


6 


I NC 


20 


on 


with 


GTB193 


10 


INC 


20 


on 


with 


GTB194 


10 


DEC 


21 


on 


with 


WH 195 


0 


INC 


22 


on 


with 


WH 196 


0 


DEC 


20 


on 


with 


WH 197 


0 


DEC 


17 


on 


with 


WH 1 98 


0 


INC 


22 


on 


with 


WH199 


1 


I NC 


20 


on 


with 


WH200 


1 


DEC 


20 


on 


with 


WH201 


2 


INC 


20 


on 


with 


WH202 


2 


DEC 


20 


on 


with 


WH203 


3 


I NC 


20 


on 


with 


WH204 


3 


INC 


20 


on 


with 


WH205 


3 


DEC 


20 


on 


with 


WH206 


3 


I NC 


21 


on 


with 


WH207 


6 


DEC 


20 


on 


with 


WH208 


6 


INC 


20 


on 


with 


WH209 


10 


INC 


20 


on 


with 
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TABLE 2 

Summary of Data Runs (cont'd) 



Run 
No . 


% 

Oi 1 


Heat 

Flux 


No. of Data 
Points 


Aux i 1 i ary 
Heater 


Convect 
Sh i e 1 


WH210 


10 


DEC 


21 


on 


with 


TXE210 


0 


DEC 


18 


on 


with 


TXE211 


0 


DEC 


12 


on 


with 


TXE212 


0 - 


DEC 


12 


on 


with 


TXE213 


0 


DEC 


14 


on 


with 


TXE214 


0 


DEC 


12 


on 


with 


TXE215 


0 


I NC 


27 


on 


with 


TXE216 


0 


DEC 


14 


on 


with 


TXE217 


1 


INC 


27 


on 


with 


TXE218 


1 


DEC 


14 


on 


with 


TXE219 


2 


DEC 


16 


on 


with 


TXE220 


2 


DEC 


16 


on 


with 


TXE221 


2 


I NC 


25 


on 


with 


TXE222 


3 


DEC 


15 


on 


with 


TXE223 


3 


I NC 


22 


on 


with 


TXE224 


3 


I NC 


24 


on 


with 


TXE225 


6 


I NC 


30 


on 


with 


TXE226 


6 


DEC 


16 


on 


with 


TXE227 


6 


INC 


25 


on 


with 


TXE228 


10 


INC 


24 


on 


with 


TXE229 


10 


DEC 


25 


on 


with 


TXE230 


10 


INC 


20 


on 


with 


THE231 


0 


INC 


20 


on 


with 


THE232 


0 


DEC 


12 


on 


with 


THE233 


1 


DEC 


14 


on 


with 


THE234 


1 


INC 


24 


on 


with 


THE235 


1 


INC 


27 


on 


with 


THE236 


2 


DEC 


14 


on 


with 


THE237 


2 


INC 


12 


on 


with 


THE238 


2 


INC 


8 


on 


with 


THE240 


3 


DEC 


14 


on 


with 


THE241 


3 


INC 


18 


on 


with 


THE242 


3 


INC 


8 


on 


with 


THE243 


3 


I NC 


5 


on 


with 


THE244 


3 


I NC 


5 


on 


with 


THE245 


6 


DEC 


14 


on 


with 


THE246 


6 


I NC 


24 


on 


with 


THE247 


6 


I NC 


18 


on 


with 


THE248 


10 


DEC 


14 


on 


with 


THE249 


10 


I NC 


24 


on 


with 
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failure 



to strictly adhere to this would result in 



premature nucleate boiling and the subsequent purging of 
the data. This initial warm-up period required a minimum 
of one hour to ensure adequate system performance. 

Upon reaching stable conditions, the initial heat flux 
for the boiling tube was set by means of the resistor 
variac. Also, the heat flux for the auxiliary cartridge 
heaters was set by a second variac. The heat flux of the 
auxiliary cartridge heaters was limited to settings that 
allowed the heating of the liquid R-114 without creating 
excess turbulence on the boiling tube. As discussed 
earlier, the boiling tube was shielded from the effects of 
the auxiliary heater by means of the convection shield. As 
the oil concentration of the r ef r i gerant-o i i mixture 
increased, the heat-flux setting for the auxiliary heater 
was additionally limited to lower heat fluxes. This 
precluded the generated oil foam from carrying over the 
sides of the convection shield and interfering with the 
operation of the boiling tube. 

After setting of the desired heat flux for the data 
run, the system was monitored for the vapor temperature, 
sump temperature and the liquid refrigerant temperature 
(two readings; left side and right side of the boiling 
tube). Control valve VC was adjusted to allow the system to 
reach steady-state conditions. Additionally, the liquid 
refrigerant temperatures were compared to ensure minimum 
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var i ance . 



After fulfilling the requirements of the steady- 



state conditions, the data-acqu i s i t i on system was initiated 
to take data. 

The data-acqui s i t i on system would scan each assigned 
channel, compute the heat- transfer coefficient and print 
the results for operator viewing. If the data were within 
the prescribed limits, they were stored on a floppy disk 
for permanent records. A representative sample of a data 
printout is contained in Appendix B. 

After each run was completed, oil was added to 
the boiling section by means of a graduated cylinder 



through va 1 ve 


VI. Table 


3 contains the 


1 i s t i ng 


of 


the required 


vo 1 ume 


of 


o i 1 


added to 


obtain 


the 


proper mass 


percent 


of 


o i 1 


for each series of 


data 


runs. Mixing 


of the 


o i 1 


and 


refrigerant 


1 i qu i d 


was 


accomp 1 i shed 


by means 


of 


the 


aux i 1 i ary 


heater , 


to 


ensure a homogeneous 


mi xture. 


Foaming of the 


oil- 


refrigerant mixtures was 


observed 


over the 


entire heat- 


flux range 


of each 


data 


run , 


a 1 though 


carryover of 



the foam over the convection shield was limited to the 
highest heat flux. As the oil concentration increased. 



carryover 




increased, but was 


1 i mi ted 


to 


the 


vertical 


section 


of 


the apparatus. 


Carryover 


o i 1 


was not 


observed 


i n 


the condenser 


section 


of 


the 


apparatus . 



Figures 3.1 to 3.3 demonstrate the foaming action of 
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the 


oil-refrigerant mixtures 


at a 


heat flux of 37 


kW/m 2 


for 


o i 1 


concentrations of 0, 


6 and 10 


percent . 




Data 


runs were made for each 


boiling 


tube tested with 


pure 


R-l 14 


and foil owed by 


runs of 


oi 1 -ref r i gerant 


mixtures 


of 


increasing oil content. 


Oi 1 -refrigerant 


mixtures 


of 


1, 2, 3, 6 and 10 


percent 


oil by mass were 


used 


for a 1 


1 tubes. 












TABLE 3 












OIL CONCENTRATIONS 


VALUES 






Percent 


Total Volume (cc) Step 


Change (cc) 




Oi 1 




of Oil 


i n 


Vo 1 ume of Oil 




0 

1 




26.8 




26.8 




2 




54.2 




27. 4 




3 




82.2 




28.0 




6 




169. 6 




87.4 




10 




295.2 




125.6 
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91 kW/m" 



Figure 3.1 



Foaming Action 



for Pure 



R- 114 
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91 kW/m 2 



Figure 3.2 



Foaming Action for R-114 with 
6 Percent Oil 
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37 kW/m 2 




91 kW/m 2 



Figure 3.3 



Foaming Action for R-114 with 
10 Percent Oil 
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I_v._ RESULTS AND DISCUSSIO N 

A. EFFECT OF CONVECTION SHIELD 

Two sets of data runs were taken on the smooth copper 
tube; one set without the convection shield to validate the 
performance in comparison with the data obtained by Reilly 
(at a saturation temperature of -2.2°C) [Ref. 5J; and a 
second data run using the convection shield and the 
auxiliary cartridge heaters. These data runs allowed for 
an accurate comparison of the data taken over a period of 
time and provided justification for the use of the 
convection shield. As can be seen from Figures 4.1, 4.2 
and 4.3, for 0, 3 and 10 percent oil concentration, 
respectively, the present data with and without the shield 
in place are essentially the same. However, as discussed 
in section II. C, the use of the shield provided consistent 
steady-state conditions over the entire range of heat 
fluxes and oil concentrations, thereby minimizing the 
efforts required by the operator. 

Figure 4.1, which contains the data run for 
pure refrigerant for the smooth copper tube at a 
saturation temperature of 2.2°C, with and without the 
shield, demonstrates the repeatability of the data. 
Also plotted on this figure are Reilly’s smooth-tube 
data for pure refrigerant at a saturation temperature of 
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ure 4.1 Performance Comparison for Pure R-114 Boiling 
from a Smootli Tube, 




(M 

□ 



O 



□ 

O 





C 2 w/M> b 



52 



Figure 4.2 Performance Comparison for R-114 with 3 Percent 
Oil Boiling from a Smooth Tube. 
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Figure 4.3 Performance Comparison for R-114 with 10 Percent 
Oil Boiling from a Smooth Tube. 



-2.2°C. This data set represents a series of data 
runs, starting with a heat flux of 100 kW/m 2 and 
decreasing the heat flux to a minimum of 1 kW/m 2 . The 
three sets of curves are nearly identical with the 
deviation between the Reilly data and the most recent data 
being attributed to the pressure effect and the possible 
aging of the tinned interface between the thermocouple 
sleeve and the inner surface of the smooth copper tube. 
The maximum discrepancy of the computed wall superheat 
between Reilly’s data and the present data is about twenty 
percent at the maximum heat-flux setting. This figure also 
shows the typical uncertainty of about 10 percent computed 
for the wall superheat values during this investigation. 

B. BOILING PERFORMANCE OF THE SMOOTH TUBE 

Figure 4.4 shows the performance of the smooth copper 
tube in pure R-114 and refrigerant-oil mixtures of 3 and 
10 mass percent oil. The behavior of the smooth copper 
tube in pure R-114 liquid demonstrates typical 
nucleate pool-boiling performance. As can be seen in 
Figure 4.4, the line between points A and B represents 
a region of constant slope (on a log-log scale) for 
the variation of heat flux with the tube wall superheat 
(Two-Tsat). This section of the curve represents natural 
convection (i.e., no bubbles were generated). 
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Upon reaching point B (see Figure 4.4). the boiling 
tube wall superheat begins to decrease even though the 
heat flux is continually increased. At this point, the 
tube begins to demonstrate the characteristics of 
incipient nucleate boiling. This region is known as the 
mixed boiling region, where the transition from natural 
convection to nucleate boiling occurs. This event is 
characterized by the increasing number of active 
nucleation sites with increasing heat flux. It should be 
noted that during all data runs, the unheated ends of the 
tubes did not show any nucleation sites. Reilly [Ref. 53 
also stated that the unenhanced ends of the boiling tube, 
at all observed heat fluxes, underwent only natural 
convection despite a small amount of heat axially conducted 
along the tube walls. It was observed that this transition 
from natural convection to nucleate boiling occurs very 
quickly, usually within a few seconds after the activation 
of the first nucleation site. 

At the heat flux represented by point C in Figure 
4.4, all the available activation sites of the boiling tube 
appear to be active and the wall superheat begins to 
increase with increasing heat flux. From point C to 
point D, where the heat flux is increased to a maximum 
of 91 kW/m 2 , new activation sites are generated (maximum 
heat flux was limited by cartridge heater rated output). 
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Figure 4.4 Heat-Transfer Performance for Smooth Tube in 
0, 3 and |0 Percent Oil . 




As the oil concentration was increased, there was a 
significant increase in the amount of foam generated. 

After obtaining a maximum heat flux at point D, the 
heat flux was decreased and is represented as the region 
between points D and E on Figure 4.4. This region of the 
curve represents decreasing wall superheat; however, more 
activation sites that were created at the higher heat-flux 
settings remain active to lower heat fluxes. This 
increased range over which the boiling tube continues to 
nucleate results in an increase in the hea t- t rans f er 
performance since the wall superheat temperature for the 
nucleate boiling at any particular heat flux is less than 
the tube wall superheat for the natural convection heating 
at the same heat flux. 

Stephan [Ref. 9] reported that the effects of adding 
oil introduces a mass-diffusion resistance to the liquid 
and therefore lowers the heat- transfer coefficient. Reilly 
[Ref. 5], in his discussion of his smooth-tube data at a 
saturation temperature of -2.2°C (Figure 4.5), stated that 
the 3 mass percent oi 1 and 10 mass percent oi 1 
concentrations demonstrated higher wall superheat 
temperature than that observed in the pure R-114 liquid. 
However, Reilly observed that the wall superheat for the 3 
percent oil concentration was higher than that for the 
10 percent oil concentration. He attributed this to the 
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Figure 4.5 Reilly’s tRef. 5] Smooth Tube Data for 
0, 3 and 10 Percent Oil. 



non-iinear physical property characteristics of the 



ref r i gerant- 


oil mixtures 


As can 


be 


seen 


i n 


Figure 


4.4, 


the general 


trend 


o f 


the data 


i s 


the 


same 


with 


t he 


excep t i on 


of the 


3 


percent 


o i 1 


concent rat i on 


wa 1 1 


superheat , 


wh i ch near 1 y 


over l ays 


the 


data 


for 


pure R 


-114 



liquid. This deviation from the Reilly data may be due in 
part to the difference in the saturation temperatures. The 



add i t ion of oil 


to the 


pure R-114 


1 i qu i d , 


up to a 


max i mum 


of 10 percent, 


shou 1 d 


delay the transition 


f r om 


na tura 1 


convection to nucleate 


poo 1 boiling 


on the 


t ube . 


With a 


re f r i gerant - o i 1 


mixture 


of 10 mass 


percent 


oil. 


it was 



observed that the surface nucleation sites did appear to 
spread more slowly as the heat flux was increased. This 
condition continued until the tube was fully engulfed in 
the nucleate boiling mode. This observation was also noted 
by Reilly [Ref. 5] and it agrees with the contention of 
Chongrungr eong and Sauer [Ref. 8] and Thome [Ref. 6] that, 
with the exception of the surface tension, the non-linear 
variations in the physical properties of the refrigerant- 
oil mixtures explain the variation in the boiling heat 
transfer of the tube. 

Figure 4.6 shows the hea t - t r ans f er coefficient of the 
smooth tube in r e f r i ge r an t - o i 1 mixtures, plotted as a 
function of heat flux. Figure 4.7 shows the same plot for 
Reilly r s data for a saturation temperature of -2.2°C 
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Figure 4.6 Boiling Heat-Transfer Coefficient for Smooth Tube. 
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Figure A. 7 Boiling Heat-Transfer Coefficient for 
Reilly's Smooth Tube Data [Ref. 53 



[Ref . 



As can be seen, the two graphs are nearly 



5 1 . 



identical, representing good repeatability of the data. 
Also, the effect of the addition of up to 6 percent oil 
(Figure 4.7) consistently represents an approximate 
decrease of 10 percent in the heat- transfer coefficient. 
As the oil concentration was increased from 6 to 10 
percent, the decrease in the hea t - t r ans f e r coefficient was 



a maximum of 


about 35 percent. These decreases in the 


heat- t ransf er 


coefficients were dependent on the comparison 


of the pure 


R-114 and the 10 percent oil concentration 


at equal heat- 


-flux settings. 



C. BOILING PERFORMANCE OF THE GEWA-T TUBE 

Figure 4.8 shows the nucleate pool-boiling performance 
of the GEWA-T (1.02 fins/mm) tube in R-114-oil mixtures. 



The magnitude 


of the tube-wall superheat was considerably 


less than the 


values obtained for the smooth tube discussed 


in Section IV. 


,B. This lower wall superheat can be 


attributed to 


the specially formed fins of the GEWA-T tube. 



The performance of the GEWA-T tube is similar to that of 
the smooth tube with the curves of the increasing heat flux 
to the point of incipient nucleate boiling (point B), 
paralleling the performance of the smooth copper tube 
(Figure 4.4). The same condition exists throughout the 
entire range of decreasing heat fluxes. 
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The region bounded by points A and B shows the 
characteristic natural -convection heating behavior. The 
onset of incipient nucleate boiling at point B is clearly 
defined for the pure R-114 liquid. However, as the oil 
concentration was increased, the point of incipient 
nucleate boiling becomes less clearly defined and the 
region of mixed boiling (point B to C) increased with 
increased heat flux. It was also observed that in the 
region from A to B, the 3 percent oil mixture showed the 
smallest wall superheat and the pure refrigerant showed the 
highest wall superheat. The reason for this behavior is 
not fully known at present, but may be due to convection 
effects within the channels of the GEWA-T tube which are 
enhanced by the foaming action in the presence of oil. 

The region bounded by the points C and D for the 
increasing heat-flux condition is indicative of fully 
developed nucleate pool boiling. The region between points 
D and E represents nucleate pool boiling for decreasing 
heat flux. Figure 4.8 shows that, as the oil 
concentration was increased, the wall superheat increased 
at a given heat flux. This is contrary to the results for 
the condition of increasing heat flux in the region of 
natural -convection heating. 

When comparing the performance of the smooth tube 
(Figure 4.4) to that of the GEWA-Tube (Figure 4.8), the 
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Figure 4.8 Meat-Transfer Performance for GEWA-T Surface. 



maximum adverse effect of the increased oil content for the 



smooth tube, occurs at the higher heat-flux settings. For 
the GEWA-T tube, the maximum decrease in the tube 
performance for increased oil content, occurs in the lower 
and intermediate heat-flux range. 

Figure 4.9 shows the he a t - t r ans f e r coefficient as a 
function of heat flux over the range of R-114-oil 
mixtures (0, 1, 2, 3, 6 and 10 percent oil). Data observed 
over the range of R-114-oil mixtures of 1, 2, 3 and 6 
percent oil are nearly coincident. The 10 percent oil 
concentration data represent an approximate 30 percent 
reduction of the hea t - t r ans f e r coefficient when compared to 
the hea t- t r ans f er coefficient obtained during the data for 
pure R-114. Figure 4.10 shows the direct comparison of the 
GEWA-T tube with the smooth tube for R-114-oil mixtures of 
0, 3 and 10 percent oil. The GEWA-T tube represents an 
overall enhancement of 4.2 over the baseline smooth copper 
tube at a heat flux of 40 kW/m 2 with pure R-114 and an 
overall enhancement of 3.8 with 10 percent oil. These 
results compare favorably with the results of Yilmaz, 
Hwalek and Westwater [Ref. 18] for a 12.3-mm 0D GEWA-T tube 
in p-xylene, in which they obtained an enhancement of 5.3 
in the pool-boiling hea t - t r ans f e r coefficient over their 
baseline smooth copper tube. Additionally, Yilmaz and 
Westwater [Ref. 12] tested a similar surface in isopropyl 
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alcohol in which an enhancement of 2.0 was observed. Marto 



and Lepere [Ref. 13] tested a 17.9-mm OD GEWA-T surface in 
both R-113 and FC-72 and reported enhancements of 2.8 and 
2.5, respectively. 

Figure 4.11 shows more clearly the degradation 
that oil causes in the boiling hea t - t r ans f er 
performance for the GEWA-T tube. This figure plots 
the heat- transfer coefficient of the GEWA-T tube in 
the oil mixture relative to the hea t - t r ans f e r 
coefficient in pure R-114 as a function of oil 
concentration. From the shape of the curves, the 
oil-caused degradation does not show a clear trend 
with the oil concentration at all practical heat 



f luxes. 


Over 


the 


range 


o f 


o i 1 


concentrations, except 


at a 


heat 


f 1 ux 


setting 


o f 


37 


kW/ra 2 , there is only 



a small degradation (up to a maximum of 35%) in the 
boiling performance. This is probably due to the 
large channels and reentrant cavities in the tube, 
al lowing the "pumping” action of the bubbles to continue 
to remove the oil-rich liquid throughout the entire 
heat-flux range without a significant reduction in 
performance. The only significant decrease in 
performance occurred at heat fluxes greater than 5 
kW/m 2 with 10 percent oil. 
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Figure 4.9 Boiling Heat-Transfer Coefficient for GEWA-T Surface. 
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Figure 4.10 Boiling Heat Transfer Coefficient Comparison for GEWA-T 
Surface and Smooth Surface . 
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Figure 4.11 Relative Effect of Oil on GEWA-T Boiling Heat-Transfer 



D. BOILING PERFORMANCE OF THE THERMOEXCEL - E TUBE 

Figure 4.12 shows the nucleate p.oo 1 - bo i 1 i ng 
performance of the Thermoexcel -E tube in R-114-oil 
mixtures. As can be seen in the figure, the region of 
natural-convection heating for the condition of increasing 
heat flux, points A to B, for the pure R-114 is very small, 
with the point of incipient nucleate boiling occurring at a 
reasonably low heat flux (1.5 kW/m 2 ) . At this point of 
incipient- nucleate boiling, the tube performs in the mixed 
boiling region as a limited number of nucleation sites 
become active. From point C to C f , the tube appears to 
operate in more of a natural-convection heating mode as no 
new nucleation sites are activated, although previously 
activated sites continue to nucleate. This phenomenon is 



probably 


due to 


the i nab i 1 i ty 


of the 


cartridge heaters 


to 


produce 


un i f or m 


heat flux at 


these 


1 owe r settings. 


At 


point C ’ 


with 


a heat -flux 


sett ing 


of 5 kW/m 2 , more 


nuc 1 ea t i on 


sites 


are generated 


and as 


the heat flux 


i s 


i ncr eased 


to 60 kW/m 2 , the tube 


becomes 


fully nuc 1 ea ted . 





For the 3 percent oil concentration, the point of 
incipient nucleate boiling, point B f , decreases as the non- 
linear physical properties of the R-114-oil mixture allow 
the onset of the mixed boiling region to occur. Also, the 
region of the mixed boiling, point B* to D # , is 
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Figure 4.12 Heat-Transfer Performance for Thermoexcel-E Surface. 



far larger than that for pure R-114, as the generation of 
new nucleation sites is retarded by the increased oil 
content. 

For the 10 percent oil concentration, on the other 
hand, the point of incipient nucleate boiling, point B tf , 
is delayed considerably. At approximately 14 kW/m 2 , the 
transition from natural -convective heating to nucleate 
boiling occurs. The exact reason for the difference between 
the 3 percent and 10 percent behavior is not known, however 
it is possibly due to the change in the physical properties 
of the R-114-oil mixture. 

The region of nucleate pool boiling for the condition 
of decreasing heat flux represented by points D to E, D f 
to E f and D Tt to E tT for 6, 3 and 10 percent oil 
concentrations, respectively, for the Thermoexce 1 -E tube, 
represent classic nucleate pool-boiling characteristics. 
As the oil concentration increases, the wall superheat also 
increases throughout the range of heat-flux settings. At 
the lower heat fluxes, the largest increase in the wall 
superheat occurred at 3 percent oil concentration when 
compared to the pure R-114, while the 10 percent oil 
concentration was observed to undergo continual degradation 
in tube performance over the entire range of heat-flux 
settings, although a more noticeable reduction occurred at 
the higher heat fluxes. At a heat-flux setting of 91 
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kW/m 2 , 



the R-114-oil mixture of 10 percent oil was observed 



to undergo a reduction in tube performance by a factor of 
3. For the practical heat-flux range of 30 to 40 kW/m 2 , 
there is an increase in wall superheat by factors of 1.4 
and 2.6 for R-114-oil mixtures of 3 and 10 percent, 
respect i ve 1 y . 

Figure 4.13 shows the heat-transfer coefficient of the 
Thermoexce 1 -E tube in ref r igerant-oi 1 mixtures, plotted as 
a function of heat flux. As can be seen in the figure, the 
heat- transfer coefficient remains nearly constant 
throughout the range of heat-flux settings for pure R-114, 
while there are slight reductions for the 1, 2 and 3 
percent oil concentrations. Additionally, it can be seen 
that at these lower oil concentrations, the hea t - t r ans f e r 
coefficients are nearly coincident. As the R-114-oil 
mixture oil concentration is increased to 6 percent, the 
tube performance increases with increasing heat flux until 
reaching a heat-flux setting of 60 kW/m 2 , at which time 
tube performance begins to deteriorate. This is contrary 
to the tube performance for 0, 2, and 3 percent oil 
concentrations in which the heat-transfer coefficient was 
observed to increase slightly over the range of heat 
fluxes. For the 10 percent oil concentration, the heat- 
transfer coefficient decreased as the heat flux was 
increased . 
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Figure 4.13 Boiling Heat-Transfer Coefficient for Thermoexcel-F Surface 



Figure 4.14 shows a direct comparison of the 
Thermoexce 1 -E tube to the smooth copper tube for R-114-oil 
mixtures of 0, 3 and 10 percent. The Thermoexce 1 -E tube 
presents enhancements of 6.6 and 4.1 over the performance 
of the smooth copper tube in pure R-114 and an R-114-oil 
mixture of 10 percent oil, respectively, at a heat-flux 
setting of 40 kW/m 2 . The results for the pure R-114 compare 
favorably with the results obtained by Yilmaz, Palen and 
Taborek [Ref. 19] for the Thermoexce 1 -E in p-xylene. 

Figure 4.15 shows the degradation that oil causes in 
the boiling hea t - t r an s f e r performance for the Thermoexce 1 -E 
tube. This figure plots the ratio of the hea t - t r ans f e r 
coefficient of the Thermoexce 1 -E tube with oil to the value 
without oil as a function of oil concentration. The shapes 
of the curves continue to demonstrate the non-linear 
variation of the physical properties of refrigerant-oil 
mixtures. However, the significant degradation of the tube 
performance is probably due to the inability of the 
generated bubbles to scavenge the interior tunnels of the 
Thermoexce 1 -E tube. This is in contrast to the GEWA-T 
tube, where more scavenging is possible due to its larger 
channel size. 

E. BOILING PERFORMANCE OF THE THERMOEXCEL -HE TUBE 

Figure 4.16 shows the nucleate pool-boiling 
performance of the Thermoexce 1 -HE tube in R-114-oil 
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Figure 4.14 Boiling Heat-Transfer Coefficient for Thermoexcel-E vs. 
Smooth Surface . 
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Figure 4.15 Relative Effect of Oil On Thermoexcel-E Boiling 
Heat-Transfer Performance. 



mixtures. In the region between points A and B/B f , 
natural-convection heating occurs for increasing heat flux. 
Although similar in performance to the Thermoexce 1 -E tube, 
the Thermoexce 1 -HE tube presented a more linear 
characteristic than that of the Thermoexce 1 -E tube, 
representing a more uniform heat distribution. At point B 
on the figure, the point of incipient nucleate boiling for 
pure R-114 liquid occurs at approximately 2 kW/m 2 . For 
R-114-oil mixture of 3 and 10 percent oil concentrations, 
the point of incipient nucleate boiling, point B f , occurs 
at approximately 8 kW/m 2 . Also in this region, the pure 
R-114 and the R-114-oil mixture of 3 percent oil 
concentration are nearly coincident in the range of heat 
flux from 500 W/m 2 to 2 kW/m 2 . 

From points B to C, the The r moex ce 1 -HE tube in the 
pure R-114 liquid achieved full nucleate boiling at 
approximately 37 kW/m 2 , with intermediate data points in 
the region of mixed boiling. These intermediate data 
points indicate that the Thermoexce 1 -HE tube was not able 
to activate all nucleation sites simultaneously for 
increasing heat-flux conditions. The R-114-oil mixture of 
10 percent oil was able to achieve activation of all 
nucleation sites upon reaching 22 kW/m 2 , without any data 
points contained in the mixed boiling region. With 3 
percent oil, the The r moe xce 1 -HE tube was able to achieve 
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Figure 4.16 Heat -Transfer Performance for Thermoexcel-HE Surface . 



full nucleate boiling at 37 kW/m 2 . 



The Thermoexcel -HE 



tube in the R-114-oil mixtures appeared to activate 
all nucleation sites simultaneously. 



The 


region from points 


C to D 


and 


from C f to D' 


represent 


the nucleate* pool 


bo i 1 i ng 


for 


increasing heat 



flux. The performance of the Thermoexcel -HE tube in this 
region appears to decrease as the oil content was 
increased. The 10 percent oil concentration increased the 
wall superheat by a factor of approximately 2 above the 
wall superheat for the pure R-114 liquid for heat-flux 
settings of 40 kW/m 2 . 

The regions between the points D and E and D’ and E 1 
represent nucleate pool-boiling for decreasing heat-flux 
conditions. As can be seen from the figure, the greatest 
reduction in performance occurs in the lower heat-flux 
settings between the pure R-114 and the 3 percent oil 
concentration. This is similar to the results for the 
Thermoexce 1 -E tube, where the tube performance decreased 
as oil content increased. There is a significant decrease 
in the The r moe x ce 1 -HE tube performance over the range of - 
heat-flux settings as the oil concent rat ion is increased to 
10 percent. 

Figure 4.17 shows the hea t - t r ans f e r coefficient of the 
Thermoexce 1 -HE tube in refrigerant-oil mixtures, plotted 
as a function of heat flux. As can be seen, the 
performance is very similar to that of the Thermoexcel -E 
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tube (Figure 4.13). 



As was observed for the The r moe xce 1 -E 



tube, the The rmoe xce 1 - HE tube undergoes a slight reduction 
in the heat-transfer coefficient as the oil concentration 
was increased from 1 to 3 percent. As the oil 
concentration is increased to 6 and 10 percent oil, the 
performance reduction is quite noticeable. With the 10 
percent oil concentration, tube performance decreases at 
both the lower and higher heat-flux settings. This is 
contrary to the findings for the Thermoexce 1 -E tube where 
the performance decreased over the entire range of heat- 
flux settings, from low heat flux to high heat-flux 
sett i ngs . 

Figure 4.18 shows a comparison of the Thermoexce 1 -HE 
tube to the smooth copper tube at R-114-oil mixtures of 0, 
3 and 10 percent oil concentrations. The The rmoe xce 1 -HE 
tube displayed an enhancement of 6.1, 4.6 and 4.1 at a heat 
flux setting of 40 kW/m 2 and R-114-oil mixtures of 0, 3 and 
10 percent oil, respectively, over the same parameters for 
the smooth copper tube. These results are nearly identical 
to those found for the Thermoexce 1 -E tube. 

Figure 4.19 shows the degradation that the oil has on 
the boiling heat-transfer performance of the Thermoexce 1 -HE 
tube. At the lower heat fluxes, the effect of the oil 
concentration is slightly more pronounced than that 
observed for the Thermoexce 1 -E tube throughout the range of 
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Figure 4.17 Boiling Heat-Transfer Coefficient for Thermoexcel-HE . 
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Figure 4.18 Boiling He at -Transfer Coefficient Comparison for 
Thermoexcel-HE Surface vs. Smooth Surface . 
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Figure 4.19 Relative Effect of Oil on Thennoexcel-HE 
Boiling Heat -Transfer Performance. 



oil concentrations, although more severe at the lower oil 
concentrations. At the higher heat fluxes, the trend 
appears to reverse. 

F. COMPARISON TO REILLY’S DATA FOR THE HIGH FLUX TUBE 

Figures 4.20, 4.21 and 4.22 provide a comparison of 
the smooth, High Flux, GEWA-T, Thermoexce 1 -E and 
Thermoexce 1 -HE tubes for 0, 3 and 10 percent oil 
concentrations, respectively. For pure R-114 liquid, the 
High F 1 ux , Thermoexce 1 -E and The rmoexce 1 -HE tubes operating 
at a heat flux of approximately 20 kW/m 2 , show nearly 
identical performance. As the heat flux was increased, the 
performance of the High Flux tube continued to 
increase while the performance of Thermoexce 1 -E and 
Thermoexce 1 -HE tubes remained nearly constant. At heat- 
flux settings lower than 20 kW/m 2 , the Thermoexce 1 -E and 
Thermoexce 1 -HE tubes showed higher heat-transfer 
coefficients than that for the High Flux tube. As the oil 
concentration was increased to 3 percent, the performance 
of all five tubes decreased, as can be seen in Figure 4.21. 
Additionally, the GEWA-T tube’s performance appeared to 
decrease the least since, as discussed earlier in section 
IV. C, the finned surface of the GEWA-T tube is not 
adversely effected by the small increase in the oil content 
due to its large reentrant cavities. From this 
observation, it can be stated that the performance of the 
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Figure 4.20 Comparison of Boiling Heat-Transfer Coefficients for 
Smooth, High Flux [Ref. 5] GEWA-T, Thermoexcel-E and 
Thermoexcel-IIE Surfaces in Pure R-114. 
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Figure 4.21 Comparison of Boiling Meat-Transfer Coefficients for 
Smooth, High Flux [Ref. 5] , GEWA-T, Thermoexcel-E 
and Thermoexcel-llE Surfaces for 3 Percent Oil Content. 
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Figure 4.22 Comparison of Boiling Heat-Transfer Coefficients for 
Smooth, High Flux [Ref. 5] , GEWA-T, Thermoexcel-E 

and Thermoexcel-HE Surfaces for 10 Percent Oil Content. 



tubes with porous coatings or small reentrant cavities is 
adversely affected by the increase in oil content due to 
the generation of an oil-rich film in the nucleation sites 
impeding the tube’ s performance. This is further evidenced 
by Figure 4.22* which depicts the performance of all five 
tubes in R-114-oil mixtures of 10 percent oil. The 
performance of all the tubes is reduced significantly, 
while the performance of the High Flux tube is severely 
decreased in the range of heat-flux settings from 22 kW/m 2 
to the maximum setting of 91 kW/m 2 . At approximately 45 
kW/m 2 , the performance of the High Flux tube, the GEWA-T, 
the Thermoexce 1 -E and the Thermoexce 1 -HE tubes are 
approximately equal, although the performance of the 
Thermoexce 1 -E and the Thermoexce 1 -HE tubes remain nearly 
constant over the range of heat-flux settings of 14 kW/m 2 
to 50 kW/m 2 . Although the performance of the High Flux 
tube exceeds that of the GEWA-T, the Thermoexce 1 -E and the 
Thermoexce 1 -HE tubes over this range, its performance 
begins to severely degrade at the midpoint of the heat-flux 
range indicated. In the heat-flux range above 45 kW/m 2 , 
the GEWA-T tube outperformed all the other tubes. This is 
attributed to the ability of bubbles generated in its 
larger channels to continually scavenge the oil-rich layer 
from the tube, and therefore remove the insulating oil from 
the tube’s channels and exterior surface. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. In pure R-114, at a heat flux of 30 kW/m 2 the 

pool-boiling heat- transfer coefficient of the 
Thermoexce 1 -E , Ther moexce 1 -HE and GEWA-T tubes 
were approximately 6.6, 6.1 and 4.2 times, 

respectively, larger than the baseline smooth tube. 

2. In R-114-oil mixtures of 3 percent oil, at a 

heat flux of 30 kW/m 2 , the poo.l -bo i 1 i n g heat- 
transfer coefficient of the Ther moexce 1 -E , Thermo- 
exce 1 -HE and GEWA-T tubes were approximately 4.8, 
4.6 and 4.0 times, respectively, larger than the 
smooth tube under the same conditions. 

3. In R-114-oil mixtures of 10 percent oil, at a heat 

flux of 30 kW/m 2 , the performance of the Thermoexcel- 
E and Thermoexce 1 -HE tubes were nearly identical, 
with an enhancement of 4.1 over the pool-boiling 
heat- transfer coefficient of the smooth tube.- The 
GEWA-T tube, under the same conditions, produced an 
enhancement of 3.8 over the smooth tube. 

4. As oil content in the refrigerant-oil mixture 

increased, the performance of enhanced surfaces 
decreased; however, the amount of this decrease was 
directly related to the size of the subsurface 
channels (i.e., the larger the reentrant cavity or 
channel, the less adverse effect of the increased oil 
content ) . 

5. As the oil content increases in the channels, the 
greater the adverse effect of the oil-rich layer has 
on the surface, insulating the surface from the 
refrigerant and therefore decreasing the performance 
of the surface. 

6. The GEWA-T surface in refrigerant-oil mixtures, with 
its larger subsurface channels, was able to 
scavenge the oil-rich layer from its cavities by 
the bubbling action generated within the channels. 

7. The performance of the Ther moexce 1 -E and Thermoexcel- 
HE tubes, with their smaller pore size and smaller 
tunnels, showed significant performance decreases 
as oil content was increased. The smaller pore and 
tunnel size under these conditions, were unable to 
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effectively disperse the oil-rich film from the sub 
-surface cavities and therefore allowing the oil-rich 
film to insulate the tube’s surface. 



B. RECOMMENDATIONS 

1. The single-tube testing should be expanded to include 
multi-tube bundles. This would provide more realistic 
data with regard to the effects of tubes operating in 
close proximity to each other. 

2. A uniform heat flux is necessary within the boiling 
tube and within the boiling section itself. There- 
fore, the present cartridge heaters need to be 
replaced with more precise heaters, capable of pro- 
ducing a uniform heat flux in the axial direction. 

3. The testing of surfaces should be expanded to study 
the effects of various refrigerant oils. 

4. Conduct experiments to establish the physical proper- 
ties of refrigerant-oil mixtures to enable further 
understanding into the boiling mechanisms involved. 

5. There needs to be additional studies undertaken on 
the effects of additional boiling temperatures on the 
mixture properties of the refrigerant-oil mixtures. 
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APPEND I. X ._.A 

THERMOCOUPLE CALIBRATION [Ref. 53 

Karasabun [Ref. 17] described the thermocouple 
calibration procedure in great detail. Reilly [Ref. 53 
stated that the data-reduction program utilized differences 
in thermocouple readings in all computations, such as the 
amount of the tube wall superheat (Two - Tsat). Reilly 
stated that the calibration of the thermocouples would be 
necessary if the thermocouples were measuring quantities 
such as absolute temperature to determine actual fluid 
properties. The procedure utilized by Karasabun and Reilly 
for the initial instrumentation for the experimental 
apparatus is as contained in Appendix A of Reference 5. 

To determine the calibration for the thermocouples, 
two thermocouples were made, one from the beginning of the 
wire roll and the second from the end. The seventh-order 
polynomial equation provided by the wire manufacturer was 
corrected by adding to it a second-order curve fit of the 
variation of the manufacturer's computed temperature for a 
given emf from a known set of reference temperatures as 
measured by a Hewlett-Packard 2Q04A quartz thermometer. The 
thermometer had a temperature resolution of -0.0001 K and 
an accuracy of to. 03 K. 
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The 



manufacturer 1 s emf 



to temperature 



convers 1 on 



equa t i on i s : 

T = a Q + a t E + a ? E 2 + a 3 E 3 + a 4 E 4 + (A.i) 

a 5 E s + a* E 6 + a 7 E 7 

where 

E = thermocouple reading in volts 

T = temperature 

a 0 = 0.100860910 

at = 25727.94369 

a 2 = -767345.8295 

a 3 = 78025595.81 

a 4 = “9247486589 

a, = 6. 97686xlO J 1 

a A = -2 . 66 1 92x 1 0 1 3 

a 7 = 3 . 94078x 1 0 1 4 

Figure A.I [Ref. 51 shows the difference in the quartz 
thermocouple readings and the actual thermocouple readings 
(discrepancy) plotted as a function of temperature. As can 
be seen from the figure, the two thermocouples agreed to 
within 0.05 K of each other. The manufacturer’s 

calibration equation needed an increase of approximately 
0.1 K to more accurately convert the emf’ s to the true 
temperature. The second-order polynomial correction was as 
foil ows ; 

DCP = b 0 + biT + b 2 T 2 ( A . 2) 

where 
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DCP = 


discrepancy (K) 


bo = 


8.6268968x10“ 2 


bi = 


3.7619902x10“ 3 


b? = 


-5.0689259x10“ 5 


T 2 = 


thermocouple readings from 
equa tion A.I in C 


Therefore , 


the temperatures computed by the da t a - r educ t i on 


progr am 


(DRP5) were emf f s converted directly to 



temperatures by equation A.i and corrected by equation A. 2. 
This calibration was. thus estimated to result . in 
temperature readings within ±0.05 K of the true 
temperature . 
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Figure A.l Thermocouple Discrepancy Correction. 



APPENDIX B 



DATA-REDUCT I ON PROGRAM CREF. 5] 



The da ta-reduct i on program below consists of the 
following sections: 

Main Program - Menu of subprogram options 

Sub Main - Takes or reprocesses data 

Sub plot - Plots data on log-log scale 

Sub Poly - Computes least squares curve fit 

Sub Plin - Plots data on 1 i near - 1 i near scale 

Sub Stats - Computes average and standard 

deviation of data 

Subprogram Main consists of the following steps: 

1. Creates data file for data and plotting file. 

2. Annotates tube type. 

3. Monitors heat flux or saturation temperature to 
establish steady-state conditions. 

4. Scans HP3497 channels and stores in data file. 

5. Converts raw emf f s to temperatures, current and 
voltages. 

6. Computes the hea t- t r ans f er rate for the tube’s 
internal cartridge heater. 

7. Computes the average wall temperature of the 

tube, the average wall superheat and 

the film temperature. 
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8. Computes the physical properties of the 
R-114 using the given correlation temperature. 

9. Computes the natural-convection hea t - t ransf er 

coefficient of R-114 for the non-boiling 

tube ends. 

10. Computes the heat loss from the non-boiling 
tube ends. 

11. Calculates the corrected heat fluxes from 
the testing tube to the R-114. 

12. Calculates the boiling hea t - t r ans f er coefficient 
of the R-114 from the tube being tested. 

13. Prints the data and stores the heat flux and tube 
wall superheat values in the plot file. 

The following is the listing of the complete da ta- reduct i on 
program (DRP5) written in Basic 3.0 for the Hewlett-Packard 
9826 computer. 
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1000! FILE NAME : DRPS 
1003! DATE: Octooer 19. 1934 

1006! REVISED: April 25. 1986 

1009! 

1012 COM /Idp/ IdP 
1015 PRINTER IS 1 
1018 CALL Select 

1021 INPUT "WANT TO SELECT ANOTHER OPTION ( 1 - Y . 0-N)?* , Isel 

1024 IF Isei-1 THEN GOTO 1018 

1027 BEEP 

1030 BEEP 

1033 PRINTER IS 1 

1036 PRINT “DATA COLLECTION/REPROCESSING COMPLETED* 

1039 END 

1042 SUB Main 

1045 COM /Idp/ Idp 

1048 COM /Cc/ C( 7) . leal 

1051 COM /Wil/ D2.Di .Do.L.Lu.Kcu 

1054 DIM Emf < 12> .T( 12) ,D1a( 1 1 ) .D2a( 1 1 ) .Dial 1 1 ) .Doa( 1 1 ) .La( 11 ) ,Lua< 1 1 ) ,Kcua( 11). 
Ett 19) , Tn$(4) [151 

1057 DATA 0. 10086091 .25727.94369.-767345. 9295. 78025595. 81 
1060 DATA -9247486589 . 6 . 97688E+1 1 ,-2.661 92E+1 3 , 3 . 94078E+1 4 
1063 READ CC*) 

1066! DATA “Smooth* , "High Flux* , ’Thermoexel-E* , *Thermoexel-HE" 

1069 DATA Smooth, High Flux, Turbo-B.Hioh Flux Mod.Turbo-B Mod 
1072 READ Tn$(+> 

1075 PRINTER IS 701 
1078 BEEP 

1081 IF Idp-4 THEN PRINTER IS 1 
1084 IF Idp-4 THEN GOTO 1774 

1087 INPUT "ENTER MONTH. DATE AND TIME (MM:DD:HH:MM:SS) * .DateS 
1090 OUTPUT 709 : * TD* :DateS 
1093 OUTPUT 709: “TD* 

1096 ENTER 709:Date$ 

1099 PRINT 

1102 PRINT " Month, date and time :’:Date$ 

1105 PRINT 

1108 PRINT USING "10X."‘N0TE: Program name : DRP4*** 

1111 BEEP 

1114 INPUT "ENTER DISK NUMBER*, Dn 

1117 PRINT USING * 1 6X . * *D isle numoer - **.ZZ*;Dn 

1120 BEEP 

1123 INPUT ’ENTER INPUT MODE ( 0-30S4A. 1 -FILE) * , Im 
1126 BEEP 

1129 INPUT ’SELECT HEATING MODE (0-ELEC: 1 -WATER) *, Ihm 
1 1 32 BEEP 

1135 INPUT ‘ENTER THERMOCOUPLE TYPE (0-NEW, 1-OLD)*, leal 
1138 IF Im-0 THEN 
1141 BEEP 

1144 INPUT "GIVE A NAME FOR THE RAW DATA FILE" ,D2_f i leS 
1147 PRINT USING *16X,**Neu file name: “ * , 1 4A* :D2_f i leS 
1150 S i ze 1 -20 

1153 CREATE BOAT D2_f ileS.Sizel 
1156 ASSIGN 'SF 1 1 e2 TO D2_fileS 
1159! 

1162! DUMMY FILE UNTIL Nrun KNOWN 
1165 Dl_f tleS-"DUMMY“ 

1163 CREATE BOAT 01 f ileS.Sizel 
1171 ASSIGN AF i lei TO D1_file$ 

1174 OUTPUT «Fi lei ;Date£ 
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1177 IF Ihm-0 THEM 
1180 BEEP 

1183 INPUT “ENTER NUMBER OF DEFECTIVE TCS ( O-DEFAULT) “. Idtc 
1186 IF Idtc-0 THEM 
1189 Ldtcl *0 
1192 Ldtc2-0 

1195 PRINT USING “16X.“"No defective TCs exist"*" 

1198 END IF 

1201 IF Idtc* 1 THEN 

1204 BEEP 

1207 INPUT “ENTER DEFECTIVE TC LOCATION* , Ldtcl 

1210 PRINT USING *16X.“*TC is defective at location **,D“:Ldtc1 

1213 Ldtc2-0 

1216 END IF 

1219 IF Idtc-2 THEN 

1222 BEEP 

1225 INPUT “ENTER DEFECTIVE TC LOCATIONS* , Ldtcl ,Ldtc2 

1228 PRINT USING *16X,“*TC are detective at locations ** ,D.4X,D* ;Ldtc1 ,Ldtc2 

1231 END IF 

1234 IF Idtc>2 THEN 

1237 BEEP 

1240 PRINTER IS 1 

1243 BEEP 

1246 PRINT “INVALID ENTRY* 

1249 PRINTER IS 701 

1252 GOTO 1180 
1255 END IF 

1253 END IF 

1261 OUTPUT 5>Fiie1 ; Ldtcl ,Ldtc2 
1264! Im-1 option 
1267 ELSE 
1270 BEEP 

1273 INPUT “GIVE THE NAME OF THE EXISTING 'DATA FILE" ,D2_f ileS 
1276 PRINT USING *16X.“*01d file name: * ' , 1 4A* : D2_f i ieS 
1279 ASSIGN »File2 TO D2_fiieS 
1282 ENTER (»File2:Nrun 
1285 ENTER 3File2:Dold$ 

1288 PRINT USING “16X.‘*This data set taken on : ** , 14A";Dold$ 

1291 ENTER »File2:Ldtc1 ,Ldtc2 
1294 IF Ldtcl >0 OR Ldtc2>0 THEN 

1297 PRINT USING “16X. "“Thermocouples were defective at locations: “* ,2(3D,4X) 

Ldtcl ,Ldtc2 

1300 END IF 

1303 ENTER »File2:Itt 

1306 END IF 

1309! IF Im-0 AND Ihm-1 THEN 1595 
1312 BEEP 

1315 INPUT “WANT TO CREATE A PLOT FILE? (0-N, 1-Y) * , Iplot 
1318 IF Iplot-1 THEN 
1321 BEEP 

1324 INPUT “GIVE NAME FOR PLOT FILE’ .P_f i IeS 

1327 CREATE BOAT P_fileS,4 

1330 ASSIGN SPlot TO P_fileS 

1333 END IF 

1336 IF Ihm-1 THEN 

1339 BEEP 

1342 INPUT “WANT TO CREATE Uo FILE? (0-N, 1-Y)“.Iuf 
1345 IF Iuf-1 THEN 
1348 BEEP 

1351 INPUT “ENTER Uo FILE NAME" ,Uf i IeS 
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1354 

1357 

1380 

1383 
1366 
1363 
1372 
1375 
1378 

1381 

1384 
1387 
1390 
1393 
1396 
1399 
1402 
1405 
1408 
1411 
1414 
1417 
1420 
1423 
1426 
1429 
1432 
1435 
1438 
1441 
1444 
1447 
1450 
1453 
1456 
1459 
1462 
1465 
1468 
1471 
1474 
1477 
1480 
1483 
1486 
1489 
1492 
1495! 
1498! 
1501 ! 
1504! 
1507 
1510 
1513 
1516 
1519! 
1522! 
1525 
1528 
1531 



CREATE BOAT IJfileS.4 
ASSIGN SUfile TO UfileS 
END IF 
BEEP 

INPUT "WANT TO CREATE Re FILE? ( 0-N. ! -Y ) * . Ire 

IF Ire-1 THEN 

BEEP 

INPUT "ENTER Re FILE NAME" , Ref i lei 
CREATE BDAT Retries. 10 



ASSIGN jRefile TO RefileS 



ENO IF 
END IF 
PRINTER IS 1 
IF Im-0 THEN 
BEEP 

PRINT USING "4X, 
IF Ihm-0 THEN 
PRINT USING “6X. 
PRINT USING "6X . 
PRINT USING *6X. 
PRINT USING “6X. 
PRINT USING *SX. 
PRINT USING *SX . 
PRINT USING *SX, 
PRINT USING *6X. 
PRINT USING "SX. 
PRINT USING *6X. 
PRINT USING "6X . 
PRINT USING "6X , 
ELSE 

PRINT USING *6X, 
PRINT USING *6X. 
PRINT USING “6X. 
PRINT USING "SX. 
PRINT USING "SX. 
END IF 
INPUT Itt 



""Select tube number"’* 

*"0 Smooth 4 inch Ref* 

""1 Smooth 4 inch Cu (Press/Slide)""" 
""2 Soft Solder 4 inch Cu""" 

"*3 Soft Solder 4 inch HIGH FLUX*"" 

* *4 Uie land Hard 8 inch""* 

"*5 HIGH FLUX 8 inch""" 

""8 GEWA-K 19 Fins/in*"" 

" *7 GEWA-K 26 Fins/in""" 

"*3 GEWA-T 19 Fins/in"*" 

""9 GEWA-T 26 Fins/in""" 

""10 THERMOEXCEL-E""" 

""11 THERMOEXCEL-HE""" 

""0 Smooth tube""* 

""1 Hign Flux""" 

"*2 Turbo-B""* 

""3 High Flux Mod"** 

"*4 Turbo-3 Mod""" 



OUTPUT SFilel lltt 
END IF 

PRINTER IS 701 

IF Itt<10 THEN PRINT USING "16X.""Tube Number: "*.D":Itt 

IF Itt>9 THEN PRINT USING "16X.""Tube Number: *".DD";Itt 

IF Ihm-1 THEN PRINT USING "16X,**Tube Type: ** . 15A* : TnSt Itt) 

SEEP 

INPUT "ENTER OUTPUT VERSION ( 0-LONG. 1 -SHORT ,2-NONE) ", Iov 
BEEP 

INPUT "SELECT (0-LIQ. 1 -VAP ,2-lLIQ+VAP) /2) " . Ilqv 



DIMENSIONS FO TESTED TUBES 
ELECTRIC HEATED MODE 

Di-Diameter at thermocouple positions 

DATA .011 1125.. 01 11 125. .01 11125. .0129540. .012446. .0129540. .0100965 
DATA .0100965. .01 157 , . 0 1 157 . . 01 157 . . 01 1 57 
READ Dial*) 

Dl-Ola(Itt) 



02-Diameter of test section to the base of fins 
DATA .015875.. 015875. .015375. .015824.. 01 5875. .015824.. 01270 
DATA .0127. .0138, .0138. .0138. .0138 
READ 02a(*> 
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1534! 

1537 1 DL'Tnside diameter of unennanced ends 

1 S40 DATA .01 27, .01 27, .0127, .01 32. .0127, .01 32. .01111 25, .01111 25. .0118, .0118, .01 

18. .01 1U 

1543 READ Dial * ) 

1546! 

1549 1 Dd"0utside diameter of unenhanced ends 

1S52 - DATA .01 5875,. 01 5375. .015875, .015824, .015875. .015824, .01270, .01270, .01331 , 
.01331, .01331. .01331 
1555 READ Doa( *) 

1558! 

1561 1 L'Cength of enhanced surface 

1 564 DATA .1016. .1016, .1016, .1016, .2032. .2032 . .2032, .2032. .2032, .2032, .2032 , .20 
32 

1567 READ laU) 

1570! 

1573 • Lu'Cength of unenhanced surface at the ends 

1576’ DATA .0254, .0254, .0254, .0254, .0762.. 0762. .0762, .0762. .0762, .0762, .0762, .07 
62 

1S79 READ Lua( * ) 

1582! 

15351 Kcu-Thermai Conductivity of tube 

1588’ DATA 398,344,344.45,344,45,344,944,393,393,393,398 

1591 READ Kcua(+) 

1594 IF Ihm-1 THEN 
1597! 

1600* Data statements for water heating mode 
1603! 

1606 DATA 0. 015375, 0.01537S, 0.0169. 0.0138. 0.0169, 0.0. 0,0.0 
1609 READ D2a( *) 

1612 DATA 0.0127,0.0127,0.0145,0.0127,0.0145.0,0,0,0.0 
1615 READ Dia(*) 

1613 DATA 0.015875.0.015375,0.0169,0.015875,0.0169,0.0,0,0,0 
1621 READ Doa(*) 

1624 DATA 0.3048,0.3048,0.3048,0.3048,0.3048.0,0,0,0.0 
1627 READ La(*> 

1630 DATA 0.02S4. 0.0254, 0.0254. 0.0254, 0.0254, 0.0, 0,0,0 
1633 READ lua(*> 

1636 DATA 398.45,398,45,398.0,0,0,0.0 
1639 READ Kcua(*) 

1642 END IF 
1645 D2“D2a(Itt) 

1648 DfDiadtt) 

1651 0o“Doa(Itt) 

1654 L'Ea(Itt) 

1657 Lu"Lua(Itt) 

1660 Kcu-Kcua( Itt) 

1663 Xn-.8 

1666 Fr-.3 

1669 IF Itt-0 THEN Cf-l.70E+9 

1672 IF ltt>0 THEN Cf -3 . 7037E+1 0 

1675 A'PI* ( 0o‘2-Di*2)/4 

1673 P-PI*Do 

1681 IF Ihm-1 THEN 
1684 BEEP 

1687 INPUT “TUBE INITIATION MODE. (1-H0T WATER. 2-STEAM, 3-COLD WATER) *, It im 
1690 IF Itim-1 THEN PRINT USING ,, 16X.*"Tube Initiate! Hot Water""* 

1693 IF Itira-2 THEM PRINT USING * 1 6)< . " "Tube Initiate! Steam""" 

1696 IF Itim-3 THEN PRINT USING "l6X.""Tube Initiate! Cold Water""" 

1699 INPUT "TEMP/VEL MODE: ( 0-T-C0NST. V-DEC: I -T-DEC.O-CONST; 2-T-INC, 0-CONST) " . 
I tv 
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1702 IF Itv-0 THEM PRINT USING “ 16X. “ “Temp/Vel Mode: Constant/Decreasing' 1 “ 

1705 IF Itv-1 THEN PRINT USING "1SX. ““Temp/Vel Mode: Deoreas ing/Cons tan t " ” 

1703 IF I tv- 2 THEN PRINT USING “16X. ""Temp/Vel Mode: Increasing/Constant"" 

1711 INPUT "WANT TO RUN WILSON PLOT? ( 1 -Y , 0-N) “ , Iwi 1 

1714 IF Ihm-1 AND Iuul-0 THEN 

1717 IF Itt-0 THEN Ci-,032 

1720 IF Itt-I OR Itt-3 THEM C1-.0S3 
1723 IF Itt-2 OR Itt-4 THEN Ci-.062 
172S BEEP 

1729 INPUT “ENTER Cl (DEF: WH-.032.HF-.059, TB-.0S2) * ,Ci 

1732 PRINT USING “16X. ““Sieder-Tate 

1735 PRINT USING “16X,“* Constant - "".Z.4D":Ci 

1738 END IF 

1741 END IF 

1744 IF Ihm-1 AND Im-1 AND Iuil-1 THEN 
1747 IF Itt-0 THEN C1-.032 
1750 IF Itt-1 OR Itt-3 THEN CI-.0S9 
1753 IF Itt-2 OR Itt-4 THEN CI-.062 
1756 ASSIGN »File2 TO + 

1759 CALL WilsonfCf.Ci) 

1762 ASSIGN SFile2 TO D2_fileS 

1765 ENTER 4File2;Nrun,Dold$.Ldtc1 ,Ldtc2, Itt 

1763 END IF 
1 771 Nsub-0 

1774 IF Idp-4 THEN Ihm-1 
1777 IF Ihm-1 THEN Nsub-8 
1780 Ntc-6 

1783 IF Ihm-0 THEN Ntc-12 

1786 J-1 

1789 Sx-0 

1792 Sy-0 

1795 Sxs-0 

1798 Sxy-0 

1801 Reoeat: ! 

1804 IF Im-0 THEN 
1807 Dtid-2.22 

1810 Ido-2 

1313 ON KEY 0.15 RECOVER 1801 
1816 PRINTER IS 1 

1819 PRINT USING *4X. “"SELECT OPTION"*" 

1822 PRINT USING *6X.*"0-TAKE DATA”"" 

1325 IF Ihm-0 THEN PRINT USING “6X."*1-SET HEAT FLUX*"" 

1828 IF Ihm-1 THEN PRINT USING ”6X.""1-SET WATER FLOW RATE""" 

1831 PRINT USING ”6X. "“2-SET Tsat""" 

1834 PRINT USING "4X.""N0TE: KEY 0 - ESCAPE’"" 

1837 BEEP 

1840 INPUT Ido 

1843 IF Ido>2 THEN Ido-2 

1846 IF Ido-0 THEN 2296 

1849! 

1852! LOOP TO SET HEAT FLUX OR FLOWMETER SETTING 

1855 IF Ido-1 THEN 

1358 IF Ihm-0 THEN 

1361 OUTPUT 709: “AR AF12 AL13 VR5" 

1864 BEEP 

1867 INPUT “ENTER DESIRED Odo’.Dadp 

1870 PRINT USING "4X. ““DESIRED Qdp ACTUAL Odp""" 

1873 Err- 1 000 

1876 FOR 1-1 TO 2 

1879 OUTPUT 709: “AS SA“ 
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1382 Sum-0 
1885 FOR Ji-1 TO 5 
1888 ENTER 709:E 
1891 Sum-Sum+E 
1894 NEXT Ji 

1897 IF 1-1 THEN Volt-Sum*5 
1900 IF 1-2 THEN Amp-E 
1903 NEXT I 

1906 Aqdp-0olt+Amp/(PI*02*L) 

1909 IF A8S( Aqdp-Dqdp) >Err THEN 

1912 IF Aqdp>Dqdp THEN 

1915 BEEP 4000,-2 

1918 BEEP 4000,. 2 

1921 BEEP 4000,. 2 

1924 ELSE 

1927 BEEP 250,. 2 

1930 BEEP 250.. 2 

1933 BEEP 250,. 2 

1936 END IF 

1939 PRINT USING “4X.MZ.3DE,2X,MZ.3DE* :Dqdp,Aqdp 

1942 WAIT 2 

1945 GOTO 1876 

1948 ELSE 

1951 BEEP 

1954 PRINT USING *4X,MZ.3DE,2X,MZ.3DE";Dqdp,Aqdp 

1957 Err- 500 

1960 WAIT 2 

1963 GOTO 1876 

1966 END IF 

1969 ELSE 

1972 BEEP 

1975 INPUT “ENTER FLOWMETER SETTING* ,Fms 

1978 GOTO 1819 

1981 END IF 

1984 END IF 

1987! 

1990! LOOP TO SET Tsat 
1993 IF Ido-2 THEN 
1996 IF Ikdt-1 THEN 2011 
1999 BEEP 

2002 INPUT “ENTER DESIRED Tsat", Dt Id 

2005! PRINT USING “4X,““ DTsat ATsat Rate Tv Rate'*'* 

2008 Ikdt-1 

2011 Oldl-O 

2014 01d2-0 

2017 Nn-1 

2020 Nrs-Nn MOD 15 

2023 Nn-Nn+1 

2026 IF Nrs-1 THEN 

2029 PRINT USING “4X.“* Tsat TIdl Tld2 Tv Tsump Tiniet 
Tout"** 

2032 END IF 

2035 IF Ihm-0 THEN OUTPUT 709: “AR AF3 ALII VR5* 

2038 IF Ihm-1 THEN OUTPUT 709:“AR AFO AL5 VR5" 

2041 FOR 1-1 TO 6 

2044 IF Ihm-0 AND I>4 THEN 2125 

2047 Sum-0 

2050 OUTPUT 709: “AS SA“ 

2053 FOR Ji-1 TO 20 

2056 ENTER 709:Eliq 



Tpi le 
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2059 Sum-Sum+El iq 

2062 NEXT Ji 
2065 Eliq-Sum/20 

2063 Tld-FNTvsv(Eliq) 

2071 IF 1-1 THEN Tldl-TId 
2074 IF 1-2 THEN TId2-Tld 
2077 IF 1-3 THEN Tv-Tld 
2030 IF 1-4 THEN Tsump-Tld 
2033 IF 1-5 THEN Tmlet-Tld 
2086 IF 1-6 THEN Tout-Tld 
2089 NEXT I 

2092 IF Ihm-1 THEN 

2095 OUTPUT 709: “AR AF20 AL20 VRS” 

2098 OUTPUT 709: “AS SA* 

2101 Sum-0 

2104 FOR Klc-1 TO 20 

2107 ENTER 709 :E 

2110 Sum-Sum*E 

2113 NEXT Kk 

2116 Emt (7) -A8S( Sum/20) 

2119 Tpile-Emt (7J/3.96E-4 
2122 END IF 

2125 Atld-<Tld1+Tld2)*.S 

2128 IF ABS( At ld-Dt Id) > .2 THEN 

2131 IF Atld>DtId THEN 

2134 BEEP 4000,. 2 

2137 BEEP 4000,. 2 

2140 BEEP 4000,. 2 

2143 ELSE 

2146 BEEP 250.. 2 

2149 BEEP 250.. 2 

2152 BEEP 250.. 2 
2155 END IF 

2153 Errl-Atld-Oldl 

2161 0 ldl -At Id 

2164 Err2-Tv-01d2 
2167 01d2-Tv 

2170 IF Ihm-0 THEN PRINT USING "4X. 5 ( MDD . DD.2X) * :D t Id. TIdl . TId2 . Tv , Tsume 
2173 IF Ihm-1 AND Idp-0 THEN PRINT USING " 4X . 7 ( MDD . DD ,2X ) “ :Dt Id . T ldl , TId2 . Tv , Ts 
umo . T inlet . Tp i le 

21 76* IF Ihm-1 AND Idp-4 THEN PRINT USING “4X . SU1DD . DD.2X ) ,3(M3D . DD.2X ) “ ; D t Id . T1 

dl . Tld2, Tv. Tsump, Tin 

2179 WAIT 2 

2182 GOTO 2020 

2185 ELSE 

2188 IF ABS( A t ld-Ot Id) > . 1 THEN 

2191 IF Atld>Dtld THEN 

2194 BEEP 3000.. 2 

2197 BEEP 3000. .2 

2200 ELSE 

2203 BEEP 800.. 2 

2206 BEEP 300,. 2 

2209 END IF 

2212 Errl-Atld-Oldl 
2215 Oldl-Atld 

2213 Err2-Tv-01d2 

2221 01d2-Tv 

2224 IF Ihm-0 THEM PRINT USING “4X . S( MDD . DD .2X) * :Dt Id. Tldl . Tld2. Tv, Tsump 
2227* IF Ihm-1 THEN PRINT USING *4X .5 < MDD . DD ,2X ) ,3(M3D.DD . 1 X) ” ;Dt Id. Tldl . Tld2, Tv 
, Tsump, Tmlet.Tpile, 
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2230 WAIT 2 
2233 GOTO 2020 
223 G ELSE 
2239 BEEP 
2242 Errl-Atld-Oldt 
2245 Oldt-Atld 
2248 Err2-Tv-01d2 
2251 01d2-Tv 

2254 IF Ihm-0 THEN PRINT USING “4X , 5 ( (1DD . DO , 2X) * ;D t Id, T Id 1 , Tid2. Tv , Tsump 

2257 IF Ihm-1 THEN PRINT USING *4X ,3( MOD . DD ,2X) “ ; Dt Id , T Idl , Tld2 . Tv , Tsump . Tinle t 

.Tpile.Tout 

22S0 WAIT 2 

22S3 GOTO 2020 

2286 END IF 
2269 END IF 
2272 END IF 

2275! ERROR TRAP FOR Ido OUT OF BOUNDS 
2278 IF Ido>2 THEN 
2281 BEEP 
2284 GOTO 1819 

2287 END IF 
2290! 

2293! TAKE DATA IF Im-0 LOOP 
2296 IF Ikol-1 THEN 2308 
2299 BEEP 

2302 INPUT "ENTER BULK OIL 7.", Bop 
2305 Ikol-1 

2308 IF Ihm-0 THEN OUTPUT 709: ‘AR AFO ALII VR5" 

2311 IF Ihm-1 THEN OUTPUT 709;'AR AFO AL5 VRS" 

2314 IF Ihm-0 THEN Ntc-!2 
2317 FOR 1-1 TO Ntc 
2320 OUTPUT 709: “AS SA“ 

2323 Sum-0 
2326 FOR Ji-t TO 20 
2329 ENTER 709 :E 
2332 Sum-Sum+E 

2335 IF I-(9-Nsub) OR I-( 1 O-Nsubf THEN Et(Ji-t)-E 
2338 NEXT Ji 
2341 Kdl-0 

2344 IF I- ( 9-Nsub ) OR I.-(IO-Nsub) THEN 
2347 Eave-Sum/20 
2350 Sum-0. 

2353 FOR Jk-0 TO 19 

2356 IF A8S (Et( Jk ) -Eave) <5 . OE-6 THEN 

2359 Sum-Sum+E t ( Jk ) 

2362 ELSE 
2365 Kdl-Kdl+1 
2368 END IF 
2371 NEXT Jk 

2374 IF I- (9-Nsub) OR I-(IO-Nsub) THEN PRINT USING *4X,’*Kdl - ’".DD*;.Kdl 
2377 IF Kdl > 1 0 THEN 
2380 BEEP 
2383 BEEP 

2386 PRINT USING ™4X.“*Too much scattering in data - repeat data set“"“ 

2389 GOTO 1816 
2392 END IF 
2395 END IF 

2398 Emt ( I) -Sum/ ( 20-Kdl ) 

2401 NEXT I 

2404 IF Ihm-t THEN 

2407 OUTPUT 709:“AR AF20 AL20 VR5“ 
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2590 NEXT I 

2593 Tw-Tua/ ( 8-Idtc ) 

2596 END IF 
2599 Tld~T ( 9-Nsub ) 

2602 Tld2-T( 10-Nsub) 

2605 Tlda-(Tld+Tld2)+.S 
2608 Tv-T(ll-Nsub) 

26! I IF I tt<3 AND Ihm-0 THEN 

2614 T ld2--99 . 99 

2617 Tv- ( T(t 0 ) + T( I I ) )/2 

2620 END IF 

2623 Tsump-T( 12-Nsub ) 

2626 IF Ihm-0 THEN 2635 
2629 Tinlet-T( 13-Nsub) 

2632 Tout-T( 1 4-Nsub ) 

2635 IF Ihm-0 THEN 
2638 Amp-ABS(Ir) 

2641 Vo 1 t-ABS ( Vr ) *25 
2644 Q-Volt*Amp 
2647 END IF 

2650 IF Itt-0 AND Ihm-0 THEN 
2653 Kcu-FNKcu ( Tui ) 

2656 ELSE 

2659 Kcu-Kcua(Itt) 

2662 END IF 
2665! 

2668! FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
2671 IF Ihm-0 THEN Tu-Tw-Q*LOG( D2/D1) /(2*PI*Kcu*L) 

2674 IF Ilqv-0 THEN Tsat-TLda 

2677 IF IIqv-1 THEN Tsat- ( T Lda+Tv) * .5 

2680 IF Ilqv-2 THEN Tsat-Tv 

2633 IF Ihm-1 THEN 
2686 Tavg-TinLet 

2689 Grad-37 . 9853+ . 1 04388*Tava 

2692 Tdrop-ABS ( Emt ( 7 ) ) * I .E+6/M 0*Grad) 

2695 Tavgc-Tinlet-Tdrop* . 5 

2698 IF ABS ( Tavg-Tavgc ) > . 0 1 THEN 

270! Tavg- ( Tavg+Tavgc) * .5 

2704 GOTO 2689 

2707 END IF 

2710! 

2713! COMPUTE WATER PROPERTIES 
2716 IF Ihm-1 THEN 
2719 Ku-FNKw(Tavg) 

2722 Muwa-FNMuui ( Tavg ) 

2725 Cpu-FNCpw ( T avg ) 

2728 Pru-FNPrw( Tavg) 

2731 Rhow-FNRhou ( Tavg ) 

2734 Tuii-Tavg 
2737! 

2740! Compute MDOT 

2743 Mdot-3.9657E-3+Fms*(3.61955E-3-Fms*<8.32006E-6-Fms*< 1 .23688E-7-Fms*4.31897 

E-10) ) ) 

2746! Mdot-Mdo t*( 1 . 0365-T inlets 1 . 96644E-3-Tinlet+5 . 252E-6 ) ) / 1 . 0037 
2749 Kdt-0 

2752 Q-Mdot*Cpw*Tdrop 

2755 Lm td-Tdrop/LOG( ( Tiniet-Tsat) /( Tinlet-Tdrop-Tsat) ) 

2758 Uo-Q/ (PI*0o*L*Lmtd) 

2761 Ru-Qo+L0G( Do/Di ) / ( 2. +Kcu ) 

2764 Tw- Tsa t+Fr+Lmtd 
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2410 OUTPUT 709: “AS SA" 

2413 Sum-0 

2416 FOR Klc-I TO 20 

2419 ENTER 709 ;E 

2422 Sum-Sum+E 

2425 NEXT KV 

2428 Emf (71 -A8S (Sum) /20 

2431 END IF 

2434 IF Ihm-0 THEN 

2437 OUTPUT 709: "AR AF12 AL13 VR5" 

2440 FOR I- I TO 2 
2443 OUTPUT 709; “AS SA" 

2446 Sum-0 
2449 FOR Ji-1 TO 2 
2452 ENTER 709 :E 
2455 Sum-Sum+E 
2458 NEXT Ji 

2461 IF 1-1 THEN Vr-Sum/2 
2464 IF 1-2 THEN Ir-Sum/2 
2467 NEXT I 
2470 END IF 
2473 ELSE 

2476 IF Ihm-0 THEN ENTER ?File2:Bop,ToldS,Emt (♦) .Vr.Ir 
2479 IF Ihm-1 THEN ENTER 8FiIe2:8op.Told$,Emt (*) ,Fms 
2482 END IF 
2485! 

2488! CONVERT emt'S TO TEMP , VOLT, CURRENT 

2491 Twa-fl 

2494 . FOR 1-1 TO Ntc 

2497 IF ldtc>0 THEN 

2500 IF I-Ldtcl OR I-Ldtc2 THEN 

2503 T(I)— 89.99 

2506 GOTO 2536 

2509 ENO IF 

2512 END IF 

2515 IF Itt<4 AND Ihm-0 THEN 

2513 IF I>4 AND I<9 THEN 

2521 T(I)— 99.99 

2524 GOTO 2536 
2527 END IF 
2530 END IF 

2533 T( I ) -FNTvsv ( Emt ( I ) ) 

2536 NEXT I 

2539 IF Itt<4 THEN • 

2542 FOR 1-1 TO 4 

2545 IF I-Ldtct OR I-Ldtc2 THEN 

2548 Twa-Twa 

2551 ELSE 

2554 Twa-Twa + T( I ) 

2557 END IF 
2560 NEXT I 
2563 Tu»a-Twa/(4-Idtc) 

2566 ELSE 

2569 IF Ihm-1 THEN 2599 

2572 FOR 1-1 TO 8 

2575 IF I-Ldtcl OR I-Ldtc2 THEN 

2578 Twa-Tua 

2581 ELSE 

2534 Tua-Twa+TI I ) 

2587 END IF 
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2767 Vu-Mdo t/ ( Rhow-*PI*0 i * 2/4 ) 

2770 Reu’Rhou^Vu^-D L/Muua 

2773 H i*Ci*Kuj/D i+Rew * . 8*Prw * ( 1 /3 . ) * ( Maua/FNMuu ( Tu i ) ) * . 1 4 
2776 Tu ic* Tavg-Q/ ( PI*Do*L*H i ) 

2773 IF ABS( Tui-Tuic) > .01 THEN 

2782 Tui * ( Twi+Tuic ) * .5 

2785 GOTO 2773 

2788 END IF 

2791 Tui-( Tui+Tid ic ) * . 5 

2794 Ho- 1 / ( 1 /Uo-Do/ ( D i*H i ) -Ru ) 

2797 END IF 

2800 END IF 

2803 IF Ihm-1 THEN 

2806 Thetab-Q/(Ho*PI*Do*L) 

2809 Tw- Tsat + The tab 

2812 ELSE 

2815 Thetab-Tu-Tsat 

2813 END IF 

2821 IF The tab< 0 THEN 
2824 BEEP 

2827 INPUT *TUALL<TSAT (0-C0NTINUE. 1 -END) ‘.lev 

2830 IF Iev-0 THEN GOTO 1804 

2833 IF lev-! THEN 3130 

2836 END IF 

2829! 

2842! COMPUTE VARIOUS PROPERTIES 
2845 Tf i lm- ( Tu+Tsa t ) * . 5 
2848 Rho-FNRho(Tfilm) 

2851 Mu-FNMu( Tf llm) 

2854 K-FNK(Tfilm) 

2857 Cp-FNCp ( Tf i lm ) 

2860 8eta-FNBeta( Tf i lm) 

2863 Hfg-FNHfg(Tsat) 

2866 Ni-Mu/Rho 
2869 Alpha-K/(Rho*Cp) 

2872 Pr-Ni/Alpha 
2875 Psa t-FNPsa t ( Tsa t ) 

2873! 

2881! COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
2884! FOR UNENHANCED END(S) 

2887 Hbar-190 

2890 Fe* ( Hbar-*P/ ( Kcu*A ) ) * ,S*Lu 
2893 Tanh-FNTanh ( Fe) 

2896 The ta*The tao*Tanh/Fe 

2899 Xx- (9.31 *8eta*Thetab-*Do '3*Tanh/ ( Fe*N i+Alpha ) ) * . 166667 
2902 Yy-(1+( ,559/Pr ) ‘ ( 9/ 1 6 ) ) * (8/27) 

2905 Hbarc’K/Do* ( . 6+ . 387*Xx/Yy )* 2 

2908 IF ABS( (Hbar-Hbarc) /Hbarc) > . 001 THEN 

2911 Hbar-(Hbar+Hbarc)*.5 

2914 GOTO 2890 

2917 END IF 

2920! 

2923! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END ( S ) 

2926 Ql- ( Hbar-*P*Kcu*A ) * .5*Thetab*Tanh 
2929 Qc-0-2 *01 

2932 As-PI*D2*L 
2935! 

2938! COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT 

2941 Odp-Qc/As 

2944 H tube-Odp/The tab 
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2947 Csf • ( Cp*Thetab/Hf g ) / ( Qdp/ ( Mu*Hf g ) *( . 0 1 4/ ( 9. 3 1 *Rho ) * . 5 ) * ( 1 /3 . ) *Pr ' I . 7 ) 

2950! 

2953! RECORD TIME OF DATA TAKING 
2956 IF Im-0 THEN 
2959 OUTPUT 709; "TD* 

2962 ENTER 709:Toid$ 

2965 END IF 
2968! 

2971! OUTPUT DATA TO PRINTER 
2974 PRINTER IS 701 
2977 IF Iov-0 THEN 
2980 PRINT 

2983 PRINT USING *10X."*Data Set Number - “• ,DDD.2X, “ ’Built Oil X - .DD.D.SX, 1 
4A“ ; J.Bop. ToidS 
2986 IF Ihm-0 THEN 

2989 PRINT USING “10X,“*TC No: 1 2 3 4 5 6 7 

8* * * 

2992 PRINT USING "10X,**Temp : “ * , 3(1 X ,MDD . DD ) “ ; T( 1 ) , T< 2) , T( 3 ) . T ( 4) , T (5 ) , T (6 ) , T< 
7) .T(8) 

2995 PRINT USING ’10X.“* Tua Tliqd Tliqd2 Tvapr Psat Tsump“‘* 

2998 PRINT USING *1 0X,2(HDD. DD, 1 X) , 1 X ,MDD. DD , 1 X ,2(1 X.MDD .DD) ,2X,MDD. D* ; Tu , Tld.T 
Id2 , Tv, Psat, Tsumo 



3001 PRINT USING ‘tOX.— Thetab Htube Qdp“** 

3004 PRINT USING ’1 OX.MDD . 3D , 1 X .HZ. 3DE, 1 X .MZ.3DE" ; The tab .Htube , Qdp 
3007 ELSE 

3010 PRINT USING *10X,” Fins Vw Tsat Tinl Tdrop Thetab q Uo 

Ho 



3013* PRINT USING “1 0X,4(2D. DD, IX) .Z.3D , 1X.DD . DD . !X,3(MZ.3DE , IX) * ;Fms, Vw.Tsat , Ti 

n le t, Tdrop. The taD.Od 

3016 END IF 

3019 END IF 

3022 IF Iov-1 THEN 

3025 IF J-1 THEN 

3028 PRINT 

3031 IF Ihm-0 THEN 



3034 

3037 


PRINT USING "10X.“‘ RUN No 
ELSE 


Oi IX 


Tsat 


Htube 


Qdp 


Thetab" “ * 


3040 

3043 


PRINT USING ™ 1 0X . “ " FMS 
END IF 


0IL7. 


TSAT 


HTUBE 


QDP 


THETAB" 



3046 END IF 

3049 IF Ihm-0 THEN 

3052 PRINT USING “12X, 3D, 4X, DD. 2X. MDD. DD. 3UX, MZ. 3DE) *:J. Bop, Tsat, Htube, Qdp. The 
tab 

3055 ELSE 

3058 PRINT USING * 1 2X. 3D . 4X . DD ,2X .HDD . DD ,3 ( 1 X ,MZ. 3DE ) “ ;Fms,Bop . Tsa t .Htube .Qdp . T 
hetab 

3061 END IF 

3064 END IF 

3067 IF Im-0 THEN 

3070 BEEP 

3073 INPUT "OK TO STORE THIS DATA SET ( 1 -Y ,0-N ) ?* ,01c 
3076 END IF 

3079 IF Ok- I OR Im-I THEN J-J+1 
3082 IF Ok — 1 AND Im-0 THEN 

3085 IF Ihm-0 THEN OUTPUT »F i le I :Bop . To idS.Emf ( * > , Vr , Ir 
3088 IF Ihm-1 THEN OUTPUT SFi le 1 ;Bop . ToidS .Emf <*) ,Fms 

3091 END IF 

3094 IF Iuf-I THEN OUTPUT Rif ilelVu.Uo 

3097 IF Ire-1 THEN OUTPUT 8Ref i le ;Fms .Reu 

3100 IF (Im-I OR Ole-I) AND Iplot-I THEN OUTPUT »PIot:Qdp, Thetab 
3103 IF Im-0 THEN 
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3106 BEE? 

3109 INPUT "MILL THERE BE ANOTHER RUN < 1-Y.Q-N)?* ,Go_on 
3112 Nrun- J 

31 15 IF Go on-0 THEN 3130 
3118 IF Go"on<>0 THEN Repeat 
3121 ELSE ~ 

3124 IF J<Nrun+1 THEN Repeat 
3127 END IF 
3130 IF Ini-0 THEN 
3133 BEEP 

3136 PRINT USING “ 1 0X. “•NOTE: data runs were stored in file ““,10A*;J- 

1 , D2_f i leS 

3139 ASSIGN -File! TO * 

3142 OUTPUT AFiIe2:Nrun-1 
3145 ASSIGN SFilel TO D1_fiLe$ 

3148 ENTER SFile! :DateS .Ldtcl ,Ldtc2. Itt 
3151 OUTPUT SFiie2:DateS, Ldtcl ,Ldtc2,Itt 
3154 FOR 1-1 TO Nrun-1 
3157 IF Ihm-0 THEN 

3160 ENTER 4File1 ;3op , ToIdS. Emt (♦) .Vr.Ir 
3163 OUTPUT $F i le2 :3op , To IdS.Emf ( *) , Vr , Ir 
3166 ELSE 

3169 ENTER SFiiel 1 Boo. To IdS.Emf (*) ,Fms 
3172 OUTPUT $F i le2l Bop , ToIdS ,Emf ( * ) ,Fms 
3175 END IF 
31 78 NEXT I 

3181 ASSIGN JFilel TO * 

3184 PURGE “DUMMY” 

3187 END IF 
3190 BEEP 
3193 PRINT 

3196* IF Iplot-1 THEN PRINT USING “10X.“*N0TE: “*.Z Z.** X-Y pairs were stored in 
p iot data f i le ” “ , 1 
3199 ASSIGN ?Fiie2 TO * 

3202 ASSIGN «PIot TO * 

3205 IF Iuf-l THEN ASSIGN ?Ufile TO * 

3208 IF Ire-1 THEN ASSIGN SRefile TO * 

321 1 CALL Stats 
3214 BEEP 

3217 INPUT “LIKE TO PLOT DATA ( 1 «Y , 0-N) ?* ,01c 
3220 IF 01c -1 THEN CALL Plot 
3223 SUBEND 
3226! 

3229! CURVE FITS OF PROPERTY FUNCTIONS 
3232 DEF FNKcu(T) 

3235! 0FHC COPPER 250 TO 300 K 

3238 Tic -T*273 . 1 5 !C TO K 

3241 K-434- . 1 1 2* Tic 

3244 RETURN K 

3247 FNEND 

3250 DEF FNHu(T) 

3253! 170 TO 360 K CURVE FIT OF VISCOUSITY 

3256 Tlc-T*273 . 1 5 !C TO K 

3259 Mu-EXP ( -4. 4636 *( 1011 .47/Tlc ) )*1 .OE-3 

3262 RETURN (1u 

3265 FNEND 

3268 DEF FNCo(T) 

3271! 130 TO 400 K CURVE FIT OF Cp 
3274 Tlc-T+273. 15 !C TO K 

3277 Cp-. 40 138+1 .65007E-3*T1c + 1 . 5 1 494E-6*TIc '2-6.67853E-IO*T1c ‘3 
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3230 Cp-Cp*1 000 
3233 RETURN Cp 
3233 FNEND 
3289 DEF FNRho(T) 

3292 Tlc-T+273. tS !C TO K 

3295 X-1-(1.8*Tk/753.9S) !K TO R 

3298 Ro-36.32+6! . 1 4641 4*X* ( t /3) +1 8. 41 801S+X+1 7 . 476338+X* .5+1 . 1 19828+X*2 

3301 Ro-Ro/. 062423 

3304 RETURN Ro 

3307 FNEND 

3310 DEF FNPr(T) 

3313 Pr-FNCp(T)*FNMu(T)/FNK(T) 

3316 RETURN Pr 
3319 FNEND 
3322 DEF FNK ( T ) 

3325! T<360 K WITH T IN C 
3328 K" . 07 1 - . 000261 +T 
3331 RETURN K 
3334 FNEND 
3337 DEF FNTanh(X) 

3340 P-EXP(X) 

3343 0-1/P 

3346 Tanh-(P-Q) /(P+Q) 

3349 RETURN Tanh 

3352 FNEND 

3355 DEF FNTvsv(U) 

3353 COM /Cc/ C(7),Icai 

3361 T-C<0) 

3364 FOR 1-1 TO 7 
3367 T-T+C( I)*V" I 
3370 NEXT I 
3373 IF Ical-1 THEN 

3376 T-T-S. 7422934E-2+T*< 9. 0277043E-3-T+( -9.325991 7E-5 ) ) 

3379 ELSE 

3382 T-T+8. 626897E-2+T* ( 3 . 761 99E-3-T+S . 06392S9E-5 ) 

3385 END IF 
3388 RETURN T 
3391 FNEND 
3394 DEF FNBeta(T) 

3397 Rop-FNRho( T+. 1 ) 

3400 Rom-FNRho( T-. 1 ) 

3403 Beta--2/( Rop+Rom) * ( Rop-Rom) /.2 

3406 RETURN Beta 

3409 FNEND 

3412 DEF FNHfg(T) 

3415 Hfg-1 .3741344E+5-T*(3.3094361E+2+T J H .2165143) 

3418 RETURN Hfg 

3421 FNEND 

3424 DEF FNPsat(Tc) 

3427! 0 TO 80 deg F CURVE FIT OF Psat 
3430 Tf-1.8*Tc+32 

3433 Pa-5. 945525+ Tf*( . 1S3S2082+Tf+< 1 .4840963E-3+Tf+9.6150671E-S) ) 

3436 Pg-Pa-1 4 . 7 

3439 IF Pg>0 THEN ! +-PSIG. — in Hg 

3442 Psat-Pg 
3445 ELSE 

3448 Psa t-Pg+29 . 92/ 14.7 
3451 END IF 
3454 RETURN Psat 
3457 FNEND 
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3460 

3463 

3466 

3463 

3472 

3475 

3478 

3481 

3484 

3487 

3490 

3493 

3496 

3499 

3502 

3505 

3508 

3511 

3514 

3517 

3520 

3523 

3526 

3529 

3532 

3535 

3538 

3541 

3544 

3547 

3550 

3553 

3556 

3559 

3562 

3565 

3563 

3571 

3574 

3577 

3530 

3533 

3536 

3589 

3592 

3S95 

3598 

3601 

3604 

3607 

3610 

3613 

3616 

3619 

3622 

3625 

3628 

3631 

3634 

3637 



i 

i 



DEF FNHsmooth(X.3op . Isat) 

DIM A( 5) .8(5) .CIS) .DCS > 

DA TA .20526 , . 25322 ..31 9048 . .55322 , . 79909 . 1 . 002S8 
DATA .74515.. 72992 . . 73 1 39 , . 7 1 225 . . 634 72 . . 64 1 97 
DATA .41 092 . . 1 7726 . .2S 1 42 . .54806 . .81916, 1 .0845 
DATA .71403. .72913. .72S65. .696691 , .665867 , .61889 
READ A<*> ,B<*> ,C(*> .DU) 

IF Bop<6 THEN I-8op 
IF Bop-6 THEN 1-4 
IF Bop- 10 THEN 1-5 
IF Isat- 1 THEN 
Hs-EXP(A( I)+8( I)*L0G(X) ) 

ELSE 

Hs-EXP(C( I)+D( I) *L0G(X) ) 

END IF 
RETURN Hs 
FNEND 

DEF FNPoly(X) 

COM /Cpiy/ A(10,10),C(10).B(4) ,Nop, Iprnt , Opo ,Ilog,Ifn,I;oin,N;oin 
X 1 -X 

Poly-B(0) 

FOR 1-1 TO Nop 
IF Ilog-1 THEN Xi-LOG(X) 

Poly-Poly+B( I)*X1 * I 
NEXT I 

IF Ilog-1 THEN Poiy-EXP(Poiy ) 

RETURN Poly 
FNEND 
SUB Poly 

DIM R( 10) ,S< 10) .Sy( 12) .Sx< 12) .Xx( 100) .Yy( 100) 

COM /Cply/ A( 1 0 . 1 0 ) ,C( 1 0 ) . 3( 4) ,N , Iprnt , Opo , Ilog , If n, I joi n , N>o in 
COM /Xxyy/ Xp(25) , Yp(2S) 

FOR 1-0 TO 4 
B( I) -0 
NEXT I 
BEEP 

INPUT “SELECT (0-FILE. 1 -KEYBOARD .2-PROGRAM) *, Im 

Im- Im-t-1 

BEEP 

INPUT “ENTER NUMBER OF X-Y PAIRS". Mp 

IF Im- 1 THEN 

BEEP 

INPUT “ENTER DATA FILE NAME*,D_f ileS 
BEEP 

INPUT "LIKE TO EXCLUDE DATA PAIRS ( 1 -Y.O-N) . led 

IF led- 1 THEN 

BEEP 

INPUT "ENTER NUMBER OF PAIRS TO BE EXCLUDED" , Ipex 
END IF 

ASSIGN 9F i 1 e TO D_fileS 

ELSE 

BEEP 

INPUT "WANT TO CREATE A DATA FILE (1 - Y , 0-N ) ?* , Yes 

IF Yes-1 THEN 

BEEP 

INPUT "GIVE A NAME FOR DATA FILE" ,D_f i leS 

CREATE 8DAT D_fileS.S 

ASSIGN ?Fi ie TO D_file5 

END IF 

END IF 
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3640 

3643 

3646 

3649 

3652 

3655 

3658 

3661 

3664 

3667 

3670 

3673 

3676 

3679 

3632 

3635 

3688 

3691 

3694 

3697 

3700 

3703 

3706 

3709 

3712 

3715 

3718 

3721 

3724 

3727 

3730 

3733 

3736 

3739 

3742 

3745 

3748 

3751 

3754 

3757 

3760 

3763 

3766 

3769 

3772 

3775 

3773 

3781 

3784 

3787 

3790 

3793 

3796 

3799 

3802 

3805 

3808 

381 1 

3814 

3817 

3320 



BEEP 

INPUT "ENTER THE ORDER OF P0LrN011IAL , • ,N 

FOR 1-0 TO N*2 

Sy(I)-0 

Sx( I)-0 

NEXT I 

IF led- I AND Im-1 THEN 
FOR 1-1 TO Ipex 
ENTER »File;X,Y 
NEXT I 
END IF 

FOR 1-1 TO Np 
IF Im-1 THEN 

IF Opo-2 THEN ENTER 9File:X.Y 
IF 0po<2 THEN ENTER SFiU;Y,X 
IF Opo-1 THEN Y-Y/X 
IF I log— 1 THEN 
IF Opo-2 THEN Xt-X/Y 
X-LOG(X) 

IF Opo-2 THEN Y-LOG(Xt) 

IF 0po<2 THEN Y-LOG(Y) 

END IF 
END IF 

IF Im-2 THEN 
BEEP 

INPUT "ENTER NEXT X-Y PAIR“,X.Y 
IF Yes-1 THEN OUTPUT ®File;X.Y 
END IF 

IF Im<3 THEN 
Xx ( I ) -X 
Yy ( I) -Y 
ELSE 

X-Xp(I-l) ■ . 

Y-Yp( 1-1 ) 

END IF 
R ( 0 ) — Y 

Sy ( 0) -Sy ( 0 ) +Y 
S ( 1 ) -X 

Sx( 1 )-Sx< 1 )+X 
FOR J-1 TO N 
R ( J ) — R ( J- 1 )*X 
Sy ( J) -Sy ( J) +R( J) 

NEXT J 

FOR J-2 TO N*2 
SU)-S( J-1 )*X 
Sx(J)-Sx(J)+S( J) 

NEXT J 
NEXT I 

IF Yes-1 AND Im-2 THEN 
BEEP 

PRINT USING •12X,D0.’“ X-Y pairs >jere stored in file * " , 1 OA" ;Np ,D_f i leS 

END IF 

Sx ( 0) -No 

FOR 1-0 TO N 

C(I)-Sy(I) 

FOR J-0 TO N 
A( I . J) -Sx( I+J) 

NEXT J 
NEXT I 

FOR 1-0 TO N-1 
CALL Divide! I) 
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3823 CALL Subtract! 1*1 ) 

382S NEXT I 

3829 B(N)-C(N)/A(N.N) 

3832 FOR 1-0 TO N-1 
3835 B(N-1-I)»C(N-1-I) 

3838 FOR J-0 TO I 

3841 B(N-1-I)-B(N-1-I)-A(N-1-I.N-J)*B(N-J) 

3844 NEXT J 

3847 B(N-I-I)-B(N-1-I)/A(N-1-I,N-1-I) 

3850 NEXT I 

3853 SPRINTER IS 701 

3856 ! PRINT B(*) 

3859 SPRINTER IS 705 
3862 IF Iprnt-0 THEN 

3865 PRINT USING ’ 1 2X, " ’EXPONENT COEFFICIENT""" 

3868 FOR 1-0 TO N 

3871 PRINT USING "15X,DD.5X,MD.7DE";I,B( I) 

3874 NEXT I 
3877 PRINT " " 

3880 PRINT USING "12X,"*DATA POINT X Y Y( CALCULATED) DISCR 

EPANCY"’" 

3883 FOR 1-1 TO Np 
3886 Yc-B(O) 

3889 FOR J-l TO N 
3892 Yc-Yc+B(J)*Xx( I) * J 
3895 NEXT J 
3898 D-Yy ( I ) - Yc 

3901 PRINT USING "1 SX .3D . 4X , 4(M0 .5DE , 1 X ) " ; I ,Xx ( I) , Yy ( I ) , Yc ,D 

3904 NEXT I 

3907 END IF 

3910 ASSIGN »File TO * 

3913 SUBEND 
3916 SUB Divide(M) 

3919 COM /Cpiy/ A(10, 10), C(10),B(4),N, Iprnt, Opo.Ilog.Ifn.Ijoin.Njoin 

3922 FOR I-M TO N 

3925 Ao-A(I.M) 

3928 FOR J-M TO N 

3931 A(I.J)-A(I,J)/Ao 

3934 NEXT J 

3937 CCI)-C(I)/Ao 

3940 NEXT I 

3943 SUBEND 

3946 SUB Subtract(K) 

3949 COM /Cpiy/ A( 1 0 , 1 0 ) ,C( 1 0 ) ,B ( 4) ,N , Iprnt ,0po, I log, If n . I join .N join 

3952 FOR I-K TO N 

3955 FOR J-K-1 TO N 

3958 A(I.J)-A(K-1 .J)-A(I.J) 

3961 NEXT J 

3964 C( I)-CCK-1 ) -C ( I) 

3367 NEXT I 
3970 SUBEND 
3973 SUB PI in 

3976 COM /Cpiy/ A( 10. 10) ,C( 10) ,B<4) ,N, Iprnt ,0po . Ilog, Ifn. Ijoi n.N join 
3979 COM /Xxyy / Xx ( 25 ) , Yy ( 25 ) 

3982 PRINTER IS 705 
3985 8EEP 

3988 INPUT “WANT TO PLOT Uo vs Vu? ( 1 -Y , 0-N ) " . Iuo 
3991 IF Iuo-0 THEN 
3994 BEEP 

3997 INPUT “SELECT (0-h/h07. same tube. 1 -h (HF ) /h( sm) “,Irt 
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4000 

4003 

4006 

4009 

4012 

4015 

4013 

402! 

4024 

4027 

4030 

4033 

4036 

4039 

4042 

4045 

4048 

4051 

4054 

4057 

4080 

4063 

4066 

4069 

4072 

4075 

4079 

4081 

4084 

4087 

4090 

4093 

4096 

4099 

4102 

4105 

4108 

41 1 1 

4114 

4117 

4120 

4123 

4126 

4129 

4132 

4135 

4138 

4141 

4144 

4147 

4150 

4153 

4156 

4159 

4162 

4165 

4168 

41 71 

4174 

4177 

4180 



BEEP 

INPUT “SELECT h/h RATIO ( 1 -FILE , 0-CQMPIJTED 1'Mhrat 

IF Ihrat-Q THEN 

BEEP 

INPUT “WHICH Tsat ( 1 -6 . 7 , 0 — 2 .2) Msat 

END IF 

Xmin-Q 

Xmax- 1 0 

Xstep-2 

IF Irt-0 THEN 

Ynun-Q 

Ymax- 1 . 4 

Ystep- .2 

ELSE 

Y nun-0 

Ymax- 1 5 

Ystep-5 

END IF 

ELSE 

Qpo-2 

Ynun-0 

Ymax-12 

Ystep-3 

Xmin-Q 

Xmax-4 

Xs tep- 1 

END IF 

IF Ihrat-1 THEN 

Ynun-0 

Ymax- 1 8 

Ystep-3 

Xmin-0 

Xmax-9 

Xstep-2 

END IF 

BEEP 

PRINT ’• IN : SP 1 : IP 2300 .2200,3300,6800 ; “ 

PRINT “SC 0,100.0, 1 00 : TL 2,0;“ 

Sf x- 1 00/ ( Xmax-Xmi n ) 

Sfy- 1 00/ ( Ymax-Ymin ) 

PRINT “PU 0,0 PD" 

FOR Xa-Xnun TO Xmax STEP Xstep 

X- ( Xa-Xmin )*Sf x 

PRINT “PA" ;X, “ ,0; XT;" 

NEXT Xa 

PRINT "PA 100, 0:PU:“ 

PRINT “PU PA 0,0 PD" 

FOR Ya-Ymin TO Ymax STEP Ystep 
Y-( Ya-Ynun)*Sfy 
PRINT “PA 0,-;Y,-YT M 
NEXT Ya 

PRINT “PA 0,100 TL 0 2 m 

FOR Xa-Xnun TO Xmax STEP Xstep 

X- ( Xa-Xmin )*Sf x 

PRINT -PA-;X,",100: XT- 

NEXT Xa 

PRINT “PA 100.100 PU PA 100.0 PD“ 

FOR Ya-Ymin TO Ymax STEP Ystep 
Y-( Ya-Ymin) *Sfy 
PRINT -PD PA 100.\Y,-YT- 
NEXT Ya 
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A 1 83 PRINT "PA 100,100 PIT 
4136 PRINT "PA 0,-2 SR 1 .5,2” 

4183 FOR Xa-Xmin TO Xmax STEP Xstep 
4192 X-<Xa-Xmin)*Sfx 
4195 PRINT "PA" ;X, " ,0: " 

4198 IF Iuo-0 THEN PRINT “CP -2 , - 1 ;LB“ ; Xa : “ “ 

4201 IF Iuo-1 THEN PRINT "CP -1 . 5 , - 1 ;LB" ;Xa ; “ " 

4204 NEXT Xa 

4207 PRINT "PU PA 0,0" 

4210 FOR Ya-Ymxn TO Ymax STEP Ystep 
4213 IF ABS(Ya)<1.E-5 THEN Ya-0 
4216 Y- ( Ya-Ymin ) *Sfy 
4219 PRINT “PA 0,";Y,““ 

4222 IF Iuo-0 THEN PRINT “CP -4 f - . 25 ;LB" : Ya ; “ “ 

4225 IF Iuo-1 THEN PRINT “CP -3 , - . 25 ; LB“ ; Ya : “ “ 

4228 NEXT Ya 

4231 XladeiS-’Oil Percent - 
4234 IF Iuo-0 THEN 

4237 IF Irt-Q THEN 

4240 Y labels- "h/h0%" 

4243 ELSE 

4246 Y labels- “h/hsmooth" 

4249 END IF 

4252 PRINT “SR 1.5,2:PU PA 50,-10 CP* : -LEN( X labe IS ) /2: “0 :LB" ;X labels: “ * 

4255 PRINT “PA -11,50 CP 0, “ ;-LEN( YlabeiS) /2*5/6: “DI 0 , 1 ;LB“ ; YlaoelS; ““ 

4258 PRINT “CP 0,0“ 

4261 ELSE 

4264 PRINT “SP0;SP2“ 

4267 PRINT “SR 1.2.2.4;PU PA -8,35:DI 0,1;LBU:PR 1 t 0.5:LBo:PR -1.0.5;LB (kU/m 

m 



4270 PRINT "PR -t.0.5:SR 1,1.5:LB2:SR 1.5.2:PR ,5,.5;LB.:PR ,5.0:LBK)“ 
4273 PRINT “PA 42,H0;DI 1,0:LBV;PR . 4 , - 1 ; LBw ; PR 1 , ,S;LB(in/s) " 

4276 PRINT "SP0:SP1" 

4279 END IF 
4282 Ipn-0 
4285 BEEP 

4288 INPUT “WANT TO PLOT DATA FROM A FILE < 1 -Y , 0-N ) ?" ,Qkp 
4291 Icn-0 



4294 

4297 

4300 

4303! 

4306 

4309 

4312! 

4315 

4313 

4321 

4324 

4327 

4330 

4333 

4336 

4339 

4342 

4345 

4348 

4351 

4354 

4357 



IF OJcp-1 THEN 
BEEP 

INPUT “ENTER THE NAME OF THE DATA FILE“ ,D_f i leS 

IF Iuo-0 THEN 

BEEP 

INPUT "SELECT (0-LINEAR, 1 -L0G(X, Y) “ * Hog 
END IF 



ASSIGN 3File TO D_fileS 
BEE? 

INPUT “ENTER THE BEGINNING RUN NUMBER", Md 
BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" 
IF Iuo-0 AND Ihrat-Q THEN 



BEEP 

INPUT "ENTER DESIRED HEAT FLUX",Q 

END IF 

BEEP 

PRINTER IS 1 

PRINT USING *4X, “"Select a symbol:“““ 
PRINT USING " 4X , “ “ 1 Star 2 Pius sxgn M 
PRINT USING " 4 X , " “ 3 Circle 4 Square" - " 
PRINT USING "4X , " "5 RomDus ,, " ,, 



Npairs 
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4360 

4363 

4366 

4369 

4372 

4375 

4378 

4331 

4384 

4387 

4390 

4333 

4396 

4399 

4402 

4405 

4408 

441 1 

4414 

4417 

4420 

4423 

4426 

4429 

4432 

4435 

4438 

4441 

4444 

4447 

4450 

4453 

4456 

4459 

4462 

4465 

4463 

4471 

4474 

4477 

4480 

4483 

4486 

4489 

4492 

4495 

4498 

4501 

4504 

4507 

4510 

4513 

4516 

4519 

4522 

4525 

4528 

4531 

4534 

4537 

4540 



PRINT USING "4X." "6 
PRINT USING " 4X , " " 7 
INPUT Sym 
PRINTER IS 705 
PRINT "PU DI" 

IF Sym* 1 THEN PRINT 
IF Sym* 2 THEN PRINT 
IF Sym*3 THEN PRINT 
Nn-4 

IF Ilog- I THEN Nn* I 

IF Md>1 THEN 

FOR 1-1 TO Old-1) 



R ign t-3 Lde-up 
Up-3ide-down 



"SM*" 

"SMo" 



ENTER 9File;Xa,Ya 
NEXT I 
END IF 

IF Ihrat-0 THEN 
Ql-Q 

IF Ilog- 1 THEN Q-LOG(Q) 
END IF 

FOR 1-1 TO Npairs 
IF Iuo-0 AND Ihrat-0 THEN 



ENTER 3Fi le;Xa,B( *) 



Ya-B(O) 

FOR K-1 TO Nn 



tr langle" 
tr iangle” " 



Ya-Ya+B(K)*Q‘K 



NEXT K 
END IF 

IF Iuo-1 OR Ihrat-1 THEN 
ENTER SFile:Xa,Ya 
IF Iuo-I THEN Ya-Ya/1000 
END IF 

IF Iuo-0 AND Ihrat-0 THEN 
IF Ilog- 1 THEN Ya-EXP(Ya) 
IF Ilog-0 THEN Ya-Ql/Ya 
IF Irt-0 THEN 
IF Xa-0 THEN 



Yo-Ya 



Ya-I 

ELSE 

Ya-Ya/Yo 
END IF 
ELSE 



Hsm-FNHsmooth(Q,Xa • Isat) 

Ya-Ya/Hsm 

END IF 

END IF 

Xx ( 1-1 ) -Xa 

Yy ( 1-1 )-Ya 

X-< Xa-Xmin)*Sf x 



Y-( Ya-Ymm) *Sfy 
IF Sym>3 THEN PRINT 
IF Sym< 4 THEN PRINT 
PRINT "PA" ♦ X . Y . " * 

IF Sym>3 THEN PRINT 
IF Sym-4 THEN PRINT 
IF Sym-5 THEN PRINT 
IF Sym-6 THEN PRINT 
IF Sym-7 THEN PRINT 
NEXT I 
BEE? 

ASSIGN SFile TO *■ 



“Sht“ 

‘SR 1 .4,2.4" 

‘SR 1 .2.1.6" 

*UC2. 4. 99. 0,-8. -4. 0.0. 3. 4.0.;" 
"UC3 .0 , 99 , -3 . -6 , -3.6.3 . 6 .3. -6 ; ‘ 
"UCO .5.3 ,99.3, -8. -6.0 .3.3: ‘ 
•UCO, -5. 3, 99. -3, 3. 6, 0.-3. -8; ‘ 
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4543 

4546 

4549 

4552 

4555 

4553 

4561 

4564 

4567 

4570 

4573 

4576 

4579 

4582 

4535 

4588 

4591 

4594 

4597 

4600 

4603 

4606 

4609 

4612 

4615 

4613 

4621 

4624 

4627 

4630 

4633 

4636 

4639 

4642 

4645 

4648 

4651 

4654 

4657 

4660 

4663 

4666 

4669 

4672 

4675 

4678 

4681 

4634 

4637 

4690 

4693 

4696 

4699 

4702 

4705 

4703 

471 1 

4714 

4717 

4720 



END IF 

PRINT "PU sir 
BEEP 

INPUT -WANT TO PLOT A POLYNOMIAL (1-Y,0-N)?*,Qlcp 

IF Okp-1 THEN 

BEEP 

PRINTER IS 1 

PRINT USING 'MX.'* “Select line type:*"* 

PRINT USING “SX ♦ “ “0 Solid line"““ 

PRINT USING “6X.““1 Dashed“““ 

PRINT USING “6X.““2,,,5 Longer line - dash""* 

INPUT Ipn 
PRINTER IS 705 
BEEP 

INPUT “SELECT (0-LINEAR, t-L0G(X, Y) ) \Ilog 

Iprnt- 1 

CALL Poly 

IF Iuo-1 THEN 

BEEP 

INPUT “DESIRE TO SET X Lower and Upper Limit ( 1 -Y , O-N)?* * Ixlim 
IF Ixlim-0 THEN 
XI 1-0 
Xu 1*7 
END IF 

IF Ixlim-1 THEM 
BEEP 

INPUT “ENTER X Lower Limit\Xll 
BEEP 

INPUT “ENTER X Upper Limif.Xul 
END IF 
END IF 

FOR Xa-Xll TO Xul STEP Xstep/25 
Icn* Icn + 1 
Ya-FNPo ly ( Xa) 

IF Iuo-1 THEN Ya-Ya/1 000 

Y* ( Ya-Ymin ) *Sfy 

X- ( Xa-Xmin ) *Sf x 

IF Y<0 THEN Y-0 

IF Y > 1 00 THEN GOTO 4637 

Pu-0 

IF Ipn* ! THEN Idf-Icn MOD 2 
IF Ipn-2 THEN Idf-Icn MOD 4 
IF Ipn-3 THEN Idf-Icn MOD 3 

IF Ipn-4 THEN Idf-Icn MOD 16 

IF Ipn-5 THEN Idf-Icn MOD 32 

IF Idf-1 THEN Pu-1 

IF Pu-0 THEN PRINT “PA“ .X, Y , “PD“ 

IF Pu-1 THEN PRINT “PA“ .X. Y * “PU“ 

NEXT Xa 
PRINT “PU“ 

GOTO 4285 
END IF 
BEEP 

INPUT “WANT TO QUIT ( I - Y , 0-N ) ?“ , Iqui t 
IF Iquit-1 THEN 4711 
GOTO 4285 
PRINT “PU SP0" 

SUBEND 
SUB Stats 
PRINTER IS 701 
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4723 


J-0 




4726 


K-0 




4729 


BEEP 




4732 


IF Iplot-I THEN ASSIGN 9F ile TO P_f lie* 




4735 


BEEP 




4738 


INPUT "LAST RUN No? ( 0-QUIT) * ,Nn 




4741 


IF Nn-0 THEN 4849 




4744 


Nn-Nn- J 




4747 


Sx-0 




4750 


Sy *0 




4753 


Sz-0 




4756 


Sxs-0 




4759 


Sys-0 




4762 


Szs-Q 




4765 


FOR 1-1 TO Nn 




4768 


J-J + 1 




4771 


ENTER SFilecQ. T 




4774 


H-Q/T 




4777 


Sx-Sx+Q 




4780 


Sxs-Sxs+Q‘2 




4783 


Sy-Sy +T 




4786 


Sys-Sys+T‘2 




4789 


Sz-Sz+H 




4792 


Szs-Szs+H‘2 




4795 


NEXT I 




4798 


Qave-Sx/Nn 




4801 


Tave-Sy /Nn 




4804 


Have-Sz/Nn 




4807 


Sdevq-SQR( ABS ( (Nn*Sxs-Sx *2) /(Nn+( Nn- 1 ) ) ) ) 




4810 


Sdev t-SQR ( ABS ( (Nn+Sys-Sy‘2)/(Nn*(Nn-1 ) ) J ) 




4813 


Sdevh-SQR ( ABS ( ( Nn*S zs-Sz *2) / (Nn*( Nn - 1 ) ) ) ) 




4816 


Sh- 1 00*Sdevh/Have 




4819 


Sq- 1 00*Sdevq/Qave 




4822 


S t- 1 00*Sdevt/Tave 




4825 


IF K-1 THEN 4843 




4828 


PRINT 




4831 


PRINT USING -1 IX, "-DATA FILE : * * J 4A" ;Fi leS 




4834 


PRINT 




4837 

» 


PRINT USING MIX, ""RUN Htube SdevH Qdp SdevQ Thetab 


SdevT"" 


4840 


K-1 




4843 


PRINT USING M 1X,DD,2(2X,D.3DE, 1X,3D.2D> ,2X,DD.3D,1X.3D.2DMJ,Have 


.Sh.Qave 


t Sq ♦ Tave.St 




4846 


GOTO 4735 




4849 


ASSIGN i?F i 1 e 1 TO * 




4852 


PRINTER IS 1 




4855 


SUBEND 




4858 


SUB Coef 




4861 


CON /Cply/ A ( 10,10),C(t0),B(4).N, Iprnt .Opo, I log, If n, I join, N join 




4864 


BEEP 




4867 


INPUT "GIVE A NAME FOR CROSS-PLOT FILE",CpfS 




4870 


CREATE BDAT Cpf$,2 




4873 


ASSIGN ?File TO CpfS 




4876 


BEEP 




4879 


INPUT "SELECT (0-LINEAR . 1 -L0G(X, Y) ) “ , Ilog 




4882 


BEEP 




4885 


INPUT "ENTER OIL PERCENT (-1 -STOP) “ ,Bop 




4888 


IF Bop<0 THEN 4900 




4891 


CALL Poly 




4894 


OUTPUT '?File;Bop.B(*) 




4897 


GOTO 4882 
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4900 ASSIGN isFiie TO * 

4903 SUBENO 

4906 SUB Wi IsoniCf *Ci ) 

4909 COM /Wil/ 02, Di ,Do,L,Lu,Kcu 

4912 DIM Emf (12) 

491 S! WLISON PLOT SUBROUTINE DETERMINE CF AND Cl 

4913 BEEP 

4921 INPUT "ENTER DATA FILE NAME" ,Fi le$ 

4924 BEEP 

4927 PRINTER IS 1 

4930 PRINT USING MX, ""Select option:""" 

4933 PRINT USING "4X."" 0 Vary Cf and Ci"" n 

4936 PRINT USING "4X,"" 1 Fix Cf Vary Ci mmm 

4939 PRINT USING "4X,"" 2 Vary Cf Fix Ci""" 

4942 INPUT "ENTER OPTION" , Icf ix 

4945 PRINTER IS 701 

4948 IF Icfix-0 THEN 4960 

4951 IF Icf ix>0 THEN BEEP 

4954 IF Icf ix- 1 THEN INPUT "ENTER Cf",Csf 

4957 IF Icfix-2 THEN INPUT "ENTER CI",Ci 

4960 PRINTER IS 1 

4963 INPUT "Uant To Vary Coef f ? ( 1 -Y , 0-N ) " , Iccoef 
4966 IF Iccoef - 1 THEN INPUT "ENTER C0EFF",R 
4969 PRINTER IS 701 

4972 IF Icfix-0 OR Icfix-2 THEN Cfa-.004 

4975 IF Icfix-1 THEN Cfa-Csf 

4978 Cia-Ci 

4981 Xn- . 3 

4984 Fr-,3 

4987 Jj-0 

4990 Rr-3. 

4993 IF Iccoef- I THEN Rr-R 

4996 PRINTER IS 1 

4999 PRINT Do,Di,Kcu 

5002 ASSIGN SFile TO FileS 

5005 ENTER $Fi lelNrun, Dates, Ldtcl ,Ldtc2, Itt 

5008 Rw-Do*L0G( D 0 /D 1 ) / (2*Kcu) 

5011 Sx-0 

5014 Sy-0 

5017 Sxy-0 

5020 Sx2-0 

5023 Sy2-0 

5026 FOR 1-1 TO Nrun 

5029 ENTER SFilelBop, ToldS.Emf ( +) ,Fms 

50321 CONVERT EMF 'S TO TEMPERATURE 

5035 FOR J-l TO 5 

5038 T ( J ) -FNTvsv ( Emf ( J ) ) 

5041 NEXT J 

5044 Tsat-(T(1 )+T(2) )*.5 
5047 Tava-T (5) 

5050 Grad- 37 . 9853+ . 1 04388*Tava 
5053 Tdrop-Emf ( 7 ) *1 . E+6/ (10. *Grad ) 

50S6 Tavgc- T ( 5 ) - Tdrop* . 5 

5059 IF ABS ( Tavg-Tavgc J > . 0 1 THEN 

5062 Tavg- ( Tavg+Tavgc) * .5 

5065 GOTO 5050 

5068 END IF 

5071 ! 

5074! Compute properties of water 
5077 Kw-FNKw( Tavg ) 
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5080 Muua*FNMuu ( Tavg ) 

5083 Cpu*FNCpu ( Tavg ) 

5086 Pru*FNPru ( Tavo ) 

5089 Rhou-FNRhou ( Tavg) 

5092! 

5095! Compute properties of Freon-114 

5098 Lmtd* Tdrop /LOG ( ( T(5 ) -Tsa t ) / (T(5) -Tdrop- Tsa t ) ) 

510t IF Jj-0 THEN 
5104 Tu*Tsat+Fr*Lmtd 
5107 The tab*Tu-Tsat 
5110 Jj-1 
5113 END IF 
5116 Tf - ( Tu+Tsa t ) * .5 
51 19 Rho*FNRho< Tf ) 

5122 Mu-FNMu(Tf) 

5125 K-FNK(Tf) 

5128 Cp-FNCp(Tf) 

5131 Beta-FNBeta( Tf ) 

5134 Hfg-FNHfg(Tsat) 

5137 Ni*Mu/Rho 
5140 A lpha*K/ ( Rho+Cp ) 

5143 Pr-Ni/Alpha 
5146! 

5149! Analysis 
5152! COMPUTE MDOT 
5155 A«PI*(0o*2-Di‘2)/4 

5158 P-P I*Do 

5161 Mdo t*3 . 9657E-3+Fms* ( 3 . 6 1 955E-3-Fms-*( 8 . 32006E-6-Fms* ( 1 .23638E-7-Fms*4.31 897 
E— 10))) 

5164 Q -Mdo t*Cpu* Tdrop 

5167! COMPUTE NATURAL -CONVECTIVE HEAT-TRANSFER COEFFICIENT 
5170! FOR UNENHANCED EN0(S) 

5173 Hbar- 1 90 

5176 Fe*(Hbar+P/(Kcu*A) ) *.5*Lu 
5179 Tanh»FNTanh ( Fe ) 

5182 Theta*Thetab*Tanh/Fe 

5185 Xx- ( 9 . 81 *Beta*The tab ♦Do * 3*Tanh/ ( Fe*N i*A Ipna ) ) * . 1 66667 
5188 Yy-(1+( .559/Pr) “ (9/16) ) * (8/27) 

5191 Hbarc~K/Do*< . 6+ . 387*Xx/Yy ) ‘2 

5194 IF ABS( (Hbar-Hbarc) /Hbaj* ) > .001 THEN 

5197 Hbar-(Hbar+Hbarc) * .5 

5200 GOTO 5176 

5203 END IF 

5206! 

5209! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 

5212 Q1 -(Hbar'*P*Kcu*A) * .5*Thetab*Tanh 

5215 Qc-Q-2*Q1 

5218 As-PI*D2*L 

5221 ! COMPUTE ACTUAL HEAT FLUX 

5224 Qdp-Qc/As 

5227 IF Ictix-0 OR Icfix>1 THEN Csf-1 /Cf * ( 1 . /Rr ) 

5230 Thetab-Csf/Cp*Hfg*(Qdp/(Mu*Hfg)*( .014/(9,31 *Rho) ) * .5 ) * ( 1 /Rr) *Pr* 1 . 7 

5233 Ho=“Qdp/Thetab 

5236 Qmega*Ho/Cf 

5239 Uo-Q/(PI*0o*L+Lmtd) 

5242 Vw-Mdot/ ( Rhoui*PI*D 1*2/4) 

5245 Reu*Rhow-*Vu*Di/Muua 
5248 Tu i~Tu+Q*Ru/ ( PI*Do*L ) 

525 1 Ga ma *Ku/ D i*Reu * . 3*Pru * ( 1 /3. ) *( Muua/FNMuu ( Tui ) ) ‘ . 1 4 

5254! PRINTER IS 1 

5257 Yu* ( 1 . /Uo-Ru ) *Gmega 
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5260 

5283 

5266 

5263 

5272 

5275 

5279 

5281 

5284 

5287 

5290 

5293 

5296 

5299 

5302 

5305 

5308 

531 1 

5314 

5317 

5320 

5323 

5326 

5329 

5332 

5335 

5338 

5341 

5344 

5347 

5350 

5353 

5356 

5359 

5362 

5365 

5368 

5371 

5374 

5377 

5380 

5383 

5386 

5389 

5392 

5395 

5398 

5401 

5404 

5407 

5410 

5413 

5416 

5419 

5422 

5425 

5423 

5431 

5434 

5437 






Xui-Omega*Oo/ ( Gama*0 1 ) 

Sx-Sx+Xw 
Sy “Sy ♦ Yw 

Sxy-Sxy + Yw*Xuj 

Sx2"Sx2+Xtd*Xw " * * 

Sy2-Sy2+Yw*Yw 
NEXT I 

ASSIGN aFiie TO ♦ 

M* ( Sx*Sy-Nrun+Sxy ) / ( Sx*Sx-Nrun*Sx2 ) 

C* ( Sy-Sx*f1 ) / Nrun 

IF Icfix-0 OR Icfix-3 OR Icfix-4 THEN 

Cic-1/ft 

Cfc-1/C 

END IF 

IF Icfix-1 THEN 
Cic-1/N 
Cfc-Cf 
END IF 

IF Icfix-2 THEN 
Cic-Ci 
Cfc-1/C 
END IF 

IF ABS( (Ci-Cic)/Cic)>.001 OR ABS( (Cf-Cfc) /Cfc) >.001 THEN 
Ci*(Ci+Cic)*.5 
Cf *( Cf +Cfc ) * > 5 
PRINTER IS 1 

PRINT USING "2X,"“ Csf - •• ,f1Z.3DE,2X, "" Ci - " " ,HZ . 3DE" :Csf , Ci 

PRINTER IS 701 

GOTO 5002 

END IF 

PRINT 

PRINTER IS 701 

PRINT USING "23X , 11 " Cf Ci - "" 

PRINT USING "3X, “"ASSUMED "" ,MZ73DE.3X .MZ.3DE" iCfa.Cia 

PRINT USING "3X. “"CALCULATED "" ,MZ.3DE.3X,MZ.3DE" ;Csf ,Ci 
PRINT 

Sum2“Sy2-2*M*Sxy-2*C*Sy +M*2*Sx2 + 2*M*C*Sx+Mrun*C ‘2 
PRINT USING “ 1 OX , " “Sum of Squares - "" ,Z.3DE“ :Suro2 
PRINT USING "10X. “"Coefficient - *" .D.DDD" :Rr 
SUBEND 
DEF FNMuw(T) 

A-247 .3/ ( T+133. 15) 

Mu*2. 4E-5*1 0 ‘A 
RETURN Mu 
FNEND 

DEF FNCpw(T) 

Cpw-4.21 1 20858-T-M 2 .26826E-3-T*( 4 . 42381 E-5 +2. 7 1 428E-7*T) ) 

RETURN Cpu*1 000 
FNEND 

DEF FNRhoid(T) 

Ro-999 .52946+TM . 0 1 269-T*< 5 . 4825 1 3E-3-T*1 .234147E-5) ) 

RETURN Ro 
FNEND 

DEF FNPrw(T) 

P ru-FNCpui ( T ) *FNMu«u ( T ) /FNKw ( T ) 

RETURN Prm 
FNEND 

DEF FNKuj(T) 

X*( T+273. 15 ) 1212. 15 

Kw— ,92247+X*r2.3395-X*( 1 ,S007-X*( .52577- . 07344*X ) ) ) 
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5440 

5443 

5446 

5449 

5452 

5455 

5453 

5461 

5464 

5467 

5470 

5473 

5476 

5479 

5482 

5485 

5488 

5491 

5494 

5497 

5500 

5503 

5506 

5509 

5512 

5515 

5518 

5521 

5524 

5527 

5530 

5533 

5536 

5539 

5542 

5545 

5548 

5551 

5554 

5557 

5560 

5563 

5566 

5569 

5572 

5575 

5578 

5581 

5534 

5587 

5590 

5593 

5596 

5599 

5602 

5605 

5608 

5611 

5614 

5617 

5620 



RETURN Kw 
FNEND 
SUB Plot 

COM /Cply/ 8(10,10). C(I0),B(4) .Nop , Iprnt , Opo .Ilog.Ifn.Ijoxn.Njoin 
DIM Bs(3) 

INTEGER lx 
PRINTER IS 1 
BEEP 

INPUT “SELECT HEATING MODE (0-ELECTRIC, 1 -WATER )“, Ihm 
Idv-0 

BEEP 

INPUT “LIKE DEFAULT VALUES FOR PLOT ( 1 -Y , 0-N) ?* , Idv 

0po-0 

BEEP 

PRINT USING “4X. “"Select 0ptxon:“*“ 

PRINT USING *6X,"“0 q versus delta-T""‘ 

PRINT USING “SX.“*1 h versus delta-T"““ 

PRINT USING ■6X,"*2 h versus q“““ 

INPUT Opo 
BEEP 

INPUT “SELECT UNITS (0-SI, 1 -ENGLISH) “, Iun 
PRINTER IS 705 
IF IdvOI THEN 
BEEP 

INPUT "ENTER NUMBER OF CYCLES FOR X-AXIS*, Cx 
BEEP 

INPUT “ENTER NUMBER OF CYCLES FOR Y-AXIS“,Cy 
BEEP 

INPUT “ENTER MIN X-VALUE (MULTIPLE OF 10)".Xmxn 
BEEP 

INPUT "ENTER MIN Y-VALUE (MULTIPLE OF 10)*.Ymxn 
ELSE 

IF Opo-O THEN 

IF Ihm- 1 THEN 

Cy- I 

Cx-2 

Xmin-1 

Ymn-10000 

ELSE 

Cy-3 

Cx-3 

Xnun- . 1 

Ymin- 1 00 

END IF 

END IF 

IF Opo- 1 THEN 

IF Ihm- I THEN 

Cy-2 

Cx-2 

Xm x n- . 1 

Ynun-100 

ELSE 

Cy-3 

Cx-2 

Xmxn- I 000 
Ymxn-1 00 
END IF 
END IF 

IF Opo-2 THEM 
IF Ihm-1 THEN 
Cy-2 
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5523 

5626 

5623 

5632 

5635 

5633 

5641 

5644 

5647 

5650 

5653 

5656 

5659 

5662 

5665 

5663 

5671 

5674 

5677 

5630 

5633 

5636 

5689 

5692 

5695 

5693 

5701 

5704 

5707 

5710 

5713 

5716 

5719 

5722 

5725 

5728 

5731 

5734 

5737 

5740 

5743 

5746 

5749 

5752 

5755 

5753 

5761 

5764 

5767 

5770 

5773 

5776 

5779 

5782 

5785 

5738 

5791 

5794 

5797 

5800 






Cx - 1 

Xm i n- 1 0000 
Ymin-1 00 
ELSE 
Cy-3 
Cx-3 
Xnun- . 1 
Ymi n- 1 00 
END IF 
END IF 
END IF 
3EEP 

PRINT -IN:SP1:IP 2200,2200.3300,6300;- 
PRINT "SC 0,100,0,100; TL 2,0;“ 

Sfx-100/Cx 

Sf^-100/Cy 

INPUT "WANT TO BY-PASS CAGE? ( 1 - Y , 0-N) * , Ibyp 
IF Ibyp* 1 THEN 6049 
PRINT “PU 0,0 PD- 
Nn-9 

FOR 1-1 TO Cx + 1 
Xat-Xmin*10- (1-1 ) 

IF I-Cx+1 THEN Nn-1 
FOR J-1 TO Nn 

IF J-1 THEN PRINT "TL 2 0" 

IF J-2 THEN PRINT "TL 1 0“ 

Xa-Xat*J 

X-LGT ( Xa/Xnun) *Sf x 
PRINT "PA";X,",0; XT;" 

NEXT J 
NEXT I 

PRINT "PA 100,0;PU; “ 

PRINT "PU PA 0,0 PD" 

Nn-9 • 

FOR 1-1 TO Cy+1 
Yat-Ymin*10* (1-1 ) 

IF I-Cy+1 THEN Nn-1 
FOR J-1 TO Nn 

IF J-1 THEN PRINT "TL 2 0" 

IF J-2 THEN PRINT "TL 1 0" 

Ya-Yat+J 

Y-LGT ( Ya/Ymin) +Sf y 
PRINT "PA 0 , ; Y, " YT" 

NEXT J 
NEXT I 

PRINT "PA 0,100 TL 0 2" 

Nn-9 

FOR 1-1 TO Cx + 1 
Xat-Xmin*10MI-1) 

IF I-Cx+I THEN Nn-1 
FOR J-1 TO Nn 

IF J-1 THEN PRINT "TL 0 2" 

IF J>1 THEN PRINT “TL 0 1" 

Xa-Xat+J 

X-LGT (Xa/Xmin) *Sf x 
PRINT "PA" ;X, " , 100: XT- 
NEXT J 
NEXT I 

PRINT "PA 100,100 PU PA 100,0 PD - 
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5303 Nn-9 

5306 FOR 1-1 TO Cy+I 
5809 Yat-Ynun*10‘ ( 1-1) 

5812 IF I-Cy+1 THEN Nn-I 

5815 FOR J-l TO Nn 

5818 IF J-1 THEN PRINT “TL 0 2* 

5821 IF J>1 THEN PRINT " TL 0 1* 

5824 Ya-Yat+J 

5827 Y-LGT (Ya/Ynun)*Sfy 

5830 PRINT “PD PR 1 00 , “ , Y , “ YT" 

5833 NEXT J 
5836 NEXT I 

5839 PRINT “PA 100,100 PU* 

5842 PRINT “PA 0,-2 SR 1.5,2" 

5845 Ii-LGT(Xmin) 

5848 FOR 1-1 TO Cx+1 
5851 Xa-Xmin+1 O' ( 1-1 ) 

5854 X-LGT (Xa/Xmin)*Sfx 
5357 PRINT *PA“;X.\0:“ 

5860 IF Ii >— 0 THEN PRINT “CP -2, -2:IB1 0:PR -2,2*,LB“ : Ii J ** 

5863 IF li<0 THEN PRINT “CP -2, -2;LB1 0 :PR 0 ,2:LB“ ; Ii : * " 

5866 Ii-Ii+1 

5869 NEXT I 

5872 PRINT “PU PA 0.0* 

5875 Ii-LGT(Ymin) 

S878 Y10-10 

5881 FOR 1-1 TO Cy+1 

5884 Ya-Ymin+1 0 * ( 1-1 ) 

5887 Y-LGT ( Ya/Yn>in)*Sfy 
S390 PRINT "PA 0,“;Y."“ 

5893 PRINT "CP -4,-.25:LB1 0;PR -2,2:LB* :Ii 1 "* 

5896 Ii-Ii+I 
5899 NEXT I 
5902 BEEP 

5905 INPUT "UANT USE DEFAULT LABELS ( 1 -Y . 0-N)?“ , Idl 
5908 IF Idi<>1 THEN 
5911 BEEP 

5914 INPUT “ENTER X-LABEL* .XlabelS 
591 7 BEEP 

5920 INPUT “ENTER Y-LABEL* , YlabelS 

5923 END IF 

5926 IF Opo<2 THEN 

5929 PRINT “SR 1 ,2:PU PA 40,-14:“ 

5932 PRINT “LB(T;PR -1.6,3 PD PR 1.2.0 PU:PR .5,-4;LBuo:PR .5,1:“ 

5935 PRINT “LB-T:PR ,5,-1 :LBsat:PR .5,1;" 

5938 IF Iun-0 THEN 
5941 PRINT “LB) (K)“ 

5944 ELSE 

5947 PRINT “LB) (F)“ 

5950 END IF 
5953 END IF 
5956 IF Opo- 2 THEN 
S959 IF Iun-0 THEN 

5962 PRINT “SR 1.5,2;PU PA 40,-14;LBq (U/m;SR 1,I.5;PR 0.5,1;LB2:SR 1.5,2:PR 
o .5, -i :LB)“ 

5965 ELSE 

5963 PRINT “SR 1.5.2:PU PA 34,-14:LBq (Btu/hr:PR ,5..5:LB.:PR .5, -.5;* 

5971 PRINT "LBf t;PR ,5,1;SR 1,1.5;LB2;SR 1 . 5 , 2 : PR ,S.-1;LB);“ 

5974 END IF 
5977 END IF 
5980 IF Opo-O THEM 
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5383 IF Iun-0 THEN 

5386 PRINT "SR 1.S.2:PU'PA -12.40:01 0.1:L3q lW/m:PR -l.0.5:SR 1.I.5:LB2:SR t 
.5. 2: PR 1 . . 5 : LB ) " 

5983 ELSE 

5992 PRINT "SR 1.5.2:PU PA -12.32:01 0.1:L3q <Btu/hr:PR -.5. .5;L3. ;PR .5. .5:" 
5995 PRINT "LBf t :SR 1,1.S;PR -1,.5;L32:PR 1,.5;SR 1.5.2:LB)" 

5998 END IF 

6001 END IF 

6004 IF 0po>0 THEN 

6007 IF Iun-O THEN 

6010! PRINT "SR 1 .S,2;PU PA - 1 2 .38 : DI 0.1:LBh (W/m:PR -1..5:SR 1.1.5 ;LB2 ;SR 1. 
5.2: PR . 5 . .5 : * 

6013 PRINT "SR 1.2,2.4:PU PA -12.37:01 0,1;LBh;PR 1.0.5;L3o:PR -1,0.5:LB (W/m 

m 

6016 PRINT "PR - 1 . . S : SR 1.1.S:LB2:SR 1.S.2:PR .5..S:LB.:PR . 5 . 0 : LBK ) " 

6019 ELSE 

6022 PRINT "SR 1.5.2:PU PA -12.28:01 0.l:LBh (Btu/hr:PR -.5. .S:LB. :PR .5. .5:* 

6025 PRINT "LBf t:PR -1,.5:SR 1.t.5:LB2:SR 1.5.2:PR .5,.S:LB.:PR .S,.S:L3F)“ 

6028 END IF 

6031 END IF 

6034 IF Idl-0 THEN 

6037 PRINT "SR 1.5.2:PU PA 50.-16 CP" : -LEN( X labels ) /2: "0 :LB" : X labels : " 

6040 PRINT "PA -14.50 CP 0 . " : -LEN< Y laDelS ) /2*5/6 : "DI 0, 1 ;L3" ; Y laDelS ; " " 

6043 PRINT "CP 0,0 DI" 

6046 END IF 

6049 Ipn-0 

6052 Xll-I.E+6 

6055 Xul— 1.E+6 

6058 Icn-0 

6061 Ifn-0 

6064 Ijoin-1 

6067 BEEP 

6070 INPUT "WANT TO PLOT DATA FROM A FILE ( 1 -Y.O-N)?* ,0V 
6073 IF Ok -1 THEN 
6076 BEEP 

6079 INPUT "ENTER THE NAME OF THE DATA FILE" ,D_f i leS 
6082 ASSIGN »File TO D_fileS 
6085 BEEP 

6088 BEEP 

6091 INPUT "ENTER THE BEGINNING RUN NUMBER". Md 
6094 BEEP 

6097 INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" .Npairs 
6100! BEE? 

6103! INPUT "CONNECT DATA UITH LINE ( 1 - Y ,0-N) ?" , Ici 

6106 BEEP 

6109 PRINTER IS 1 

6112 PRINT USING "4X. “"Select a symool:""" 

6115 PRINT USING "6X.""1 Star 2 Plus sian’"" 

6118 PRINT USING "6X.""3 Circle 4 Square""" 

6121 PRINT USING "6X.""5 Ronous""" 

6124 PRINT USING "6X.“"S Right-side-up triangle*"" 

6127 PRINT USING "6X."*7 Up-side-down triangle""" 

6130 INPUT Sym 

6133 PRINTER IS 705 
6136 PRINT "PU DI“ 

6139 IF Sym-1 THEN PRINT *SM*“ 

6142 IF Sym - 2 THEN PRINT *SM+“ 

6145 IF Sym-3 THEN PRINT "SMo" 

6148 IF Md> 1 THEN 
6151 FOR 1-1 TO (Md-1) 
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6154 ENTER ?Fiie: Ya.Xa 

6157 NEXT I 

6160 END IF 

6163 FOR 1-1 TO Npairs 

6166 ENTER SFileJ Ya.Xa 

6169 IF 1-1 THEN Ql-Ya 

6172 IF I-Npairs THEN Q2-Ya 

6175 IF Opo-1 THEN Ya-Ya/Xa 

6178 IF Opo-2 THEN 

6181 Q-Ya 

6184 Ya-Ya/Xa 

6187 Xa-Q 

6190 END IF 

6193 IF Xa<Xl 1 THEN Xll-Xa 
6196 IF Xa>Xu 1 THEN Xul-Xa 
6199 IF Iun-1 THEN 
6202 IF 0 po< 2 THEN Xa-Xa*1 .8 
6205 IF 0po>0 THEN Ya-Ya+.l76t 
6208 IF 0 po- 0 THEN Ya-Ya*.317 
6211 IF Opo-2 THEN Xa-Xa*.317 
6214 END IF 

6217 X«LGT(Xa/Xmin)*Sfx 
6220 Y-LGT (Ya/Ymin)*Sfy 
6223 K;-0 

6226 CALL Symb(X.Y,Syn»,Icl,Kj) 

6229 GOTO 6268 

6232 IF Sym;3 THEN PRINT "SM* 

6235 IF Sym<4 THEN PRINT “SR 1.4, 2. 4* 

6238 IF Icl-O THEN 
6241 PRINT “Pfl“.X,Y,“* 

6244 ELSE 

6247 PRINT “PA* ,X, Y, “PD" 

6250 END IF 

6253 IF Sym>3 THEN PRINT “SR 1.2. 1.6* 

6256 IF Syn-4 THEN PRINT “UC2 . 4 , 99 . 0 . -8 . -4 , 0 .0 . 8 . 4 . 0 : “ 
6259 IF Sym-S THEN PRINT “UC3.0.99.-3.-6. -3.6.3.S.3, -6; 
6262 IF Sym-6 THEN PRINT “UCO .5 .3. 99.3 . -8 . -6 . 0 .3 ,3: “ 
6265 IF Sym-7 THEN PRINT "UCO ,-S .3,99 . -3, 3 .6. 0 , -3. -8: “ 
6268 NEXT I 
6271 PRINT "PU“ 

6274 BEEP 

6277 INPUT "WANT TO LABEL? ( 1 - Y , 0-N ) ’ , I lab 
6280 IF Ilab-1 THEN 
6283 PRINT “SP0:SP2* 

6286 BEEP 

6289 IF Klab-0 THEN 
6292 X lab-5 
6295 Y lab-85 



6298 


INPUT "ENTER INITIAL X.Y LOCATIONS - 


tXlab . 


Ylab 


6301 


Xtt-Xlab-S 






6304 


Ytt-Y lab+8 






6307 


PRINT “SR 1,1.5" 






6310 


BEEP 






6313 


INPUT “SELECT LABEL TYPE (Q-FILE NAME ♦ 1- 


TUBE 


6316 


IF Ifntt-0 THEN 






6319 


PRINT “SM:PA" ,Xtt. Ytt, *LB 7. 


Heat 


File 


6322 


Ytt-Ytt-3 






632S 


PRINT “SM:PA“ ,Xtt, Ytt. “LB 7. 


Heat 


File 


6328 


ELSE 






6331 


PRINT “SM:PA“ .Xtt . Ytt , "LB 7. 


Heat 


Tube 



, If ntt 
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Oil Flux Type 



6334 Ytt-Ytt-3 
S337 PRINT *PA*.Xtt.Ytt.“LB 
6340 END IF 

6343 IF Sym-I THEN PRINT "SM*“ 

6346 IF Sym-2 THEN PRINT *SM + “ 

6349 IF Sym-3 THEN PRINT ,, SI1o'‘ 

6352 K Lab" 1 
6355 END IF 
6358 Kj-1 
6361 CALL SyroblXlab. Ylab.Syro.Icl.Kj) 

6364 PRINT "SR 1,1.S:SM“ 

6367 IF Syn»<4 THEN PRINT "PR 2,0* 

6370 BEEP 

6373 INPUT "ENTER BOP’.Bop 
6376 IF Ifntt-1 THEM 
6379 BEEP 

6382 INPUT "ENTER TUBE TYPE* . T typeS 
6385 END IF 

6388 IF 8op< 1 0 THEN PRINT "PR 2.0 :LB" :Bop ; * * 

6391 IF Bop>9 THEN PRINT "PR .S,0:L8" ;Bop;** 

6394 Ihf-0 

6397 IF 01 >02 THEN Ihf-I 

6400 IF Ihf-0 THEN PRINT *PR 4,0:L8Inc“ 

6403 IF Ihf-I THEN PRINT "PR 4,0;LBDec" 

6406 IF Ifntt-0 THEN 

6409 PRINT "PR 2,0;LB*;D_f ileS; ** 

6412 ELSE 

6415 PRINT *PR 2.0:L3*;TtypeS:** 

6418 END IF 

6421 PRINT "SPOlSPI ;SR 1.5.2* 

6424 YLao-YIab-S 
6427 END IF 
6430 BEEP 

6433 ASSIGN SFile TO ♦ 

6436 XI L —X 11/1 .2 
6439 Xul-Xul*1 .2 
6442! GOTO 3040 
6445 END IF 
6448 PRINT "PU SM* 

6451 BEEP 

6454 INPUT "WANT TO PLOT A POLYNOMIAL ( 1 -Y.O-N)?* ,Go_on 

6457 IF Go_on* 1 THEN 

6460 BEEP 

6463 PRINTER IS 1 

6466 PRINT USING "4X . “ "Select line type:*** 

6469 PRINT USING "SX.“"0 Solid line"*" 

6472 PRINT USING “6X."*1 Dasned""* 

6475 PRINT USING "6X,**2,,,S Longer Line - dash*** 

6478 INPUT Ipn 
6481 PRINTER IS 705 
6484 BEEP 

6487 INPUT "SELECT ( 0-LlN . 1 -LOG( X . Y 1 ) * , I log 

6490 Iprnt-1 

6493 CALL Poly 

6496 IF Ifn-0 THEN 

6499 Njoin-I 

6502 BEEP 

6505 INPUT "ENTER NUMBER OF FILES TO JOIN?*,Njoin 
6508 END IF 
6511 I jo i n-0 

6514 IF Ifn<N;oin THEN Ijoin-I 
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6517 

6520 

6523 

6526 

6523 

6532 

6535 

6538 

6541 

6544 

6547 

6550 

6553 

6556 

6559 

6562 

6565 

6568 

6571 

6574 

6577 

6580 

6583 

6586 

6589 

6592 

6595 

6598 

6601 

6604 

6607 

6610 

6613 

6616 

6619 

6622 

6625 

6628 

6631 

6634 

6637 

6640 

6643 

6646 

6649 

6652 

6655 

6658 

666 1 

6664 

6667 

6670 

6673 

6676 

6679 

6682 

6685 

6633 

6691 

6694 



IF Ifn-0 OR I;otn-1 THEN 

FOR Ij-0 TO 3 

Bs( I j ) -Bs( I; ) +B( I j ) 

NEXT Ij 
If n-If n+1 
ENO IF 

IF N;om-Ifn THEN 

FOR I j-0 TO 3 

B( I;)-8s( I j)/Njoin 

Bs(I;)-0 

NEXT Ij 

Ifn-0 

ELSE 

GOTO 6067 
END IF 
BEEP 

INPUT “ENTER Y LOWER AND UPPER LIMITS" , Y 1 1 , Yul 
FOR Xx-0 TO Cx STEP Cx/200 
Xa*Xmin*1 0 'Xx 

IF Xa<XI I OR Xa>Xul THEN 6655 

Icn-Icn+1 

Pu-0 

IF Ipn-I THEN Idf-Icn MOD 2 

IF Ipn-2 THEN Idf-Icn MOD 4 

IF Ipn-3 THEN Idf-Icn MOD 8 

IF Ipn-4 THEN Idf-Icn MOD 16 

IF Ipn-5 THEN Idf-Icn MOD 28 

IF Idf-1 THEN Pu* 1 

IF Opo-O THEN Ya-FNPoly < Xa? 

IF Opo-2 AND I log-0 THEN Ya-Xa/FNPoly (Xa) 

IF Opo-2 AND Ilog*l THEN Ya*FNPo ly (Xa ) 

IF Opo-1 THEN Ya-FNPoly ( Xa ) 

IF Ya<Ynun THEN 6655 
IF Ya<Yl 1 OR Ya>Yul THEN 6655 
IF Iun-1 THEN 
IF 0po<2 THEN Xa-Xa*1 .3 
IF 0po>0 THEN Ya-Ya*. 1 761 
IF Opo-O THEN Ya-Ya*. 31 7 
IF Opo-2 THEN Xa-Xa*.317 
END IF 

Y-LGT(Ya/Ynun)*Sfy 

X-LGT ( Xa/Xiru n ) *Sf x 

IF Y<0 THEN Y-0 

IF Y>1 00 THEN GOTO 6655 

IF Pu-0 THEN PRINT “PA* . X . Y . “PD" 

IF Pu— 1 THEN PRINT *PA“ ,X, Y , "PU" 

NEXT Xx 
PRINT "PU" 

GOTO 6067 
END IF 
BEEP 

INPUT “WANT TO PLOT REILLY'S SM00TH-TU8E DATA? ( 1 -Y , 0-N) " , Irly 

IF Opo-O OR Opo-1 THEN 

XI 1-3 

Xul-20 

END IF 

IF Opo-2 THEN 
XI 1-10000 
Xul-100000 
END IF 
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6637 

6700 

6703 

6706 

6709 

67 1 2 

6715 

6713 

6721 

6724 

6727 

6730 

6733 

6736 

6739 

6742 

6745 

6748 

6751 

6754 

6757 

6760 

6763 

6766 

6769 

6772 

6775 

6778 

6781 

6784 

6787 

6790 

6793 

6736 

6799 

6802 

6305 

6303 

681 1 

6314 

6817 

6320 

6323 

6826 

6329 

6332 

6335 

6333 

6841 

6344 

6347 

6350 

6353 

6856 

6359 

6362 

6365 

6863 

6371 

6374 



IF IrLy»1 THEM 
Y 11*20 
Yul-70 
BEEP 

INPUT -ENTER LOUER AND UPPER Y-LIMITS FOR PLOTTING " , Y 1 1, Yul 
FOR Xx-0 TO Cx STEP Cx/200 
Xa-Xmi n-* 1 0 * Xx 

IF Xa<X 1 1 OR Xa>Xul THEN 6757 
Xl-LOG(Xa) 

IF Opo-Q THEN Yl— 2 • SA02837E- T+X 1* < 4 . 9720 1 5 1 -X I *2 . S 1 34787E- T ) 

IF Opo-1 THEN Yl— 2. 5402337E-1 +X1*( 3.97201 51 -XL*2.S134787E-1 ) 

IF Opo- 2 THEN YI— 3.7Q73801E-1+Xl*(8.7259190E-1-XI*6.3826842E-3> 
Ya-EXP(Yl) 

Y-LGT(Ya/Ynun)*Sfy 

X- LGT ( Xa/Xmin) *Sf x 
Ipu-0 

IF Y < Y 1 1 THEN Ipu-1 

IF Y> Yul THEN GOTO 6757 

IF Ipu-0 THEN PRINT "PA" , X, Y* "PD“ 

IF Ipu-1 THEN PRINT "PA" Y ♦ “PU“ 

NEXT Xx 
PRINT "PU" 

END IF 
BEEP 

INPUT "WANT TO PLOT RQHSENQW CORRELATION? ( 1 -Y , 0-N> \ Xrohs 

IF Irohs-1 THEN 

Yll-15 

Yul -30 

BEEP 

INPUT "ENTER Tsat (Deg CJ’.Tsat 

Csf - . 0040 

BEEP 

INPUT "ENTER Csf (DEF-0. 004) " .Csf 
Tf- Tsa t+2 

FOR Xx-0 TO Cx STEP Cx/200 
Xa-Xm i rv*1 0* Xx 

IF Xa<Xll OR Xa>Xul THEN 6389 

XI- LQG(Xa) 

IF 0po<2 THEN Tf- Tsa t+Xa/2 
Rho-FNRho ( Tf ) 

K-FNK( Tf ) 
rtu-FNMu ( Tf ) 

Cp-FNCp( T f ) 

Hfg-FNHf g( Tsat ) 

N i -Mu/Rho 
Pr-Cp*Mu/K 

Omega-Csf*Hfg/Cp*( ( .014/(9.31 *Rho) ) ‘ .5/(Mu*Hfg) ) * ( 1 . /3 ) *Pr ‘ 1 . 7 

IF Opo-Q THEN Ya- ( Xa/Omega ) ‘ 3 

IF Opo-1 THEM Ya- ( Xa/Omega ) *3/Xa 

IF Qpo-2 THEN Ya-Xa 4 (2 . /3 ) /Omega 

IF Qpo-2 THEN 

T f c-Tsat+Xa/Ya* .5 

IF ABS( Tf-T fc) > .01 THEN 

Tf- ( T f +Tf c) * .5 

GOTO 6814 

END IF 

END IF 

Y-LGT (Ya/Ymin)+Sfy 
X-LGT (Xa/Xmin )*Sf x 
Ipu-0 
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7060 PRINT 

7063 PRINT USING “ T 0 K . " Fms Tin Tev Vw'2 

7066 IF K-0 THEN 
7063 PRINTER IS 701 

7072 PRINT 

7075 PRINT * Month . date and time: ":Date$ 

7073 IF Itt-0 THEN PRINT USING *10X.“Tube Type: 

7081 IF Itt-I THEN PRINT USING "1 OX. “‘Tube Type: 

7084 IF Itt-2 THEN PRINT USING "10X,“Tube Type: 

7087 PRINT 

7090 PRINT USING ‘10X,** Fms Tin Tev Vu - 2 

7093 PRINTER IS 1 
7096 K-! 

7099 END IF 
7102 BEEP 

71 OS INPUT ‘ENTER FLOUMETER READING*. Fms 
7108 OUTPUT 709; ‘AR AFO AL4 VR1 “ 



Tdrop'“" 



Wieland Smooth 

High Flux " “ 
Turbo-B'”" 



Tdrop"” 



7111 FOR L-0 TO 4 
7114 OUTPUT 709: “AS SA‘ 

7117 IF L>0 AND L<4 THEN 7141 
7120 S-0 

7123 FOR 1-0 TO 9 
7126 ENTER 709:E 
7129 S-S+E 
7132 NEXT I 

7135 IF L-0 THEN Emf(0)-ABS(S/10) 

7138 IF L-4 THEN Emt ( 1 )-ABS(S/1 0) 

7141 NEXT L 

7144 OUTPUT 709: ‘AR AF20 AL20 VR1* 

7147 OUTPUT 709: ‘AS SA‘ 

7150 Etp-0 

7153 FOR 1-0 TO 9 

7156 ENTER 709:Et 

71S9 Etp-Etp+£t 

7162 NEXT I 

7165 Etp-Etp/10 

7168 T m-FNTvsv( Emt ( 1 ) ) 

7171 Tev-FNTvsv ( Emt ( 0 ) ) 

7174 Grad-37 . 9853+. 1 04388* T in 

7177 Mdot-3.96S7E-3+Fms*(3.61955E-3-Fms*<8.S2006E-6-Fms*< 1 .23638E-7-Fnw4 . 3 1 897 

E— 10))) 

Vw-Mdot/ ( 1 000+PI*D i ‘ 2) *4 
Tdrop-£tp*1 .E+6/ t 1 0*Grad) 

PRINT USING ‘1 0X.3( DO .DD.4X) . 1X.Z.DD.4X.MZ. 40" ;Fms . Tin, Tev, Vw‘2, Tdrop 
BEEP 

INPUT "WANT TO ACCEPT THIS DATA SET? ( 1-Y.0-N) ‘,0k 
J-J + 1 

IF Ok-O THEN 
J-J-1 
GOTO 7102 
ELSE 

OUTPUT SFileiFms.Emt (*) ,Etp 
PRINTER IS 701 

PRINT USING "10X.3(DD.0D.4X) .IX, Z.DD.4X.MZ. 4D‘:Fms. Tin. Tev, Vu‘2. Tdrop 
PRINTER IS 1 
BEEP 

INPUT ‘WILL THERE 8E ANOTHER DATA SET? ( 1 -Y.O-N) ‘ ,Go_on 
IF Go on- 1 THEN 7102 
END IF 

ASSIGN AFile TO * 



7180 

7183 

7186 

7139 

7192 

7195 

7198 

7201 

7204 

7207 

7210 

7213 

7216 

7219 

7222 

7225 

7228 

7231 

7234 
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7.237 PRINTER IS 701 
7240 PRINT 

7243 PRINT USING * 1 0X . * "NOTE : "\ZZ."’ data sets are storea in file 
i lei 

7246 PRINTER IS 1 
7249 SUBEND 
7252 SUB Uoort 
7255 PRINTER IS 1 
7258 BEEP 

7261 INPUT “Enter Uo File Name"*FileS 
7264 BEEP 

7267 INPUT "Number of Data Runs" ,Nrun 

7270 INPUT “Do You Want a PLot File?<1 -Y ,0-N) tt t Iplot 

7273 BEEP 

7276 IF Iplot-1 THEN 

7279 INPUT “Give Plot File Name 0 ,P_f i leS 

7282 CREATE BDAT P_fileS.4 

7285 ASSIGN 8Plot TO P files 

7288 END IF 

7291 PRINTER IS 701 

7294 PRINT 

7297 PRINT 

7300 PRINT USING "1GX."" Water Vei Uo""" 

7303 ASSIGN SFiie TO FiieS 

7306 IF Iplot-1 THEN ASSIGN SFilef TO P file* 

7309 FOR 1-1 TO Nrun 

7312 ENTER *File:Uu*,Uo 

7315 IF Iplot-1 THEN OUTPUT SFilel iVw.Uo 

7318 PRINT USING "1SX.D.DD.SX.MZ.3DE" : Vw.Uo 

7321 NEXT I 

7324 ASSIGN SFiie TO > 

7327 ASSIGN *File1 TO * 

7330 PRINT USING "10X,""N0TE: "",ZZ,"" data sets are stored in file 
n. FiieS 



7333 IF Iplot-1 THEN 

7336 PRINT USING "1 OX. ""NOTE: "".ZZ*"" X-Y Pairs are stored in file 
n, P_f ileS 



7339 

7342 

7345 

7348 

7351 

7354 

7357 

7360 

7363 

7366 

7369 

7372 

7375 

7378 

7381 

7384 

7387 

7390 

7393 

7396 

7399 

7402 

7405 



END IF 

PRINTER IS 1 
SUBEND 
SUB Select 



COM /Idp/ Idp 
BEEP 

PRINTER IS 1 

PRINT USING "4X. “"Select option:""" 



PRINT USING " 6X . " " 0 
PRINT USING " 6X , " " 1 
PRINT USING "6X . "" 2 
PRINT USING °6X* " " 3 
PRINT USING °6X , " " 4 
PRINT USING *6X f " " 5 
PRINT USING "6X . " " 6 
PRINT USING "SX • " " 7 
PRINT USING "SX,"" 8 
PRINT USING "SX, * " 9 
PRINT USING "SX .""10 
PRINT USING "6X.""11 
INPUT Idp 

IF Idp-0 THEN CALL Main 
IF Idp-1 THEN CALL Plot 



Tatcing data or re-processing old data" " M 
Plotting data on Log-Log """ 

Plotting data on Linear""" 

Malce cross-plot coefft file""" 
Re-circulate water""" 

Purge " " “ 

T-Drop correction""" 

Print Uo File""" 

Modify X-Y file""" 

Move""" 

Comb""" 

Fixup""" 



" " , 1 5 A " ; J , F 



" " , ISA" ;Nru 
"" * ISA" ;Nru 
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7408 IF Idp-2 THEN CALL PI in 
7411 IF Idp-3 THEN CALL Coef 
7414 IF Idp-4 THEN CALL Main 
7417 IF Idp-S THEN CALL Purg 
7420 IF I dp-6 THEN CALL Tdcn 
7423 IF I dp-7 THEN CALL Uoprt 
7426 IF Idp-8 THEN CALL Xymod 
7429 IF Idp-9 THEN CALL Move 
7432 IF Idp-10 THEN CALL Comb 
7435 IF Idp-11 THEN CALL Fixup 
7438 SUBEND 
7441 SUB Xymod 
7444 PRINTER IS 1 
7447 BEEP 

7450 INPUT “ENTER FILE NAME", Fi lei 
7453 ASSIGN GFilel TO FileS 
7456 BEEP 

7459 INPUT “ENTER NUMBER OF X-Y PAIRS", Np 
7462 BEEP 

7465 INPUT “ENTER NEH FILE NAME\FiIe2S 
7468 CREATE BDAT File2S.S 
7471 ASSIGN $File2 TO File2S 
7474 BEEP 

7477 INPUT “ENTER NUMBER OF X-Y PAIRS TO BE DELETED* .Ndel 
7480 IF Ndel-0 THEN 7492 
7483 FOR 1-1 TO Ndel 
7486 BEEP 

7489 INPUT “ENTER DATA SET NUMBER TO BE DELETED* ,Nd( I ) 

7492 NEXT I 

7495 FOR J-l TO Np 

7498 ENTER $Fi lei :X.Y 

7501 FOR 1-1 TO Ndel 

75d4 IF Nd( I ) - J THEN 7516 

7507 NEXT I 

7510 OUTPUT «File2:X.Y 

7513 PRINT J.X,Y 

7516 NEXT J 

7519 PRINTER IS 701 

7522 ASSIGN AFilel TO * 

7525 ASSIGN $File2 TO * 

7528 SUBEND 
7531 SUB Move 
7534! FILE NAME: MOVE 
7537! 

7540 DIM Bop (66) ,A<66) ,B(66) ,C(66) ,D(66) ,E(66) ,F(66) ,G(66) ,H(66) , J(66) ,K<66) ,L( 
66) , M(66) 

7543 DIM Told$(66) [ 14] ,N(66) ,Vr(66) ,Ir<66) 

7546 BEEP 

7549 INPUT “OLD FILE TO MOVE" ,D2_f i leS 
7552 ASSIGN «File2 TO D2_file$ 

7555 ENTER *File2;Nrun, Dates. Ldtcl ,Ldtc2, Itt 

7558 FOR 1-1 TO Nrun 

7561 ENTER *Fi le2;Bop( I) .To ld$< I) 

7564 ENTER JFile2:A( I) ,3(1) ,C( I ) ,D( I) ,£( I) ,F( I) ,G( I) ,H< I ) , J( I) .K( I) .L( I ) ,M< I ) .N 
(I) 

7567 ENTER »File2:Vr( I) ,Ir(I) 

7570 NEXT I 

7573 ASSIGN ?File2 TO * 

7576 BEEP 

7579 INPUT 'SHIFT DISK AND HIT CONTINUE" ,0k 
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7532 BEEP 

7535 INPUT “INPUT BDAT SIZE", Size 

7533 CREATE BDAT D2_f i leS.Si ze 

7591 ASSIGN SFilel TO D2_fileS 

7594 OUTPUT SFilel :Nrun, DateS, Ldtcl ,Ldtc2. Itt 

7597 FOR I- I TO Nrun 

7G00 OUTPUT OFilel :Bop(I) .ToIdSCI) 

7303 OUTPUT SFiiel :A( I) ,B( I) ,C( I) .D( I) .£( I) ,F( I) ,G( I) ,H( I) . J< I) ,K( I) ,L( I) ,f1< I) . 
N ( I) 

7606 OUTPUT SFilel ;Vr( I) , Ir< I) 

7609 NEXT I 

76 1 2 ASSIGN SFiiel TO * 

7615! RENAME "TEST* TO D2_fileS 

7613 BEEP 2000. .2 
7621 BEEP 4000.. 2 
7624 BEEP 4000.. 2 

7627 PRINT "DATA FILE MOVED* 

7630 SUBEND 
7633 SUB Comb 

7636! FILE NAME: COMB 
7639! 

7642 DIM Emt (12) 

7645 BEEP 

7648 INPUT ‘OLD FILE TO FIXUP*. D2 fiieS 
7651 ASSIGN 8File2 TO D2 fiieS 
7654 D1_f ileS“*TEST’ 

7657 CREATE BDAT D1_fileS.30 

7660 ASSIGN SFiiel TO D1_fileS 

7663 ENTER SFile2:Nrun. DateS. Ldtcl ,Ldtc2.Itt 

7666 IF K-0 THEN OUTPUT *F i lei iNrun, DateS.Ldtc I ,Ldtc2 , Itt 

7669 FOR I-l TO Nrun 

7672 ENTER SFile2:3op.ToldS.Erat (*> .Vr.Ir 
7675 OUTPUT SFilel :8op . ToldS.Emt < *) , Vr . Ir 
7678 NEXT I 

7631 ASSIGN SFile2 TO * 

7634! RENAME 'TEST” TO D2_fileS 
7637 BEEP 2000.. 2 
7690 BEEP 4000.. 2 
7693 BEEP 4000,. 2 
7696 BEEP 

7699 INPUT "UANT TO ADD ANOTHER FILE < 1 -Y . 0-N) T . Oka 

7702 IF Oka- 1 THEN 
7705 K-l 

7703 BEEP 

7711 INPUT "GIVE NEW FILE NAME" ,Nf i leS 
7714 ASSIGN SFile2 TO NfileS 
7717 GOTO 7663 
7720 END IF 

7723 ASSIGN ?Fiie2 TO * 

7726 SUBEND 

7729 SUB Fixup 

7732! FILE: FIXUP 

7735! DATE: February 18, 1986 

7738 ! 

7741 DIM Emt (12) 

7744 BEEP 

7747 INPUT "OLD FILE TO FIXUP’ .D2_fi le$ 

7750 ASSIGN SFiIe2 TO D2_FUeS 
7753 D1 Files- “TEST" 

7756 CREATE BDAT D1_fiieS,20 
7759 ASSIGN SFiiel TO D1_FileS 
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7762 ENTER '«Fi.le2lNrun, Dates. Ldtcl .LdtcZ.Itt 
7765 Nruno-27 

7763 OUTPUT $Filet :Nruno. Datei. Ldtcl ,Ldtc2. Itt 
7771 FOR 1-1 TO Nrun 

7774 ENTER 9Fil«2:Bop,ToldS,Enf (*) .Vr.Ir 
7777 IF 1-1 THEN 7783 

7780 OUTPUT SFilel :8op,ToldS,Emt (*) .Vr.Ir 
7783 NEXT I 

7786 ASSIGN 3File2 TO * 

7789 ASSIGN SFilel TO * 

7792! RENAME "TEST" TO D2_file$ 

7795 BEEP 2000.. 2 
7798 BEEP 4000.. 2 
7801 BEEP 4000,. 2 
7804 SUBEND 
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APPENDIX C 



EXAMPLES OF REPRESENTATIVE DATA RUNS 
The two following data runs (TXE222 and GTB185) are 
representative samples of the data taken over the course of 
the investigation. The first data run for the Thermoexce 1 -E 
surface was for a decreasing heat-flux condition and an oil 
concentration of three percent. The second data run was for 
the GEWA-T surface was also for a decreasing heat-flux 
condition but with an oil concentration of one percent. 



Month, date and time :05:0S: II :41 : 20 

NOTE: Program nan* : DRP4 
Olsk numoer • 13 
Old file name: TXE222 



Th 


13 data set taken on i 02: T8: 


13:42:58 




Tube Nuraoer: TO 








Data Set 


Number - t 3ulk Oil V, • 


3.0 


02: 


13:13:48:44 


TC No: 


12 3 4 


S 


6 


7 3 


Temp : 


8.7! 9. 45 9.00 3.48 


9.70 


3.12 


10.21 3.59 


Tea 


Tliqd Tliqd2 Tvaor Psat 


Tsump 






3.74 


2.23 2. T 7 2.06 -1 .27 


-18.9 






Thetab 


Htuoe Qdp 








6.54 1 


! . 45SE+04 9.5l9E*04 






- 


Data Set 


Numoer * 2 3u i k Oil /C - 


3.0 


02: 


13: 13:49:57 


TC No: 


12 3 4 


5 


6 


7 3 


Temp : 


8. SI 9.01 8.42 3.00 


9.68 


3.29 


10.33 8.51 


Tea 


Tliqd Tliqd2 Tvaor Psat 


Tsurao 






8.S5 


2.24 2.22 2.09 -1.73 


-18.3 






Thetab 


Htuoe Qdp 








S.3T9 


1 .S08E-04 9.S3TE+04 








Data Set 


Numoer • 3 3uik Oil 2 - 


3.0 


02: 


18: 13:50:51 


TC No: 


12 3 4 


5 


6 


7 8 


Temp : 


8.68 3.38 3.33 7.39 


9.34 


8.07 


10.22 8.64 


Tea 


Tli od Tliqd2 Tvapr Psat 


Tsump 






8. S3 


2.26 2.19 2.07 -1.74 


-18.8 






Thetao 


Htuoe Qdp 








S.30S 


I.StlE+04 9.S26E+Q4 








Data Set 


Nuraoer - 4 8uilc Oil 7. - 


3.0 


02: 


18:14:02:37 


TC No: 


12 3 4 


5 


6 


7 8 


Temo : 


6.93 7.05 6.71 6.48 


7.23 


6.24 


7.63 6.56 


Twa 


Tliqd Tliqd2 Tvapr Psat 


Tsumo 






S . oto 


2.25 2.21 2.14 -1.73 


-18. 5 






Thetab 


Htuoe Qdp 








4.433 


1.414E+04 S.2S9E+04 








Data Set 


NumOer - 5 Bulk Oil 7. - 


3.0 


02: 


18: 14:03: 13 


TC No: 


12 3 4 


5 


6 


7 8 


Temo : 


6.32 7.18 6.30 6.56 


7.35 


6.25 


7.65 6.52 


Tea 


Tliqd Tliqd2 Tvapr Psat 


Tsurao 






S. 7 1 


2.23 2.25 2.18 -1.69 


-18.4 






Thetao 


Htuoe Qdp 








4.444 


1 . 408E+04 6.2S8E-04 








Data Set 


Nuraoer - 6 3uik Oil 2 • 


3.0 


02: 


18:14:13:56 


TC No: 


12 3 4 


5 


6 


7 8 


Temo : 


S.S0 5.65 5.41 5.26 


5.58 


4.90 


5.84 5.06 


Tea 


Tliqd Tliqd2 Tvaor Psat 


Tsump 






S. 28 


2.32 2.28 2.35 -1.66 


-16.3 






Thetao 


Htuoe Qdp 








2.98 T 


1 . 304E+04 3.388E+04 
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Oata Set 


Numoer 


• 7 Bulk Oil 7. - 


3.0 


02: 


13:14:14:29 


TC No: 


1 


2 3 4' 


5 


6 


7 B 


Temo : 


5.45 


5.60 S.37 S. 3 1 


5.61 


4.38 


5.32 5.02 


Twa 


T 1 iqd 


Tliqd2 Tvapr Psat 


T sumo 






5.26 


2.29 


2.26 2.35 -1.68 


-16.2 






Thetab 


Htube 


QdP 








2.387 


1.299E* 


■04 3.3316*04 








Data Set 


Numoer 


• 8 Bulk Oil Z - 


3.0 


02: 


13:14:20:33 


TC No: 


1 


2 3 4 


5 


6 


7 B 


Temp : 


4.25 


4.33 4.22 4.16 


4.49 


4.07 


4.58 4.06 


Twa 


T1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






4.20 


2.24 


2.16 2.64 -1.77 


-15.3 






Thetab 


Htube 


Qdp 








1.998 


1 . 1 76E*04 2.3496*04 








Oata Set 


Number 


- 9 Bulk Oil /. • 


3.0 


02: 


18:14:21 :05 


TC No: 


1 


2 3 4 


5 


6 


7 3 


Temp : 


4.33 


4.41 4.26 4.23 


4.52 


4.10 


4.65 4.11 


Twa 


T 1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






4.25 


2.29 


2.23 2.71 -1.71 


-15.3 






Thetab 


Htuoe 


Qdp 








2.004 


1 . 1696*04 2.3446*04 








Oata Set 


Numoer 


• 10 Bulk Oil 7. - 


3.0 


02: 


18:14:26:17 


TC No: 


1 


2 3 <4 


5 


6 


7 B 


Temp : 


3.63 


3.58 3.59 3.54 


3.30 


3.51 


3.35 3.47 


T ua 


T1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






3.59 


2.33 


2.27 2.95 -1.65 


-16.4 






Thetao 


Htube 


Qdp 








1.293 


1.1376*04 1 . 4696+04 








Oata Set 


Numoer 


• 11 Bulk Oil Z - 


3.0 


02: 


18:14:26:48 


TC No: 


1 


2 3 4 


5 


6 


7 B 


T emp : 


3.56 


3.61 3.55 3.51 


3.75 


3.44 


3.80 3.39 


Twa 


T 1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






3.53 


2.28 


2.21 2.34 -1.72 


-16.4 






Thetab 


Htube 


Qdp 








1.285 


1 . 1 426*04 1.4676*04 








Oata Set 


Number 


• 12 Bulk Oil 7. - 


3.0 


02: 


18:14:33:1 1 


TC No: 


1 


2 3 4 


5 


6 


7 3 


Temp : 


3.06 


3.07 3.02 3.00 


3.12 


2.98 


3.14 2.90 


Tua 


T1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






3.01 


2.43 


2.11 3.43 -1.69 


-17.1 






Thetab 


Htuoe 


Qdp 








.739 


1 . 1206 


*04 3.283E*03 


• 






Oata Set 


: Numoer 


- 13 Bulk Oil 7. - 


3.0 


02: 


18:14:33:43 


TC No: 


1 


2 3 4 


5 


6 


7 B 


Temp : 


3.00 


3.01 2.94 2.92 


3.05 


2.90 


3.08 2.81 


Twa 


T 1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






2.94 


2.36 


2.05 3.41 -1.77 


-17.2 






Thetab 


Htuoe 


Qdp 








*’ .738 


1 . 1366 


*04 3.386E*03 








Oata Set 


: Numoer 


- 14 Bulk Oil Z - 


3.0 


02:18:14:40:43 


TC No: 


1 


2 3 4 


5 


6 


7 3 


Temp : 


2.70 


2.74 2.73 2.71 


2.30 


2.71 


2.79 2.61 


Tua 


T 1 iqd 


Tliqd2 Tvapr Psat 


Tsump 






2.71 


2.41 


2.03 3.79 -1.75 


-17.9 






Thetab 


Htube 


Qdp 








.488 


1.0786*04 5.2S7E+03 









Oata Set Number - IS Bulk Oil Z • 3.0 02:18:14:41:19 

TC No: 1 2 3 4 S S 7 3 

Temp : 2. 66 2.70 2. 65 2.67 2.74 2.66 2.71 2.54 

Twa Tliod Tliqd2 Tvapr Psat Tsump 

2.65 2.39 2.00 3.78 -1.77 -18.0 

Thetab Htube Qdp 

.453 1 . 1 45E*Q4 5.2446*03 
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Month, date and time :05 : 06 : 1 l : 47 : 44 



MOTE: Program name : DRP4 
Oisk numoer • 10 
Old file name: GTB18S 



This data 


set taken on *. 11:25: 


11:35:33 


Tuoe Mumoer: 9 






Data Set 


Mumoer 


1 Bulk Oil 7. - 


1.0 


11:25:15:46:12 


TC Mo: 


1 


2 3 4 


S 


6 7 8 


Temo : 


3.27 


3.61 7.33 7.87 


7.93 


7.32 9.85 7.81 


Tuia 


T l iqd 


Tliqd2 Tvaor Psat 


T3umo 




7.97 


2.23 


2.23 2.24 -1.73 


-16.6 




Thetao 


Htuoe 


Odo 






5 . 739 


1.581E+04 9 .073004 






Data Set 


Mumoer 


- 2 Bulk Oil 7. - 


1.0 


1 1 : 25 : 1 S : 46 : 48 


TC Mo: 


1 


2 3 4 


S 


6 7 3 


Temo : 


8.20 


8.65 7.86 7.90 


7.38 


7.85 9.88 7.34 


Twa 


Tl iqd 


Tliqd2 Tvaor Psat 


Tsumo 




8.00 


2.28 


2.27 2.28 -1.68 


-16.6 




Thetao 


Htuoe 


Qdo 






5 . 728 


1.582E< 


►04 9.062E+04 






Data Set 


Mumoer 


- 3 Bulk Oil 1 - 


1.0 


11:25:15:54:05 


TC Mo: 


1 


2 3 4 


5 


6 7 8 


Temo : 


6.93 


7.13 6.62 6.52 


6.64 


6.50 7.78 6. 54 


Tuia 


Tl iqd 


Jliqd2 Tvaor Psat 


Tsumo 




6.66 


2.21 


2.22 2.77 -1.75 


-17.0 




Thetao 


Htuoe 


Odo 






4.448 


1.348E- 


►04 S.995E-04 






Data Set 


Mumoer 


• 4 Bulk Oil 7. - 


1 .0 


11:25:15:54:55 


TC Mo: 


1 


2 3 4 


S 


6 7 8 


Temo : 


6.9S 


7. IS 6.61 6.62 


6.62 


6.49 7.77 6.50 


Twa 


Tl iqd 


Tliqd2 Tvaor Psat 


Tsumo 




6.66 


2.22 


2.23 2.81 -1.74 


-16.3 




Thetao 


Htuoe 


Odo 






4.430 


1.3S3E+04 5.993E+04 






Data Set 


Numoer 


- S Bulk Oil 7 - 


1.0 


11:25: 16:03:48 


TC No: 


1 


2 3 4 


S 


6 7 8 


Temo : 


6.09 


6.19 S.32 S.31 


5.77 


S . 74 6.37 S.66 


Twa 


Tl iqd 


Tliqd2 Tvaor P3at 


Tsumo 




5.32 


2.24 


2.23 3.18 -1.73 


-17.4 




Thetao 


HtuDe 


Odo 






3. 533 


1.036E 


■►04 3.712E-04 






Data Set 


Mumoer 


- 6 Bulk Oil 7. - 


J.O 


11:25:16:04:22 


TC Mo: 


1 


2 3 4 


5 


6 7 8 


Temo : 


6. 1 1 


6.23 5.34 5.34 


5.30 


S. 77 6.42 5.73 


Twa 


Tl iqd 


Tliqd2 Tvaor P3at 


Tsumo 




S.8S 


2.27 


2.26 3.22 -1.69 


-17.4 




Thetao 


HtuOe 


Odo 






3.S88 


1.034E+04 3.7T2E+04 






Data Set 


: Mumoer 


7 Bulk Oil 7. - 


1.0 


11:25:16:10:15 


TC Mo: 


1 


2 3 4 


s 


6 7 8 


Temo : 


S . 4 1 


S.S2 S.19 S.32 


S .25 


S.16 S.S0 S.12 


T wa 


Tl iqd 


Tliqd2 Tvaor Psat 


Tsumo 




5.24 


2.24 


2.20 3.61 -1.74 


-17.7 




Thetao 


HtuOe 


Odo 






3.020 


7.320E< 


'03 2.211E+04 






Data Set 


Mumoer 


8 Bulk Oil 7. • 


1.0 


11:25:16:12:11 


TC Mo: 


1 


2 3 4 


5 


6 7 8 


Temo : 


S.2S 


S.32 4.98 S. 1 3 


S.06 


4.95 S.29 4.90 


Twa 


Tl iqd 


Tliqd2 Tvaor Psat 


Tsumo 




5.04 


2.14 


2.05 3. 56 -1.38 


-I 7.7 




Thetao 


Htuoe 


Odo 






2.949 


7 . 480E+03 2.20SE*04 






Data Set 


Mumoer 


- 9 3ulk Oil 7. - 


1 .0 


11:25:16:17:23 


TC No: 


1 


2 3 4 


S 


6 7 8 


Temo : 


5. 3S 


5.14 4.81 S.OS 


4.34 


4.81 4.95 4.78 


Twa 


Tl iqd 


Tliqd2 Tvaor Psat 


Tsumo 




4.30 


2.2S 


2.22 3.45 -1.73 


-17.7 




Thetao 


Htuoe 


Odo 






2.663 


5.11 4E+Q3 1.262E+04 






Data Set 


Number 


- 10 Bulk Oil 7. - 


1.0 


11:25:16:18:04 


TC Mo: 


1 


2 3 4 


S 


6 7 8 


Temo : 


S.10 


5.20 4.83 S.09 


4.99 


4.85 S.OO 4.85 


Twa 


Tl iqd 


Tliqd2 Tvaor P3at 


Tsumo 


4.9S 


2.30 


2.27 3.44 -1.67 


-17.6 




ThetaO 


HtuOe 


Qdo 






2.6S7 


S.103E+03 1.3S6E+04 
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Data Set 


Number 


- 1 1 Bulk Oil 2 “ 1-0 


TC No: 


1 


2 3 4 5 


Temp : 


4.75 


4.38 4.42 4.34 4.75 


Twa 


Tl iqd 


TUqd2 Tvapr Psat Tsump 


4.66 


2.30 


2.29 3.15 -1.66 -16.3 


Thetab 


Htube 


Qdp 


2.360 


3.273E+03 7.72SE+03 


Data Set 


Number 


- 12 Bulk Oil 7 - 1-0 


TC No: 


1 


2 3 4 5 


T amp : 


4.76 


4.38 4.42 4.34 4.74 


Twa 


Tl iqd 


Tiiqd2 Tvapr Psat Tsump 


4.66 


2.30 


2.29 3.13 -1.66 -16.0 


Thetab 


Htuoe 


Qdp 


2.360 


3.264E+03 7.702E-03 


Data Set 


Numoer 


- 13 Bulk Oil 7 - 1.0 


TC No: 


1 


2 3 4 5 


Temp ' 


4.45 


4.54 4.09 4.60 4.39 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat Tsump 


4.37 


2.27 


2.27 2.95 -1.69 -14.9 


Thetab 


Htube 


Qdp 


2.094 


2.361E+Q3 4.944E+03 


Data Set 


: Numoer 


- 14 Bulk Oil 7 - 1.0 


TC No: 


1 


2 3 4 5 


Tamp : 


4.46 


4.55 4.12 4.61 4.43 


* Twa 


Ti iqd 


Tliqd2 Tvapr Psat Tsump 


4.39 


2.30 


2.29 2.94 -1.66 -14.8 


Thetab 


Htube 


Qdp 


2.096 


2.352E+03 4.931E+03 



11:25:16:24:54 
6 7 3 

4.61 4.62 4.53 



11:25:16:25:33 
6 7 3 

4.61 4.60 4.53 



1 1 : 25 : 16:31 :21 
6 7 3 

4.34 4.36 4.28 



11:25:16:31:59 
6 7 8 

4.38 4.38 4.32 



Data Set 


Number 


• 15 Bulk Oil V. - 


1.0 


11:25:16:37:53 


TC No: 


1 


2 3 4 


5 


6 7 B 


Tamp : 


4.19 


4.29 3.34 4.38 


4.07 


4.06 4.00 3.99 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




4.09 


2.23 


2.23 2.33 -1.73 


-14.3 




Thetab 


Htuoe 


Qdp 






1.863 


1.S82E*03 2.946E+03 






Data Set 


Number 


- 16 Bulk Oil 7 - 


1.0 


11:25:16:38:28 


TC No: 


1 


2 3 4 


5 


6 7 8 


Temp : 


4.20 


4.30 3.83 4.36 


4.07 


4.05 3.98 3.99 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




4.09 


2.23 


2.22 2.32 -1.74 


-14.3 




Thetab 


Htuoe 


Qdp 






1.362 


l.582E« 


>03 2. 946E+03 






Data Set 


Number 


- 17 Bulk Oil 7. - 


1.0 


11 :2S.M6:44:S6 


TC No: 


-1 


2 3 4 


5 


6 7 B 


Tamp : 


3.97 


4.04 3.67 4.17 


3.77 


3.31 3.67 3.74 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




3.35 


2.19 


2.19 2.71 -1.78 


-15.2 




Thetao 


Htube 


Qdp 






1.660 


1 . 1 46E- 


>03 1.901E+03 






Data Set 


Number 


- 18 Bulk Oil 7. - 


1.0 


11:25:16:45:33 


TC No: 


1 


2 3 4 


5 


6 7 3 


Temp : 


3.97 


4.04 3.68 4.19 


3.78 


3.82 3.66 3.72 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




3.85 


2.19 


2.19 2.71 -1.77 


-15.2 




Thetab 


Htuoe 


Qdp 






1 .660 


1.147E*03 1.904E+Q3 






Data Set 


Numoer 


- 19 Bulk Oil 7. - 


1.0 


11:25:16:51:19 


TC No: 


1 


2 3 4 


5 


6 7 8 


Temp : 


3.69 


3.73 3.50 3.77 


3.40 


3.46 3.25 3.32 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




3.5! 


2.19 


2.18 2.53 -1.78 


-15.4 




Thetab 


Htube 


Qdp 






1.326 


3.S73E+Q2 1.137E+03 






Data Set 


Number 


- 20 Bulk Oil 7. - 


1.0 


11:25:16:53:16 


TC No: 


1 


2 3 4 


5 


6 7 8 


Tamp : 


3.77 


3.31 3.57 3.84 


3.46 


3.54 3.32 3.38 


Twa 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




3.58 


2.25 


2.25 2.62 -1.71 


-15.7 




Thetab 


Htube 


Qdp 






1.335 


8.630E+Q2 1. IS2E+03 
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APPENDIX D 



UNCERTAINTY ANALYSIS 

The uncertainty of the heat- transfer coefficient at 60 
kW/m 2 and 8 kW/m 2 for data runs of the GEWA-T tube (GTB185 
with one percent oil concentration) and of the Thermoexcel- 
E (TXE222 with three percent oil concentration) are 
analyzed below. The method of analysis is based on Kline- 
McClintock [Ref. 20] method of uncertainty analysis. 

The heat- transfer coefficient is: 

h = (q c )/(T w0 - Ts a t ) (D.l) 



and 




Two — T g a T = Twt 


-[Qc In (D 7 /D,)3/<2 7T k L) - T s A T (D.2) 



where 



h = 


heat-transfer coefficient 


qc = 


heat flux corrected for end loss 


TuO = 


average outer tube wall temperature 


T a a t = 


saturation temperature 


T„ , = 


average inner tube wall temperature 


Qc = 


heat input corrected for end losses 


D, = 


tube internal diameter 


d 2 = 


tube external diameter 


k = 


tube wall thermal conductivity 


L = 


length of tube’s heated surface 
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To begin the analysis, a dummy variable is assigned to the 
conduction term as follows: 

DELT = C Q c In <D?/D,)]/< 2 7T k L) (D.3) 

In accordance with Kline and McClintock, the 

uncertainty of the heat-transfer coefficient is: 




The uncertainty for the conduction term is estimated since 
the logarithmic error is negligible when compared to the 
other error of the analysis. The uncertainty of the 
conduction term is as follows: 



<5delt 

DELT 




where Qc = qc 7T D 2 L 

The uncertainty of this term is as follows: 



(D. 6) 



<5q c = 

Qc 





I / 2 

(D. 7) 



Table 4 contains a listing of the computed terms of 
equations D.i through D.7 for selected data runs for the 
GEWA-T and the Thermoe xce 1 -E tubes at specified heat-flux 
setting of 60 kW/m 2 and 8 kW/m 2 . The algorithm in the data- 



141 



TABLE 4 

UNCERTAINTY ANALYSIS 

File Name GTB185 TXE222 



Heat Flux 


60 kW/m2 


8 kW/m2 


60 kW/m2 


8 kW/m2 


Oi 1 % 


1% 


1% 


3% 


3 % 


«J 3.g_ 

Qc 


0.015 


0.015 


0.015 


0. 015 


6 _ k 










k 


0.03 


0.03 


0.03 


0. 03 


<5jL 










L 


0.0025 


0.0025 


0. 0025 


0.0025 


DELT 


0.486 


0. 064 


0. 486 


0.064 


5DELT 










DELT 


0.069 


0.520 


0. 069 


0.526 


T u , ( C) 


6.845 


4. 678 


6.904 


2. 964 


Tut ( C) 


0. 49 


0.065 


0.49 


0.065 


T s a T ( C) 


2. 215 


2. 295 


2. 265 


2. 205 


T sat (K) 


0. 1 


0. 1 


0. 1 


0. 1 


5 Tu i 

T« o “Tg ft t 


0.111 


0.028 


0.111 


0. 088 


5 DELT 

Two “T sat 


0. 008 


0.014 


0.008 


0.046 


_ 6Tsat 

Two “Tsat 


0.023 


0. 042 


0. 023 


0. 135 


q c 


0. 02 


0. 02 


0.02 


0.02 


oh 










h 


0. 115 


0.056 


0. 115 


0. 169 


h ( W/ m 2 K) 


13480 


3273 


14140 


11200 
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reduct i on 



program used this approach and assumed that all 



the eight wall thermocouples were used to calculate 
heat- transfer coefficient. 



the 
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